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Introduction 

O. Behnke, M Cacciari, M. Corradi, A. Dainese, A. Geiser, A. Meyer, M. Smizanska, U. Uwer, and C. Weiser 



Perturbative QCD is expected to provide reliable predictions for the production of bottom and (to a lesser ex- 
tent) charm quarks since their masses are large enough to assure the applicability of perturbative calculations. 
Anyway a direct comparison of perturbative QCD predictions to heavy flavour production data is not straight- 
forward. Difficulties arise from the presence of scales very different from the quark masses that reduce the 
predictivity of fixed-order theory, from the non-perturbative ingredients needed to parametrize the fragmenta- 
tion of the heavy quarks into the observed heavy hadrons and from the limited phase space accessible to present 
detectors. Moreover a breakdown of the standard collinear factorization approach can be expected at low-x. 
The study of heavy quark production in hadronic interactions and in ep collisions at HERA has been therefore 
an active field, in the effort to overcome these difficulties and to get a deeper understanding of hard interactions. 

Besides its intrinsic interest, a precise understanding of heavy quark production is important at LHC 
because charm and beauty from QCD processes are relevant backgrounds to other interesting processes from 
the Standard Model (e.g. Higgs to bb) or beyond. Moreover, theoretical and experimental techniques developed 
at HERA in the heavy quark field, such as heavy-quark parton densities or 6-tagging, are also of great value for 
future measurements at the LHC. 

The present status of the heavy-quark production theory is critically reviewed in the first contribution. 
The second contribution summarises the present heavy flavour data from HERA and gives an outlook of what 
can be expected from HERA-II. The potential of the LHC experiments for charm and beauty physics is reviewed 
in the 3rd contribution. Then the relevance of saturation and low-x effects to heavy quark production at HERA 
and at the LHC are discussed. The non-perturbative aspects of heavy-quark fragmentation and their relevance to 
HERA and LHC are discussed in the next contribution. Finally, a comparison of different theoretical predictions 
for HERA and LHC based on different approaches is presented. 
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Abstract 

We review some of the main theoretical aspects of heavy quark production at HERA 
that will be important for understanding similar processes at the LHC. 

coordinators: M. Cacciari, and E. Laenen 

1 Introduction 

The value for the LHC physics program of heavy quark production studies at HERA consists not only of 
measured quantities such as parton distributions, heavy quark masses etc. but at least as much of the theoretical 
ideas on heavy quark production that were developed and refined in the course of these studies. The strong 
experimental interest in heavy quark observables at HERA has led to a significantly increased understanding of 
the benefits and limitations of finite order calculations. It has stimulated theorists to deepen their insight into 
the issue of when a heavy quark should be treated as a parton, and it has provoked novel proposals to explain 
the hadronization of heavy quarks. In what follows we review and critically assess some of these ideas. 

2 Heavy quark production 

The study of heavy quarks, historically plagued by low production rates and large uncertainties, has now entered 
the regime of 'precision physics'. On the one hand, the larger centre-of-mass energies of the colliders running 
now (Tevatron, HERA) and in the near future (LHC) lead to a much more copious production yield. On the 
other hand, technological advances such as the introduction of microvertex detectors based on semiconductor 
devices allow for much better tagging of the produced heavy flavours, and hence better measurements. Needless 
to say, an equally substantial improvement of the theoretical calculations has been needed in order to match 
this progress and therefore deliver predictions with an accuracy at least as good as that of the experimental 
measurements. Properly testing and constraining the theoretical calculations will in turn help in refining the 
predictions for the LHC. 

One example for which a good theoretical accuracy at the LHC is desirable is in calculating the total Z 
boson production rate, a process which can be used as a luminosity candle and which we would hke to have 
under control at the one per cent level. One chaimel contributing to this process is gluon-gluon fusion followed 
by bottom-antibottom annihilation, gg ^ bb ^ Z. This chaimel provides about 5% of the total Z yield [1]: 
hence, it must be under control at the 20% level in order to achieve the sought-for final 1% accuracy. 
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As it turns out, it is more efficient and more reliable to rewrite this in terms of a perturbatively calculated 
parton distribution function (PDF) for the bottom quark, i.e. as the effective process bb Z. The theoreti- 
cal tools that we use to construct such heavy quark parton distribution functions must therefore be tested by 
employing them in other theoretical predictions, to be compared to the available experimental data. In the 
following section we shall list a number of examples where this is done. 

From the point of view of 'standard' perturbative QCD calculations, the situation has not changed since 
the beginning of the '90s: fully massive next-to-leading order (NLO) calculations were made available for 
hadron-hadron [2, 3, 4, 5, 6], photon-hadron [7, 8, 9] (i.e. photoproduction) and electron-hadron [10, 11, 12, 13] 
(i.e. Deep Inelastic Scattering, DIS) collisions. These calculations still constitute the state of the art as far as 
fixed order results are concerned, and they form the basis for all modem phenomenological predictions. 

Over the years, and with increasing experimental accuracies, it however became evident that perturbative 

QCD alone did not suffice. In fact, real particles - hadrons and leptons - are observed in the detectors, not 
the unphysical quarks and gluons of perturbative QCD. A proper comparison between theory and experiment 
requires that this gap be bridged by a description of the transition. Of course, the accuracy of such a description 
will reflect on the overall accuracy of the comparison. When the precision requirements were not too tight, 
one usually employed a Monte Carlo description to 'correct' the data, deconvoluting hadronization effects and 
extrapolating to the full phase space. The final 'experimental' result could then easily be compared to the 
perturbative calculation. This procedure has the inherent drawback of including the bias of our theoretical 
understanding (as implemented in the Monte Carlo) into an experimental measurement. This bias is of course 
likely to be more important when the correction to be performed is very large. It can sometimes become almost 
unacceptable, for instance when exclusive measurements are extrapolated by a factor of ten or so in order to 
produce an experimental result for a total photoproduction cross section or a heavy quark structure function. 

The alternative approach is to present (multi)differential experimental measurements, with cuts as close 
as possible to the real ones, which is to say with as little theoretical correction/extrapolation as possible. The 
theoretical prediction must then be refined in order to compare with the real data that it must describe. This 
has two consequences. First, one has to deal with differential distributions which, in certain regions of phase 
space, display a bad convergence in perturbation theory. All-order resummations must then be performed in 
order to produce reliable predictions. Second, differential distributions of real hadrons depend unavoidably on 
some non-perturbative phenomenological inputs, fragmentation functions. Such inputs must be extracted from 
data and matched to the perturbative theory in a proper way, pretty much like parton distribution functions of 
light quarks and gluons are. 

In the following sections we review the state of the art of theoretical calculations of heavy quark produc- 
tion in a number of high energy processes, pointing out similarities and differences. In particolar, resummations 
aimed at improving the theoretical description of heavy quark production at large transverse momentum or large 
photon virtuality in DIS (Section 3), small centre-of-mass energy (Section 5) and large centre of-mass energy 
(Section 6) are described in some detail. 



3 Collinear resummations and heavy quark PDFs 

Perturbative calculations of heavy quark production contain badly converging logarithmic terms of quasi- 
coUinear origin in higher orders when a second energy scale is present and it is much larger than the heavy 
quark mass m. Examples are the (square root of the) photon virtuality in DIS and the transverse momen- 
tum pt in either hadroproduction or photoproduction. Naming generically E the large scale, we can write 
schematically the cross section for the production of the heavy quark Q as 



aQ{E,m) = (70 



1 + 

V n=l 



k=0 



(1) 



where uo stands for the Born cross section, and the coefficients Cnk can contain constants as well as functions 
of m and E, vanishing as powers of m/E when E m. 

Resummation approaches bear many different names, (ZM-VFNS, ACOT, FONLL, BSMN to name but 
a few) but they all share the goal of resumming leading (a" ln"(£^^/m^), LL) and sometimes also next-to- 
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leading (a" ln'^~^{E'^ /m?), NLL) logarithmic terms to all orders in the cross section above. This is achieved 
by discarding power suppressed m/E terms, and factoring all the logarithms into a resummation factor, to be 
obtained via Altarelli-Parisi evolution of an initial condition set at the heavy quark mass scale, 

(T^^^(E, m) = aoC{E, ^i)f{^i, m) = aoC{E, fi)E{n, mo)/(mo, m) , (2) 

where and no represent artificial factorization scales, to be taken of order E and m respectively. The 'prod- 
ucts' between the various functions actually hide convolution operations with respect to momentum fractions, 
not expUcitly shown as arguments. C{E, fi) is a perturbatively calculable coefficient function, which does not 
contain large logarithms thanks to the choice /i ~ E'. The function /(/i,m) can represent either a parton 
distribution or a fragmentation function for a heavy quark, and contains the resummation of the coUinear log- 
arithms. Due to the large heavy quark mass, its initial condition fino, m) can be calculated in perturbation 
theory [14, 15]: this is the distinctive feature that sets heavy quark parton and fragmentation functions apart 
from hght flavour ones, whose initial conditions are instead entirely non-perturbative and must be fitted to data. 

Once a massless but resummed result, valid in the E ^ m region, is obtained, one would like to 
interpolate it with a fixed order cross section, valid instead in the iiJ ~ m region, so as to retain predictivity 
over the whole E range. 

The differences between the various approaches are then to be found essentially in two points: 

• the perturbative order to which the initial condition /(/xq, m) is evaluated, and the perturbative accuracy 
of the evolution; 

• the way the matching with the fixed order calculation is performed. 

We summarize below the features of the most commonly used implementations. 

3.1 ACOT - Aivazis, CoDins, Olness, 'Ring 

This approach was the first to try to improve the prediction of the heavy quark structure functions (Q^, m^) 
and i^l (Q^' "^b) large » m^, m^, by moving potentially large logarithms ln((5^/m^) into heavy quark 
parton densities [16, 17]. A general all-order analysis of factorization for the total inclusive i<2(Q^) in this 
context was presented in [18]. 

3.2 SimpHfied ACOT and ACOT(x) 

author: S. Kretzer 

The original ACOT prescription [16, 17] has been simpUfied in [19] along lines suggested in [20, 18]. In a 
nutshell, diagrams with initial state heavy quark legs can be treated as if they represented massless quarks. More 
generally, the diagrams can be manipulated by power suppressed terms provided that higher order diagrams are 
regularized consistently. ACOT(x) [21, 22] explores this freedom to improve on the threshold behaviour of 
partonic heavy quark schemes by enforcing the physical pair-production threshold on a term-by-term basis. 
Heuristically, it comes down to a simple re-scaling of Bjorken-x, i.e. in LO 

Ff oc c(x)|^=^3.(i+4„2/Q2) . (3) 

Physical arguments -mostly kinematic- have been given in [21, 23, 22], here we will establish the connection 
with the FONLL terminology of Section 1.3.3 while focusing on the inclusive DIS process. Much of the 
following has appeared before, in one form or another, in the literature [16, 17, 18, 24, 25, 26, 27, 28, 19]. 

We formulate ACOT(x) as an explicit manipulation of resummed terms of the perturbative series. We 
follow [24] in notation and add an O (ct^) fixed order (FO) calculation to an all order collinearly resummed 
(RS) result. In RS heavy quark mass dependence other than logarithmic is neglected. When we remove double- 
counting terms from FO -i- RS the zero mass Umit (FOMO) of the FO calculation will be required as an auxiliary 
quantity. Just as in RS, only asymptotic mass logarithms are retained in FOMO. We write therefore, as usual, 

^ACOT^Q^ m) = FO + (RS - FOMO) x G (4) 
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where G is an arbitrary operation which behaves Uke G = l + O (^^^.In [24] G was chosen to be an overall 

multiplicative factor. More generally, it can be seen as an operation which only modifies, with 0{m^ /Q"^) 
power-suppressed terms, perturbative coefficients beyond those which have been explicitly calculated, and 
which are therefore unknown anyway. Any choice for G with this behaviour is therefore legitimate. 

To motivate the ACOT(x) choice for G we first re-write more exphcitly the three terms given above in 
the case of inclusive DIS: 

FO = as9^H{Q,m) (5) 
FOMO = as(^g^Pi^hn^+g^Cg^ (6) 

RS = c{x) + as{g0Cg + c(g)Cq) (7) 

where H{Q, m) is the massive coefficient function for the FO gluon fusion process, Cg and Cq are the gluon 
and light quark coefficient functions (the MS scheme is implied), and g and c are the gluon and charm (i.e. 
heavy quark) parton distribution functions (both the coefficient functions and the PDFs depend, of course, 
on the factorization scale ji ~ Q). is the leading order AltarelU-Parisi spUtting vertex. The symbol 
(§) = d^/^ ... denotes a threshold-respecting convolution integral. One can convince oneself that the standard 
convolution 0, with a; — > x the lower limit of integration, only differs by (8) by power-suppressed terms, 
(8) = + 0(mV<5^)- 

The combined result (4) reads now 



^ACOT^Q^ m) = FO + (RS - FOMO) x G 



= as g + 



c(x) -asg^ Pj°) \n^ + asC^Cq 



xG, (8) 



and we recognize the Kramer-Olness-Soper simplified ACOT framework of [19]^ if we set G = 1. Different 
choices for G can still be made, but natural demands are that: 

• In kinematic regions where FO represents the relevant physics (i.e. Q ~ m), G should efficiently sup- 
press uncontrolled spurious higher order terms in the square bracket of eq.(8). 

• For computational efficiency, the simple c{x) term alone should provide an optimized effective O (a^) 
approximation. 

The ACOT(x) scheme implements these requests by making an implicit choice for G which corresponds to 
writing 



^ACOT(x) (g^ ^) ^ po + (RS - FOMO) X G 



= as g + 



c(x) -asg In + a, c (8)C, 
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(9) 



Further details on ACOT(x) can be found in [21, 22, 23]. These articles also contain a more intuitive perspective 
of ACOT(x). Moreover, [22] describes a PDF set that is consistent with ACOT(x) applications. 



3.3 BSMN - Buza, Smith, Matiounine, van Neerven 

In Refs. [29, 30, 31, 32, 33] the treatment of heavy quarks as a parton was fully explored through next-to- 
next-leading order (NNLO), based on a precise two-loop analysis of the heavy quark structure functions from 
an operator point of view. This analysis yielded a number of results. One result is important beyond the 
observable at hand: the authors obtained the complete set of NNLO matching conditions for parton evolution 
across flavor thresholds. They found that, unhke at NLO, the matching conditions are discontinuous at the 
flavor thresholds. These conditions are necessary for any NNLO calculation at the LHC, and have already been 
implemented in a number of evolution packages [34, 35]. 

'See Eqs. (7), (8) there. General choices for G correspond to the discussion above these equations. 
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Furthermore, their two-loop calculations explicitly showed that the heavy quark structure functions in 

such a variable flavor approach are not infrared safe: one needs to either define a heavy quark-jet structure 
function, or introduce a fragmentation function to absorb the uncancelled divergence. In either case, a set of 
contributions to the inclusive light parton structure functions must be included at NNLO. 

A dedicated analysis [36] for charm electroproduction showed that even at very large one could not 
distinguish the fixed order NLO calculation of [10] and the NNLO VFNS calculations of [31], given the exper- 
imental data available in the year 2000. This demonstrates the possibihty that the large logarithms ln(Q^/m^) 
together with small coefficients can in the end have Uttle weight in the overall hadronic cross section. 



3.4 FONLL - Fixed Order plus Next-to-Leading Log resummation 

This approach was developed for improving the large-pT differential cross section for heavy quark production 
in hadron-hadron collisions [37]. It was successively extended to photoproduction [38], and in a second phase a 
matching to the fixed order NLO calculations was performed [24, 39]. The FONLL acronym refers specifically 

to the matched version. 

From the point of view of perturbative logarithms, it contains a NLO-accurate initial condition and full 
NLL evolution. It therefore reproduces the full NLL structure of the NLO calculation, and resums to all orders 
the large logarithms with NLL accuracy. 

The matching with the fixed order result is performed according to the following master formula (see 
eq.(16) of [24]): 

^FONLL^^^^ rn) = FO + (RS - FOM0)G{m,pT) , (10) 

where FO stands for the NLO fixed order massive calculation, FOMO for its m/px — * hmit (where however 
In pT/m terms and non- vanishing terms are kept), and RS for the massless, resummed calculation^. The 
RS— FOMO subtraction is meant to cancel the terms which are present in both RS and FO. This difference starts 
therefore at order with respect to the Born cross section: at large px it resums correctly the NLL terms, at 
small pt it only contains spurious terms, which are suppressed by the function G{m,pT) = Pt/{Pt + c^m^), 
with c = 5, in order to ensure a physically correct behaviour. The choice of the suppression factor was 
motivated in [24] by the observation that the massless limit starts to approach the massive hadroproduction 
calculation at 0{al) only for pT > 5m. Below this value the massless limit returns unreliable results, and its 
contribution must therefore be suppressed. It is important to reahze that G{m,pT) only affects terms which 
are beyond the control of perturbation theory, and therefore it does not spoil the NLO+NLL accuracy. The 
choice to control such terms by means of an ad-hoc function might seem a somewhat unpleasant characteristic 
of this approach. However, it simply portraits the freedom one has in performing the matching, and does not 
represent a shortcoming of the approach: different matching procedures will simply make other imphcit or 
explicit choices for G{m,pT)- 

For the sake of making comparisons with other approaches easier, the formula (10) can be rewritten with 
some more details as follows: 

+ j ^i^j ^W>kxiPT)Dk^Q - Y ^i^j CTij^QxiPT, m;m-^0)\ G{m,pT) • (H) 

\ijkec+n ijeC J 

A number of ingredients needing definition have been introduced. The kernel cross sections aij^Qx{pT,m) 
are the massive NLO calculations for heavy quark production of Refs. [2, 3, 4, 5, 6]. When convoluted with 
the PDFs for light flavours Fi (i G £) they yield the FO term in eq. (10). The aij^Qx{pT, in;m ^ 0) terms 
represent the m ^ limit of the massive NLO cross sections, performed by sending to zero m/pT terms 
while preserving ln(px/m) contributions and non- vanishing constants. When convoluted with light flavour 



^This term might also be referred to as a 'zero-mass variable flavour number scheme' (ZM-VFNS) contribution. However this name, 
while by itself completely general, has been used in the past for specific approaches with different overall perturbative accuracies. We 
shall therefore avoid its use. It will be understood that 'RS' in this approach has full NLL accuracy. 



6 



PDFs they give FOMO. Finally, (^fj^kxiPT) the massless MS'-subtracted NLO cross section kernels given 
in [40]. In addition to the light flavour PDFs, they are also convoluted with the perturbatively-calculated 
parton distribution functions for the heavy quarks (i G Ti) and with the fragmentation functions describing the 
transformation of a parton into a heavy quark, D^^q [15], to give the term RS. 

The formula given above returns the differential cross section for heavy quark production, evaluated with 
NLO + NLL accuracy. In order to obtain the corresponding cross section for an observable heavy meson it must 
still be convoluted with the proper scale-independent non-perturbative fragmentation function, extracted from 
experimental data, describing the heavy quark — > heavy hadron transition: 

<^H {PT,m)=aQ [pT,m)L>Q_^H . (12) 

Phenomenological analyses of charm- and bottom-flavoured hadrons production within the FONLL approach 
have been given in [41, 42, 43, 44, 45]. 

3.5 GM-VFNS - General mass variable flavour number scheme 

This approach also combines a massless resummed calculation with a massive fixed order one, for predicting 
Pt distributions in hadron-hadron collisions. One difference with respect to FONLL is that this approach does 
not include the perturbative NLO parton-to-heavy-quark fragmentation functions D^^q. Rather, it directly 
convolutes a properly MS subtracted cross section (with mass terms also included, hence the 'general mass' 
name) with non-perturbative fragmentation functions for heavy mesons Dq^h^, fitted at LEP in a pure MS 
scheme. The cross section can be schematically written as 

cth [PT,m) = ^ l<il<j aij^Qx(pT,m)DQ_^jj + 2^ ^i^j f^ij-.kxKPTjJ^k-^H ' (l^) 

ijec ijkec+n 

where the 'massive-but-subtracted' cross section kernels aij-^QxiPT, m) are defined by 

^ij^QxiPT, m) = aij^QipT, m) - (7ij^Qx{pT, m;m^O) + af-^QxiPr) ■ (14) 

The new kernels aij^QxiPTjm) defined by this operation (of the form FO-FOMO-i-RS) can be convoluted 
with an evolved MS-subtracted fragmentation function, but they also retain power suppressed m/pT terms. It 
should also be noted that the sum in the second term of (13) only runs over contributions not already included 
in the first. 

Recalling the way the perturbative parton-to-heavy-quark D^^q fragmentation functions are defined 
in [15], setting 

D'^Zr = Dk^QD^^H, keC + n, (15) 

and comparing eqs.(13) and (11), it can be seen that the GM-VFNS master formula is a reshuffling of the 
FONLL one, up to higher-orders terms. 

Two comments are worth making. The first is that due to the absence of the perturbative D^^q terms, 
eq. (13) cannot reproduce the NLO heavy quark production cross section: even the normalization must be 
extracted from the experimental data. Eq. (11), on the other hand, can reproduce the heavy quark spectrum, 
and only the heavy quark heavy meson transition is fitted to data. The second remark concerns the higher 
order power suppressed terms: since GM-VNFS implicitly makes a different choice for the G{m,pT) function, 
the results from the two approaches might differ considerably in the pr ^ m region since, while formally 
suppressed, such terms can be numerically important. 

An example of a phenomenological apphcation of the GM-VFNS scheme is given below. 

3.6 Hadroproduction of heavy mesons in a massive VFNS 

authors: B.A. Kniehl and I. Schienbein 

Various approaches for next-to-leading-order (NLO) calculations in perturbative QCD have been applied 
to one-particle-inclusive hadroproduction of heavy mesons. The general-mass variable-flavor-number scheme 
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(GM-VFNS) devised by us in Ref. [46, 47] is closely related to the conventional massless variable-flavor- 
number scheme (ZM-VFNS), but keeps all rn^ /p^ terms in the hard-scattering cross sections, where m is 
the mass of the heavy quark and px the transverse momentum of the observed meson, in order to achieve 
better accuracy in the intermediate region px > m. The massive hard-scattering cross sections have been 
constructed in such a way that the conventional hard-scattering cross sections in the MS scheme are recovered 
in the limit pr ^ oo (or m ^ 0). The requirement to adjust the massive theory to the ZM-VFNS with MS 
subtraction is necessary, since all commonly used PDFs and FFs for heavy flavors are defined in this particular 
scheme. In this sense, this subtraction scheme is a consistent extension of the conventional ZM-VFNS for 
including charm-quark mass effects. It should be noted that our implementation of a GM-VFNS is similar to 
the ACOT scheme [16, 17], which has been extended to one-particle-inclusive production of B mesons a few 
years ago [48]. As explained in the second paper of Ref. [46, 47], there are small differences concerning the 
coUinear subtraction terms. Furthermore, in Ref. [48], the resummation of the final-state coUinear logarithms 
has been performed only to leading logarithmic accuracy. The field-theoretical foundation of a GM-VFNS has 
been provided a few years ago by a factorization proof including heavy-quark masses [18]. Therefore, it is 
possible to extract improved universal parton distribution functions (PDFs) [22] and fragmentation functions 
(FFs) [49] from fits employing massive hard-scattering cross sections. From this perspective, it is important to 
compute massive hard-scattering cross sections in a given massive scheme for all relevant processes. Explicit 
calculations in the original ACOT scheme have been performed in Ref. [50, 5 1] for inclusive and semi-inclusive 
deep-inelastic scattering (DIS). Furthermore, our calculation in Ref. [46, 47] for hadronic collisions completes 
earUer work in the GM-VFNS on L'-meson production in 77 and jp collisions [52, 53, 54], and it is planned 
to extend our analysis to the case of heavy-meson production in DIS. 

Next, we show predictions for the cross section da/dpT of £)*-meson production obtained in the GM- 
VFNS and the ZM-VFNS. The cross section has been scaled with p^ in order to arrive at a flatpT distribution, 
which is useful for visualizing the heavy-quark mass effects. The hard-scattering cross sections are convoluted 
with the (anti-)proton PDFs and FFs for the transition of the final-state parton into the observed D* meson. We 
use the CTEQ6M PDFs [55] and the FFs for D* mesons from Ref. [56]. As in the experimental analysis, the 
theoretical results are presented for the average -|- D*~)/2. We consider da/dpx at \/5 = 1.96 TeV as a 
function of pT with y integrated over the range —1.0 < y < 1.0. We take the charm mass to be m = 1.5 GeV 
and evaluate Q;i"^^^(/Xi?) with n/ = 4 and scale parameter = 328 MeV, corresponding to af\'mz) = 
0.1181. The results are presented in Fig. 1 for the GM-VFNS (black lines) and the ZM-VFNS (red hues) in 
comparison with CDF data [57]. The solid lines have been obtained with = /up = /Xp = rriT- The upper 
and lower dashed curves represent the maximum and minimum cross sections found by varying fi^ , //t , and /i^ 
independently within a factor of 2 up and down relative to the central values requiring for their ratios to satisfy 
the inequalities 0.5 < /^f/mRj A*f/mR) Mf//^f — 2- seen, for large values of px, the predictions of the 

GM-VFNS nicely converge to the corresponding results in the ZM-VFNS. Both approaches lead to reasonable 
descriptions of the data, but the inclusion of the positive mass effects clearly improves the agreement with the 
data. It should be noted that the mass effects are largest for the upper curves of the uncertainty band, which have 
been obtained with the smaller value of the renormalization scale implying a larger as{iJ,R). At pT = 5 GeV, 
one observes an increase of the massless cross section by about 35%. A more detailed comparison of the 
GM-VFNS with CDF data [57] including D^, D+, and D+ mesons can be found in Refs. [58, 59]. 

Residual sources of theoretical uncertainty include the variations of the charm mass and the employed 
PDF and FF sets. A variation of the value of the charm mass does not contribute much to the theoretical un- 
certainty. Also, the use of other up-to-date NLO proton PDF sets produces only minor differences. Concerning 
the choice of the NLO FF sets, we obtain results reduced by a factor of 1 .2-1 .3 when we use the NLO sets 
from Ref. [60], which is mainly caused by a considerably different gluon FF. A more detailed discussion can 
be found in Ref. [56]. 

4 A case study in collinear resummation: &-quark PDF from Z + b production at LHC 

authors: S. Diglio, F. Maltoni, F. Petrucci, A. Tonazzo and M. Verducci 
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Fig. 1: QCD predictions for one-particle-inclusive production of D* mesons at the Tevatron Run II in comparison with CDF data [57]. 
The results are shown for the average D* — [D*^ + D*^)/2. The solid lines have been obtained with /ir = /if = /il? = niT- The 
upper and lower dashed curves represent the maximum and minimum cross sections found by varying /^r,, /if, and /.ip independently 
within a factor of 2 up and down relative to the central values while keeping their ratios in the range 0.5 < /.ip/ziR, /ip/z^iR, /if/mf ^ 2. 



4.1 Introduction 

The discovery of new physics at LHC will probably rely on the detailed understanding of standard-model 
background processes. Outstanding among these is the production of weak bosons {W, Z) in association with 
jets, one or more of which contains a heavy quark {Q = c, h). The prime example is the discovery of the top 
quark at the Fermilab Tevatron, which required a thorough understanding of the T^+jets background, with one 
or more heavy-quark jets. The discovery of single-top-quark production via the weak interaction will require 
an even more sophisticated understanding of this background [61, 62]. 

For many processes involving production of heavy quarks, there are two ways (schemes) to perform 
the calculation in QCD: the fixed-flavor-scheme (FFS) and variable-flavor-scheme (VFS). The main practical 
difference between the two approaches is simple: in the VFS the heavy-quark can also be in the initial state, 
and in that case is assumed to be massless, while in the FFS it appears only as a final state (massive) particle. 
QCD factorisation tells us that if calculations could be performed at arbitrary high order, the two schemes 
would be equivalent. At fixed order, on the other hand, differences arise and one should choose that describing 
more effectively the kinematics of the process of interest. This freedom has sometimes created intense and 
fruitful debates among the QCD practitioners (see, e.g., Ref. [63] for a detailed comparison of Higgs boson 
production in association with bottom quarks). Here we just recall the main two reasons for using a heavy- 
quark distribution function. First, it resums coUinear logarithms of the form In Q /niQ to all orders, where Q is 
the scale of the hard scattering and mg is the mass of the heavy quark. Second, it simplifies the leading-order 
process, which often makes a higher-order calculation feasible. There are many processes in the standard model 
and in models beyond it, such as SUSY, that are better described using a bottom in the initial state. In Table 1, 
we give a non-exhaustive list of processes that will be relevant for QCD, EW and SUSY studies at the LHC, 
and the QCD order at which they are known. 

At present the b distribution function is derived perturbatively from the gluon distribution function [17, 
18, 55, 34]. Recently, direct, albeit not very precise, measurements of FI have become available that are 
compatible with the perturbative determination [64, 65]. In the light of its phenomenological importance, a 
better direct determination of the h distribution function is certainly desirable. 

To this aim it has been proposed to use the associated production of a photon and a &-jet via gh —>■ at 
the LHC [66]. This measurement suffers from two main limitations. The first is the large contamination from 
charm which has a much larger cross section due to both the pdf and the electromagnetic coupling. The second 
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NLO 


Charged Higgs 


gb tH- 


SUSY Higgs discovery 


NLO 



Table 1: Process relevant for SM measurements and SUSY discoveries at the LHC which entail the use of bottom in the initial state. 
All of them are known at least at NLO accuracy. 



is that the theoretical prediction at NLO for an isolated photon is uncertain, due to necessity of introducing a 
photon fragmentation function, which is at present poorly known. 

In this note we follow the suggestion of Ref. [67] to use Z production in association with a 6-jet to extract 
information on the fe-pdf. At leading order, it proceeds via gb —>■ Zb, as shown in Fig. 2. This process is known 
at NLO, including 7/Z interference effects. The advantages of using a 7/Z decaying into leptons with respect 
to a real photon are noticeable. The NLO cross section is theoretically very well known and, apart from the 
PDF's, free of non-perturbative inputs. In addition, the competing process gc Zc is suppressed by the ratio 
of the couplings of the charm and the bottom to the Z, and makes the 6-pdf determination much cleaner. 




Fig. 2: Leading Order Feynman diagrams for associated production of a Z boson and a single high-pT heavy quark (Q = c, b). 

The DO Experiment at Tevatron has recently measured the cross-section ratio a{Z + b) /a{Z + jet) [68], 
and their result is consistent with the NLO calculation. 

As pointed out in [67], the measurement of this process at the LHC should be even more interesting 
because the contribution of the leading order process, sensitive to the b content of the proton, is more relevant 
than at the Tevatron. In addition, the total cross-section is larger by a factor 50, and the relative contribution of 
background processes, mainly Z + c, is smaller. These features are summarised in table 2, taken from Ref. [67]. 

Besides the possibility of extracting the 6-pdf, Z + b represents also a benchmark and in some cases 
a background to the search of the Higgs boson, when it is produced in association with a single high-pT- b 
quark [63] : the dominant leading-order subprocess for the production of a Higgs boson via its coupling to the 
bis bb ^ h; however, if the presence of a single b with high px is demanded, the dominant process becomes 
gb hb, with cross-sections of the order of tens of fb. The h can then decay to the same final states as the Z; 
in particular, the decay h —>■ is enhanced in some models [69, 70, 71]. 

A preliminary analysis on the potential of the ATLAS experiment to measure the Z+b-jet production at 
the LHC is presented in the following. 
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Cross sections (pb) 


Tevatron 


LHC 


Process 


ZQ inclusive 


gb^ Zb 
gb Zbb 


13.4 ±0.9 ±0.8 ±0.8 
6.83 


in4n+70+™ +30 
49.2 


gc Zc 
gc — ^ Zee 


20.3liJ ± O-l^il 
13.8 


1390 ± 100l?°l^° 
89.7 




Zj inclusive 


qq -> Zg, gq Zq 


-10-10+44+9+7 
iUiU_4o_2_i2 


1 c^87n+900+60 +300 
± dO / u_goo -300-500 



Table 2: Next-to-leading-order inclusive cross sections (pb) for Z-boson production in association with heavy-quark jets at the Tevatron 
= 1.96 TeV pp) and the LHC (^s = 14 TeV pp). A jet lies in the range pr > 15 GeV/c and |?7| < 2 (Tevatron) or |77| < 2.5 
(LHC). ZQ indicates events containing a heavy quark, Zj events which do not contain a heavy quark. 

4.2 A study of LHC measurement potential 

A sample of Z+jet events generated using the PYTHIA Monte Carlo [72] was processed with a fast simulation 
of the ATLAS detector, the ATLFAST package [73]. Only decays of the Z boson to ix^ ii~ were taken into 
account. The signal was defined as the sample events containing a b quark with pT > 15 GeV/c and |r?| < 2.5. 
The background samples containing respectively a c quark within the same cuts, or a jet originating from a 
hght quark or a gluon in the same range, were considered separately. The NLO cross-sections computed in [67] 
were used for the signal and for these two classes of background, while the cross-section given by PYTHIA 
was taken for the other types of events. 

The experimental selection of Z+jet events with Z IJ^'^ required the detection of two muons of 
opposite charge with pT > 20 GeV/c and \r]\ < 2.5 and one hadronic jet. The presence of two Mgh-pr 
muons ensures the possibility to have high trigger efficiency on this type of events. In addition, to reject the 
contribution from virtual photons, the invariant mass M^^ of the muon pair was required to be close to the Z 
mass (80 GeV/c^ < M^^ < 105 GeV/c^). About 50% of signal events are retained after applying these cuts, 
the loss being equally due to the r] acceptance and to the pr cut. 

The selection of events where the jet originates from a b quark was based on two different tagging 
methods, as described in the following. Their complementarity is still to be studied in detail, however the 
comparison of two independent selections will be important to control the systematic uncertainties. 

The first method to select Z + b events was based solely on the presence of a third muon. Hadrons 
containing a b quark give origin to prompt muon decays in about 12% of the cases. The efficiency of this 
method, therefore, cannot exceed this value, however the background is also expected to be small. The "third 
muon", considered to be the muon from the b hadron decay, will in general be softer and closer to a jet than the 
muons from the Z decay. The distribution of the transverse momentum of the third muon in Z + j events is 
shown in Fig. 3. Different thresholds on the third muon px were considered for the final selection. 

The second analysis used an inclusive method for 6-tagging, based on the presence of secondary vertices 
and of tracks with high impact parameter with respect to the primary vertex, originated from the decay of the 
long-lived b hadrons. The ATLFAST package reproduces the ATLAS 6-tagging capabilities by applying the 
tagging efficiency on b jets and a mis-tag rate on non-6 jets on a statistical basis, according to the values set 
by the user to reproduce the actual detector performance. The efficiency of the inclusive 6-tagging on signal 
events, after the selection described above, is about 30%. The mistagging probabihty is about 4% on c-quark 
jets and 0.5% on light jets. 

The overall efficiency on signal events, the expected number of signal and background events with an 
integrated luminosity of 10 fb~^ and the expected purity of the selected samples are reported in table 3. With 
the fast simulation, the soft muon tagging capabihties are optimistic, in that full efficiency and no mis-tag are 
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Fig. 3: Distribution of tlie transverse momentum of the tliird muon in a Z+]et sample, for signal events (left) and for events with no b 
quark (right). 

assumed for the lepton identification; more realistic assumptions will be made when the study is carried on 
with the full detector simulation. The efficiency on signal events achieved with the inclusive 6-tagging method, 
where the results of the fast simulation are more realistic, is higher than with the soft muon tagging, while the 
purity of the selected sample is still quite good. Consistent results were obtained with a full simulation of the 
ATLAS experiment, on a small statistics sample. 

A better determination of the signal component in the selected sample will eventually be achieved by 
exploiting the information on the transverse momentum of the 6-jet or of the third muon. 



Cut 


Efficiency 


jyPT>'i-5 GcV, |-»?|<2.5 


^Pt>15 GcV, |»?|<2.5 


other 


Purity 


pt^> i GeV/c 


4% 


13990 


6270 





69% 


pi^>5 GeV/c 


3% 


11090 


5210 





69% 


> 6 GeV/c 


2.5% 


8430 


4180 





67% 


incl. 6-tag 


14% 


49500 


17400 


49600 


43% 



Table 3: Expected efficiency, statistics and purity in a data sample corresponding to an integrated luminosity of 10 fb~ , using the 
soft muon tagging with different thresholds on the muon transverse momentum and the inclusive 6-tagging. Nb denotes the number of 
expected signal events as defined in the text, Nc the number of selected events with a c jet with px > 15 GeV and \r]\ < 2.5, Mother the 
selected events from other processes. The statistical error on efficiencies and purities, due to the limited size of the simulated sample, 
is at the level of 1-2%. 

Given the large statistics of the available data samples, the measurement will be limited by systematic 
effects. 

The possibility to control the systematic effects directly from data samples has been explored, in partic- 
ular the evaluation of ^-tagging performance and of the residual background. 

The 6-tagging efficiency can be checked using 6-enriched samples. Based on previous experience at 
Tevatron and LEP, we can expect a relative uncertainty of about 5%. 

The background in the selected sample is mainly due to mis-tagged jets from c and light quarks. This 
can be controlled by looking at the number of 6-tagged jets in data samples that in principle should contain no 
6-jets at first order. VF-f-jet events, for example, will be available with large statistics and with jets covering the 
full pt range of the signal. It can therefore be expected to estimate the background from mis-tagging with a 
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Fig. 4: Systematics due to mis-tagging of 6-jets as evaluated from T4^+jet events. Left: relative error on background level per jet px 
bin. Right: pr distribution of jets in event selected as Z + 6; the error band on the background contribution represents the systematic 
uncertainty, as derived from the previous plot. 

relative uncertainty at the level of few percent, as shown by the plots in figure 4. 
4.3 Conclusions and outlook 

Z boson production in association with a 5-jet can provide information on the 6-pdf. 

A preliminary study of the Z + b channel using a fast simulation of the ATLAS detector has shown 
that this type of event will be observed with very high statistics and good purity at the LHC. Given the large 
statistics of the samples, the precision of the Z + b cross-section measurement will be limited by systematic 
effects. Some possibilities to evaluate systematic uncertainties directly form the data have been considered. An 
overall accuracy on the measurement at the level of 5% can be expected. 

The availability of large samples opens interesting possibilities for the study of differential distributions: 
for instance, measuring the cross-section as a function of the r] and px of the Z boson would allow for the 
measurement of the b PDF as a function of the momentum fraction carried by the quark inside the proton. 
These items are an additional topic for further studies. 

5 Soft-gluon resummation 

authors: T. O. Eynck and E. Laenen 

QCD factorizes long- and short distance dynamics in inclusive cross sections with initial state hadrons 
into non-perturbative, but universal parton distribution functions, and perturbatively calculable hard scattering 
functions. Large remnants of the long-distance dynamics occur near the threshold edge of phase space in the 
form of logarithmic distributions that are singular at the edge. Resummation [74, 75] of these effects organizes 
them to all orders in perturbation theory, and thereby extends the predictive power of QCD. 

Threshold resummation is now a well-established calculational scheme with systematically improvable 
accuracy. It allows organization of all subleading powers of the logarithmic enhancements, and can be consis- 
tently matched to finite order perturbation theory. Resummed expressions, which take the form of exponentiated 
integrals over functions of the running coupling, require however a prescription for their numerical evaluation 
to handle a Landau pole singularity of the coupling. But for this intrinsic ambiguity (which must cancel against 
ambiguities in power corrections), threshold resummation is just as systematically improvable as the standard 
coupling constant expansion. 

As stated earlier, the more differential a cross section, the better suited it is for phenomenology, because 
one may incorporate detector-specific acceptance cuts and thereby reduce the need for extrapolation. Therefore 
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we should like to better understand the behavior of threshold-resummed expressions for double-differential 
cross sections. A study for the inclusive threshold-resummed heavy quark structure function can be found in 
Ref. [76]. Here we examine the differential structure function for the reaction 



which we write as 



We define the invariants 



S 



(PF^{S,TuUi) 
dTi dUi 

= {p + qf = S'-Q\ T^ = {p-p^f-m\ 
= {q-pif-rri^, 84 = S' + T1+U1. 



(16) 



(17) 



The invariant mass squared of the final state X' is given by 

Mx, =m'^ + S4 



(18) 



(19) 



so that the elastic (threshold) limit for the subprocess (16) is approached by S4 0. It may be converted to the 
double-differential structure function in terms of the heavy quark transverse momentum and rapidity, e.g. 
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where e.g. [11] 



Pt 



d{p^Y dyQ dTi dUi 



S'TjUj + Q^T^ + Q^S'Ti 
g,2 



m 



(1/2) 



(20) 



(21) 



At the parton level one may define invariants equivalent to those in (18), which we will denote by using lower 
case. The order-by-order perturbation theory expansion for the partonic version of this distribution oj{s4, ti,ui) 
and its all-order resummation have the following schematic forms 
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(We have also computed 52(A); by including ever more gi functions in the exponent in Eq. (22) we can increase 
the parametric accuracy of the resummation.) The symbol U represents, in this case, the logarithmically 
singular plus-distributions 

In-i(p)' 
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with p = Si/rn?, or, after a Laplace transform J dpex.p{—Np) by In* N. The conversion to momentum space 
then requires a numerical inverse Laplace transform. For the case at hand one needs to compute 
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Fig. 5: Expandability of the resummed expressions for d^F^/dTidSi witli NLL exponent (ratio to LO) 

with c the intercept of the contour with the real N axis, and 4>g{N) the gluon density in moment space. We 
chose a toy density for the gluon PDF, and the minimal prescription [77] to perform the N integral. In Fig. 5 
we evaluate this expansion as a function of the recoil mass 5*4, and compare it to its finite order expansions. We 
keep the variable Ti fixed at the average of its minimum and maximum allowed value. Clearly, for reasonable 
values of ^4 the resummed result is akeady well-approximated by its 2nd and 3rd order expansions. 

Another way to evaluate the resummed expression is in terms of towers of [78] L = \nN. 



LO = hoo{as) 



h 



fc=i 



(26) 



where the indicated coefficients Ckj can be determined exactly. More accuracy here means including more 
subleading towers. This method is equivalent, but not identical to the minimal prescription method. In practice, 
one need only include the first 4 terms in each tower, the higher terms are vanishingly small. The ambiguities 
mentioned earlier are shifted to far-subleading towers in this approach. To exhibit the convergence of terms in 
the towers, it will be useful and illustrative to exhibit contributions of a particular order in the strong coupling 
and the large logarithms. We will employ the notation 



n'^lo 



/L 



N 



(27) 



for finite order results, where k indicates the order in the strong coupling, the subscript N denotes moments, 
and / expresses if only the leading term {I = 1, L^^at), or also the next-to-leading term {I = 2, L^~^) is 
included, etc. In Fig. 6 we see also in this approach a rapid convergence toward the tower-resummed result. 
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Fig. 6: Tower resummation at N'*LO - A;Ljv, k € {1, 2, 3} (N*LO - 2Ljv and N^LO - 3Ljv almost coincide). 



A more complete study of the relevance of threshold resummation for electroproduction of heavy quarks 
at HERA still awaits. We note that even if the size of the corrections does not cause much concern for the 
perturbative analysis of an observable, threshold resummation or its finite order approximations, often lead to 
a reduction of scale dependence [79], indeed also seen in Ref. [76]. 
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6 kf factorization 

authors: S.P. Baranov, H. Jung, A.V. Lipatov andN.P. Zotov 
6.1 Introduction 

The transverse momenta of the partons initiating a hard scattering process, like heavy quark production via 
79 ~^ QQ or 99 ~^ QQ lepto- (hadro-) production, respectively, is mainly generated by the QCD evolution, 
which can reach large values, in DGLAP up to the factorization scale, in BFKL/CCFM/LDC even larger. 

The typical transverse momenta of the gluons in a process gg ^ X for different masses M of the system 
X are shown in Fig. 7 as a function of the momentum fraction x of one of the gluons for LHC energies. The 
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Fig. 7: Average tranverse gluon momenta kt in processes gg ^ X for different masses M of the system X as a function of the 
momentum fraction of one of the gluons x. The thin Unes indicate the RMS spread of the distributions. In (b) is shown the definition 
of x,kt and M for a gluon induced process. 




transverse momenta can become large, so that they cannot be neglected. A theoretical approach, formulated 
for small x, which takes into account the tranverse momenta is the fct-factorization [80, 81] or semi-hard [82] 
approach. 

In fcj -factorization the cross section for any process pp ^ X can be written as: 

<^ = J dxidx2 J dktidkt2A{xi,kti,q)A{x2,kt2,q)^{xi,X2,ktikt2,q) (28) 

with A{x, kt, q) being the un-integrated (/c^ -dependent) parton density function uPDF, q defines the factoriza- 
tion scale and a is the partonic cross section. The off-shell matrix-elements a are calculated in [80, 81]. 

The effects of finite transverse momenta are present independent of the evolution scheme: uPDFs can be 
defined also for the DGLAP evolution. A more detailed discussion on these effects can be found in [83, 84]. 

It is interesting to note, that the fct-factorization approach (in LO a^) agrees very well with calculations 
performed in the collinear approach in NLO ag, which is shown in [85]. The main effect comes from a 
more realistic treatment of the kinematics already in LO, which otherwise has to be covered in NLO. The kt 
factorization approach, however, is strictly vaUd only at small x, where the virtuality of the exchanged gluons 
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can be identified with its tranverse momentum k'^ ~ —k^. The full expression for the virtuaUty is [86]: 

= (29) 

I — X I — X 

with m being the recoiling mass of the hadronic system except the hard scattering process, taking into account 
the history of the evolution process. For finite x the mass effects can be substantial. 



6.2 Open bb production and correlations at the LHC 

Heavy quark production in the fc^ -factorization approach at HERA and the Tevatron was considered already 
in many papers (see, for example, [82, 87, 88, 89, 90]). In Ref. [91] the /c^ -factorization approach was used 
for a more detailed analysis of the DO and CDF experimental data. The effect of the initial gluon tranverse 
momenta on the kinematics of the bb production at the LHC were investigated [92]. The renormalization and 
factorization scales were set equal to either the initial gluon virtualities, = = q^^ or = m^y, as is 
in the standard collinear QCD, and the quark mass of m5 = 4.75 GeV was used. 

In Fig. 8a we show the transverse momentum distributions of B mesons at LHC energies. The calculation 
was performed in the range \'q^\ < 1 and the Peterson fragmentation with e = 0.006 using the KMS [93] 
parameterization for the un-integrated gluon density (see [83, 84]). The prediction for the azimuthal correlations 
between the muons coming from B meson decays are shown in Fig. 86 with > 6 GeV and \'q>^\ < 2.5. The 
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Fig. 8: Prediction for B-meson production at the LHC using the KMS un-integrated gluon density. In a thepr distribution of S-mesons 
is shown. In b the azimuthal fijj, correlation coming from the B decays is shown. 

azimuthal correlations indicate an important theoretical difference between the collinear and /ct -factorization 
approaches. In the collinear approximation at parton level and leading order, the b quarks are be produced 
exactly back-to-back, which is clearly unphysical when the gluon is evolved up to a large enough scale. Only 
starting with NLO a significant deviation from the back-to-back scenario is found. Thus the NLO calculation 
has to correct for the wrong kinematics in LO together with higher order corrections, leading to large K factors. 
In the fct-factorization, the transverse momenta of the gluons are correctly treated already in LO. In the kt 
- factorization approach the NLO corrections are therefore expected to be much smaller, since here only the 
purely dynamical corrections have to be applied, whereas the kinematics are already correctly treated in LO. 



6.3 Quarkonium production and polarization at the LHC 

Since the initial gluons have non-zero transverse momenta, they are off-shell, and they have a longitudinal 
component in their polarization vector. Typically, the kt values of the two colfiding gluons are much different, 
as the parton evolution is equivalent to the random walk in the In ] kf \ plane, not in kt plane. Roughly speaking, 
the kt of one of the gluons can be neglected in comparison with that of the other. So, in the initial state we 
have one nearly on-sheU (transversely polarized) gluon and one off-shell (longitudinally polarized) gluon. After 
the interaction, they convert into one on-shell gluon and a heavy vector meson. Simple heficity conservation 
arguments show that the polarization of vector meson must be longitudinal, in contrast with the ordinary parton 
model, where the initial gluons are both on-shell. This effect has been already studied for the HERA [94] 
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and Tevatron [95] conditions. Fig.9a shows the predictions for the LHC energy obtained with KMS [93] 
parameterization for un-integrated gluon densities. The calculations are restricted to the pseudorapidity interval 
—2.5 < rjr < 2.5 and assume ATLAS "fiGfiS" trigger cut, which means one muon with pt > 6 GeV and 
another muon with > 3 GeV. 
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Fig. 9: Predictions of different theoretical approaches for quarkonium production. In (a) the fraction of longitudinally polarized T 
mesons is shown: solid histogram - collinear parton model, singlet + octet; dashed - A;* -factorization with KMS u.g.d.. In (6) the ratio 
of the production rates Xbi Ixhi is shown: solid histogram - collinear parton model, singlet + octet; dashed - kt factorization with KMS 
u.g.d. 




Important effects are also seen in the production of P-wave bottomium states with different spins Xb\ 
and Xbi- At the Tevatron energies, this process has been considered in Ref. [96], and the predictions for the 
LHC are presented in Fig.9fe. The P-wave states are assumed to be detected via the decay Xfc ^ It" + 7> 
with an additional requirement that the energy of the decay photon be greater than 2 GeV. The ratio of the 
production rates cr(xi)/c(X2) is qualitatively different in the A;t-factorization and the colUnear parton model. 
The underlying physics is clearly connected with gluon off-shellness. In the colhnear parton model, the relative 
suppression of xi states becomes stronger with increasing because of the increasing role of the color-octet 
contribution: in this approach, the leading-order fragmentation of an on-shell transversely polarized gluon into 
a vector meson is forbidden. In contrast with that, in the /cffactorization approach, the increase in the final state 
■pT is only due to the increasing transverse momenta (and virtualities) of the initial gluons, and, consequently, 
the suppression motivated by the Landau- Yang theorem becomes weaker at large pt- 



6.4 Associated Higgs + jets production at the LHC 

The dominant mechanism for Higgs production at the LHC is gluon-gluon fusion, and the calculations can be 
significantly simplified in the large top mass limit (Mu < 2Mtop) [97]. The differential cross section of the 
inclusive Higgs production pp ^ H + X inihe fct-factorization approach has been calculated in [98, 99] and 
can be written as: 

da(pp^H + X) f a'^ilJ.'^) GpV^ r 2 2l2 2^a x 

dyn J 2887r x,X2mls ^ ^ ^ ^^^^ 

,2 .,2x Ar^ 1,2 ..2nj,,2 ji,2 d{Aip) 



)A{X2, 



27r ' 

where Gp is the Fermi coupling constant, A{x,'k'^, /j?) is the un-integrated gluon distribution, Aip the az- 
imuthal angle between the momenta kiT and k2T, and the transverse momentum of the produced Higgs boson 
is Pt = kiT + k2T- It should be noted, that this process is particularly interesting in /c^-factorization, as the 
transverse momenta of the gluons are in the same order as their longitudinal momenta (~ 0(10 GeV)) [100]. 

The total inclusive Higgs production cross section at the LHC energies (^/s = 14 TeV) is plotted in 
Fig. 10(a) as a function of the Higgs mass in the mass range itlh = 100 — 200 GeV. The solid line is obtained 
by fixing both the factorization and renormahzation scales at the default value /j, = mn with the J2003 (set 
1) CCFM un-integrated gluon distribution [101]. In order to estimate the theoretical uncertainties, we take 
H = ^ruH and vary the scale parameter ^ between 1/2 and 2 about the default value ^ = 1. The uncertainty 
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band is presented by the upper and lower dashed Unes. We find that our central values agree very well with 
recent NNLO results [102]. 




100 150 200 0.0 1.0 2.0 3.0 
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Fig. 10: Higg production at the LHC. In (a) the total cross section for Higgs boson production as a function of Higgs mass is shown: 
the solid curve corresponds to the default scale ji = tuh, upper and lower dashed curves - /tt = ruH/'i and n = 2mif/2 respectively. 
In (6) the jet-jet azimuthal angle distribution in the Higgs+jet+jet production at y/s = 14 TeV. The kinematical cut |pjet t| > 20 GeV 
was applied for both jets. Solid and dashed lines correspond to the J2003 (set 1) and J2003 (set 2) [101] u.g.d. respectively. 

To demonstrate the capabilities of the /cj -factorization approach, we calculate the azimuthal angle A0 
distribution between the two final jets transverse momenta in the Higgs-f-jet-t-jet production process. Our results 
are shown in Fig. 10(6). The dip at = it/2 comes from the cos(A(/9) in eq.(30). In the approach presented 
here, the kt of the initial gluons is approximately compensated by the transverse momenta of the jets [103]: 
k.T ^ —PTjeU and, consequently, Acj) ^ Aip applying a cut-off |pjetT| > 20 GeV. This dip is characteristic 
for the loop-induced Higgs coupling to gluons in the framework of fixed-order perturbative QCD calculations 
[102]. Thus, we illustrate that the features usually interpreted as NNLO effects are reproduced in the kf 
factorization with LO matrix elements. 

However, we see a very large difference (about one order of magnitude) between the predictions based 
on the J2003 gluon densities set 1 and set 2 [101], showing the sensitivity to the shape of the un-integrated 
gluon density. 

6.5 Conclusions 

The finite kt of the initial state gluons significantly modifies the kinematics of the gluon-gluon fusion process 
and leads to nontrivial angular correlations between the final state heavy quarks. The longitudinal polarization 
of the initial off-shell gluons manifests in the longitudinal polarization of J/tp and T mesons at moderate pr 
and, also, affects the production rates of P-wave quarkonia. 

The predictions in fcrfactorization are very close to NNLO pQCD results for the inclusive Higgs pro- 
duction at the LHC, since the main part of high-order collinear pQCD corrections is already included in the 
fcj -factorization. In the A:^ -factorization approach the calculation of associated Higgs-njets production is much 
simpler than in the collinear factorization approach. However, the large scale dependence of our calculations 
(of the order of 20 — 50%) probably indicates the sensitivity to the unintegrated gluon distributions. 

7 Baryon charge transfer and production asymmetry of A°/ A° in Iiadron interactions 
author: O.I. Piskounova. 

7.1 Introduction to the QGSM 

The phenomenon of nonzero asymmetry of baryon production with nonbaryonic beams (7r,/x,e) was mentioned 
and explained in a few theoretical papers. Baryon charge can be transferred from proton or nucleus targets 
through the large rapidity gap with the string junction. In baryonic beam reactions {p,A, etc.) this effect is dis- 
played in the valuable baryon/antibaryon spectrum asymmetry at y = 0. Every theoretical discourse on baryon 
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charge transfer appeals to the value of the intercept, asj{0), that is an intercept of the Regge-trajectory of imag- 
inary particles which consists only of string junctions from baryon and antibaryon. Practically, the models that 
can account for this effect are only non perturbative QCD phenomenological models: the Dual Parton Model 
(DPM) [104] and the Quark Gluon String Model (QGSM) [105] as well as the DPMJET Monte Carlo expansion 
of these two models. Both analytical models are similar and they were based on the common Regge asymptotic 
presentation of constituent quark structure functions and string (quark) fragmentation functions. Here we are 
considering QGSM. In the comparison to the other models, QGSM accounts for many Pomeron exchanges. 
This approach works very well to give us the correct description of particle production cross sections at very 
high energies. The QGSM procedures of constructing of quark/diquark structure functions and fragmentation 
functions were presented in many previous publications. We take into consideration the ir-p reaction that gives 
similar asymmetries as the 7-p reaction. The spectra in this reaction are more sensitive to the baryon excess in 
the region of positive xp than the spectra of baryons in p-p collisions. 



7.2 Comparison with Experimental Data 

The asymmetry A{y) between the spectra of baryons and antibaryons is defined as: 

dN^"/dx-dN^''/dx 
^""^ " dN^ydx + dN^"/dx' ^ ^ 

The EHS and the NA49 experiments have presented evidence for a nonzero baryon production asymmetry in 
proton-proton fixed target interactions, measuring at y = values of the order of 0.5 - 0.3. In pion-proton 
interactions (E769) we can see the y dependence of the asymmetry and the measured asymmetry, which was 
multiplied by a factor of 2 in order to be compared with the pp asymmetry. 




Fig. 11: Asymmetry in A° and A" production and QGSM curves: asj(0)=0.5(dashed line) and asj(0)=0.9 (solid line). 

The data from these experiments can be presented in one plot for A{Ay), where Ay is the rapidity 
distance from interacting target-proton (see Fig. 11). It is seen that the points are situated on the same hne. 
If we add the data from proton-nucleus experiments (HERA-B and RHIC) they are still approximately on this 
hne. Only the STAR asymmetry point at ^/s = 130 GeV can be interpreted as a sign that the curve is bent. And 
the result of the HI experiment at HERA [106] calls certainly for a steeper curve. What means do we have in 
QGSM to describe this dependence? In Ref. [107] it was shown that the data of the E769 and HI experiments 
can not be described with the same value of asj{0): the points at lower energy correspond to q;5j(0)=0.5, 
while the HI point requires a5j(0)=0.9. 
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7.3 Summary 

The purpose of this contribution is to show the band of asymmetries that can be predicted for the LHC exper- 
iments between the two possibilities given above for asj(0). The results are shown in Fig. 11. The soUd Une 
represents the case of asj{0)=0.9. This curve fits the data at low energies (small Ay) due to varying the energy 
splitting between string junction and diquark configuration: 0.1*SJ+0.9*DQ. What we had to tune also was the 
fragmentation parameter af=0.15 instead of 0.2 accepted in previous papers. Also the curve for q;5j(0)=0.5 is 
shown in Fig. 11 with a dashed line. This line certainly doesn't fit the HI point and gives a negligible asym- 
metry at the energy of the LHC experiments. Finally, we have the prediction for strange baryon asymmetries 
at the LHC within the range: 0.003 < yl < 0.04. The same procedure has to be applied to charmed baryon 
asymmetry to get the predictions at LHC energy. 
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Experimental overview of heavy quark measurements at HERA 
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Abstract 

Experimental results on heavy flavour production at HERA are reviewed in the con- 
text of their relevance for future measurements at the LHC. 

1 Introduction 

Measurements of heavy flavour production at HERA can have significant impacts on the preparation and under- 
standing of heavy flavour measurements at the LHC, and on the understanding of background processes to new 
physics discoveries [1]. The purpose of this contribution is to summarize the current status of heavy flavour 
measurements at HERA, and provide an outlook on how they might improve in the near future. The relation 
of these measurements to measurements at the LHC will be covered in more detail in subsequent contribu- 
tions [2, 3, 4]. Since the top quark is too heavy to be produced at HERA with a significant rate, the term "heavy 
flavour" refers to b and c quarks only. The dominant diagram for heavy flavour production at HERA is shown 
in Fig. 1. The theory of heavy quark production at HERA is covered in the theoretical review section [5]. 



e± 




Fig. 1: Feynman graph for the production of a heavy quark pair via the leading order boson-gluon-fusion (BGF) process. 

The interest in heavy flavour production arises from several aspects. 

• Tagging a heavy flavour particle, e.g. inside a jet, establishes that this jet arises from a quark rather than 
a gluon. The number of possible QCD diagrams is thus reduced, and specific QCD final states can be 
studied more precisely than in inclusive measurements. This is even more true when both quarks of a 
heavy flavour quark pair can be tagged. 

• The charm and beauty masses {rrih, rric ^ ^qcd) provide energy scales which are large enough to allow 
perturbative calculations using a "massive" scheme [6, 7]. All QCD-processes involving heavy quarks 
should thus be reliably calculable. However, these mass scales often compete with other scales occurring 
in the same process, such as the transverse momentum (px) of the heavy quarks, or the virtuality of the 
exchanged photon, Q^. Since the perturbative expansion can not be optimized for all scales at once, 
additional theoretical uncertainties enter which reduce the reliability of the predictions. If one of the 
competing scales (px, Q^) is much larger than the quark mass, approximations in which the heavy quarks 
are treated as massless [8, 9, 10, 11, 12, 13, 14] can improve the reliability. Mixed schemes [15, 16, 17] 
are also possible. Understanding and resolving these difficulties should contribute to the understanding 
of multi-scale problems in general. 

• Tagging the final state also constrains the initial state. Therefore, heavy flavour measurements can be used 
to measure or constrain parton density functions. In particular. Fig. 1 illustrates the direct sensitivity to 
the gluon density in the proton. Alternatively, in appropriate kinematic ranges, the initial state splitting 
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of the gluon or photon into a heavy quark pair can be absorbed into the parton density definition. If the 
mass can be neglected, the same diagram (or higher order variants of it) can be reinterpreted as a way to 
measure the heavy flavour content of the proton or of the photon. These can in turn be used to calculate 
cross sections for other processes, such as Higgs production at the LHC. 
• The production of "hidden" heavy flavour states (quarkonia) yields further insights into the interplay of 
(perturbative) heavy quark production and (non-perturbative) binding effects. 

At HERA, the fraction of charm production vs. inclusive QCD processes is of order 10% in the per- 
turbative QCD regime. Reasonably large samples can therefore be obtained despite the partially rather low 
tagging efficiency. Beauty production is suppressed with respect to charm production by the larger b mass, 
and by the smaller coupling to the photon. The resulting total cross section is two orders of magnitude smaller 
than the one for charm. Beauty studies at HERA are thus often limited by statistics. Kinematically, beauty 
production at HERA is similar to top production at LHC {nii,/ ^shera ~ rtit/ ^/slhc)- On the other hand, in 
the "interesting" physics region beauty is produced as copiously at the LHC as charm is at HERA. 

2 Open charm production 

Charmed mesons are tagged at HERA in different ways. A typical mass distribution for the "golden" channel 
Z)'^7r+, —>■ K^n^ (+ c.c.) is shown in Fig. 2 [20]. Despite the low branching ratio, this channel 
yields large statistics charm samples of high purity. Fig. 3 [21] shows a corresponding D* production cross 
section in photoproduction for different kinematic variables. In general, D* production is well described by 
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Fig. 2: Total inclusive D*^ sample obtained by ZEUS for the HERA I data period in the golden decay channel D* 

{K-t:+)7v+. 



next-to-leading order QCD predictions, although the data often prefer the upper edge of the theoretical error 
band. Some deviations are observed in particular regions of phase space. For instance, there are indications 
that forward (i.e. in the direction of the proton) charm production might be slightly larger than theoretical 
expectations (Fig. 3b). Also, there are regions of phase space which effectively require four-body final states 
which are not covered by NLO calculations (see Fig. 5 in [1]). In order to describe such phase space regions, 
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Fig. 3: D*^ single differential cross sections in photoproduction as function of the D*^ transverse momentum and pseudorapidity. The 
measurements are compared to NLO calculations in the massive (NLO), massless (NLL) and mixed (FONLL) scheme. 



either NNLO calculations, or parton shower extensions to NLO calculations such as MC@NLO [18, 19] will 
be needed. 

Other ways to tag charm include the reconstruction of a secondary vertex in a high resolution Micro- 
Vertex-Detector (MVD) in addition to the reconstruction of a charmed meson mass (Fig. 4) [22], or the recon- 
struction of inclusive multiple impact parameters resulting from the finite charm meson lifetime. A resulting 
cross section for Z)+ production is shown in Fig. 5. 

Since the charm mass of approximately 1.5 GeV is not very much above the threshold at which per- 
turbative calculations are believed to produce reliable results, the generally good agreement of perturbative 
QCD predictions with the data is highly nontrivial, and encouraging concerning the validity of corresponding 
predictions for the LHC. 
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Fig. 4: D+ mass peak in HI before (left) and after (right) a cut on the decay length significance. 
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. 5: Cross sections for D production in DIS. A leading order -I- parton shower QCD prediction (AROMA) is shown for comparison. 
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3 Open beauty production 

Open beauty production is detected at HERA using essentially three different methods, related to the large b 
mass or long b lifetime. 

• If a jet is built out of the fragmentation and decay products of a 6 meson/quark, the transverse momentum 
of the decay products with respect to the jet axis will be of order half the b mass. This is significantly 
larger than for decay products of charm particles, or the transverse momenta induced by non-perturbative 
fragmentation effects, which are both of order 1 GeV or less. This distribution of this transverse momen- 
tum, called pl^^ can thus be used to measure the beauty contribution to a given jet sample. 

• Due to the CKM-suppressed weak deacy of the b quark, the lifetime of b hadrons is longer than that 
of charmed particles. Furthermore, the larger decay angle due to the larger mass results in a higher 
significance of the decay signature. Detectors with a resolution in the 100 fi region or better can thus 
separate the beauty contribution from charm and light flavour contributions. 

• Again due to their mass, b hadrons take a larger fraction of the available energy in the fragmentation pro- 
cess. Furthermore, they produce decay products with sizeable transverse momentum even when produced 
close to the kinematic threshold. Simple lower cuts on the transverse momentum of such decay products 
therefore enrich the beauty content of a sample. Applying such cuts on two different decay products 
(double tagging) often sufficiently enriches the beauty content such that the remaining background can 
be eliminated or measured by studying the correlation between these decay products. 




Fig. 6: Distributions of the impact parameter S of the muon track (left) and the transverse muon momentum relative to the axis 
of the associated jet (right) in HI. Also shown are the estimated contributions of events arising from b quarks (dashed line), c quarks 
(dotted line) and the light quarks (dash-dotted line). 

An example for an analysis using the first two methods with muons from semileptonic b decays is shown 
in Fig. 6 [23]. Some cross sections resulting from this type of analysis are shown in Fig. 3 of [1]. In general, 
reasonable agreement is observed between the data and corresponding NLO QCD predictions, although, as in 
the charm case, the data tend to prefer the upper edge of the theoretical error band. In some regions of phase 
space, e.g. at low pi^ or high ry^ differences of up to two standard deviations are observed. More precise 
measurements (section 8) will be needed to decide whether these deviations are really significant. 

An example for an analysis using the 2nd method only is shown in Fig. 7 [24], while an example for an 
analysis using the third method is shown in Fig. 8 [25]. 

Fig. 9 shows a summary of the data/theory comparison for all HERA beauty results as a function of 
Q^. For the measurements sensitive to b quarks with p^ ~ mi, or lower (black points) there is a trend that the 
"massive" NLO QCD predictions [7] tend to underestimate the b production rate at very low Q^. Depending on 
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Fig. 7: Upper left: significance Si = S/a{5) distribution per event for events that contain one selected track associated to the jet axis. 
Lower left: significance S2 ~ S/a{5) distribution per event of the track with the second highest absolute significance for events with 
> 2 selected tracks associated to the jet. Right: 5*1 and 5*2 distributions after subtracting the negatvie bins in the Si and S2 distributions 
from the positive. 



the chosen set of structure functions and parameters, a "mixed" prediction (VFNS) [16, 17] might describe the 
data better. For the higher px measurements (red/grey points), no clear trend is observed. Note that theoretical 
errors, which are typically of order 30%, are not shown. Fig. 10 shows a similar compilation for all HERA 
measurements in photoproduction (Q^ < 1 GeV), now as a function of the b quark px- A similar trend is 
observed towards low pT (but note that several measurements appear in both figures). Again, more precise 
measurements are needed to determine whether these trends can be confirmed. 
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4 Quarkonium production 

Inelastic heavy quarkonia, like open charm and beauty production, are produced at HERA via the process of 
photon-gluon fusion. The two charm or beauty quarks hadronize to form a charmonium or bottomonium state. 

A number of models have been suggested to describe inelastic quarkonium production in the framework 
of perturbative QCD, such as the color-singlet model (CSM) [33, 34], the color-evaporation model [35, 36] 
and soft color interactions [37]. Comprehensive reports on the physics of charmonium production are available 
[38, 39]. 

Recently the ansatz of non-relativistic quantum chromodynamics (NRQCD) factorization was intro- 
duced. In the NRQCD approach non-perturbative effects associated with the binding of a qq pair into a quarko- 
nium are factored into NRQCD matrix elements that scale in a definite manner with the typical relative velocity 
V of the heavy quark in the quarkonium. This way, colour octet quark anti-quark states, carrying different 
angular momenta and color charges than the quarkonium, can contribute to the charmonium production cross 
section. Theoretical calculations based on the NRQCD factorization approach [40, 41, 42] are available in 
leading order [43, 44, 45, 46, 47, 48]. In the NRQCD factorization approach the size of the color octet con- 
tributions, which are described by long distance matrix elements (LDME), are additional free parameters and 
have been determined in fits to the Tevatron data [49]. The NRQCD factorization approach incorporates the 
color singlet model i.e. the state qq[\,^ Si] which is recovered in the limit in which the long distance matrix 
elements for other qq states tend to zero. 

At HERA, cross sections measurements for photoproduction of J/ijj and ipi'^S) and for electroproduc- 
tion of J/ijj mesons have been performed [52, 53, 54, 55]. Bottomonium data are not available due to statistical 
limitations of the data. 

For J and 'ip{2S) photoproduction, calculations of the color-singlet contribution are available to next- 
to-leading order perturbation theory [50, 51]. Calculations which include the color octet contributions as pre- 
dicted by NRQCD are available in leading order. 

Figure 11 shows the measurements of the J/ip photoproduction cross section by the HI collabora- 
tion [52] and the ZEUS collaboration [53] which are in good agreement with each other. The variable z 
(left figure) denotes the fraction of the photon energy in the proton rest frame that is transferred to the J/ip. 
Reasonable agreement is found between the HERA data and the NRQCD factorization ansatz in leading order 
(LO, CS-i-CO). The uncertainty indicated by the open band is due to the uncertainty in the color-octet NRQCD 




Fig. 1 1 : Differential chiarmonium photoproduction cross sections as measured by HI and ZEUS in comparison to calculations from LO 
NRQCD factorization (open band), NLO Color singlet contribution (shaded band) and LO color-singlet contribution (dotted line). 

matrix elements. In contrast, the shaded band shows the calculation of the color-singlet contribution (NLO, CS) 
which is performed to next-to-leading order in ag [50, 51]. This NLO, CS contribution alone describes the data 
quite well without inclusion of color-octet contributions. Comparison between the NLO,CS prediction (shaded 
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Fig. 12: Differential Charmonium photoproduction in comparison with a prediction using the color singlet-model and kt factorization 
as implemented in the Monte Carlo generator CASCADE. 



band) and the LO,CS prediction (dotted line) shows that the NLO corrections are crucial for the description of 
the HERA photoproduction data. 

Charmonium production cross sections have also been calculated in the kt factorization approach [56, 
57, 58]. In these calculations the color-singlet model is used to describe the formation of the charmonium 
state. Figure 12 shows a comparison of the HI data with the predictions from the kt factorization approach as 
implemented in the Monte Carlo generator CASCADE [59]. Good agreement is observed between data and 
predictions for z < 0.8. At high z values, the CASCADE calculation underestimates the cross section. The 
CASCADE predictions for the the pf,^^ dependence of the cross section fit the data considerably better than the 
LO,CS calculation in the collinear factorization approach (dotted curve in fig. 11). 

In fig 13 the differential cross sections for electroproduction of J /ip mesons as measured by HI [54] and 
ZEUS [55] are shown as a function of z and compared with predictions from the color singlet model (shaded 
band), with the NRQCD calculation [60] (CS-i-CO, open band), and also with calculations in the kt factorization 
approach (dotted Une) as provided by [58] and implemented in the Monte Carlo generator CASCADE (dash- 
dotted line). 

In the left figure the data are seen to agree well with the predictions using the color singlet model (shaded 
band and lines) while the full NRQCD calculation (open band), including color-octet contributions is wrong 
in shape and normaUzation. The agreement deteriorates when the cut > 1 GeV is applied (right fig. 13). 
This cut is justified, however, as towards small perturbation theory becomes increasingly unreliable due to 

collinear singularities for the contributions e + g ^ e + cc[n] + g with n=^S^^ and '^Pj^^ [60]. 

In conclusion, NRQCD, as presently available in leading order, does not give a satisfactory description 
of the HERA data. In contrast, the color singlet model shows a reasonable description of the HERA data, 
when implemented in calculations to next-to-leading order perturbation theory or in calculations in which the 
fct-factorization approach is used. 

5 Charm and Beauty contributions to structure functions 

To a good approximation, except at very high Q^, the cross section for inclusive deep inelastic electron scat- 
tering off the proton at HERA can be described in terms of a single proton structure function F2 (for formula 
see [1]). This structure function only depends on the photon virtuality, Q^, and on the Bjorken scaUng variable 
X. Assuming that the electron scatters off a single quark in the proton (0th order QCD, quark-parton model) x 
can be reinterpreted as the fraction of the proton momentum carried by the struck quark. This is a reasonable 
approximation for the light quark content of the proton. 

For heavy quarks, the situation is a bit more complicated. Due to the heavy quark mass, on-shell heavy 
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Fig. 13: Differential cross sections da/dz without (left) and with (right) a cut on p*^ > 1 GeV. The data are compared to the NRQCD 
calculation (CS+CO, open band), the color-singlet contribution (CS, shaded band), with a prediction in the kt factorization approach 
assuming the CSM (dotted line) and with the Monte Carlo generator CASCADE (dash-dotted line). 



quarks can not exist within the proton. Rather, the dominant process for heavy quark production is the 1st order 
QCD BGF process depicted in Fig. 1. However, this process (and other higher order processes) still contributes 
to electron scattering, and hence to F2. This can be interpreted in two ways. 

In the massive approximation, heavy quarks are treated as being produced dynamically in the scattering 
process. The heavy quark contribution to F2, frequently denoted as Fg'^ and therefore indirectly measures 
the gluon content of the proton. If is large enough such that the quark mass can be neglected {Q^ 3> 
(2m(3)^), the splitting of the gluon into a heavy quark pair can be reinterpreted to occur within the proton. Fg'^ 
and then measure the occurrence of virtual heavy quark pairs in the proton, or the "heavy quark structure" 
of the proton. 

For charm production, the condition » (2mQ)^ is valid for a large part of the HERA phase space. 
For beauty, it is only satisfied at very large Q^. This is also the region most interesting for physics at the LHC. 

Similar arguments hold for the heavy quark structure of the photon. 

As an example. Fig. 14 [61, 62, 63, 24, 31] shows ^2*^ as measured by the ZEUS and HI collaborations. 
A different representation of these results is shown in Fig. 6 of [1]. There, also i^f^ is shown. Good agreement 
is observed with QCD predictions. Parametrizations of heavy quark densities of the proton at LHC energies 
should therefore be valid within their respective errors. 
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Fig. 14: F2 results as a function of x in bins of Q^, from the HI and ZEUS D*^ analyses and from the Hi inclusive lifetime tagging 
measurements The data are compared to a NLO prediction using the ZEUS NLO fit results for the proton parton densities. 

6 Charm fragmentation 

The large cross section for charm production at HERA allows measurements of charm fragmentation which are 
very competitive with e+e" measurements. As this topic is covered very nicely in [1] and [2] it is not treated 
further here. 

7 Quark-antiquark correlations 

Heavy quarks are always produced in pairs. An interesting way to check QCD is thus to verify whether the 
kinematic correlations between the quark pair are correctly described by QCD. 

Fig. 15 shows different interpretations of the same higher order beauty production process. These 
different interpretations partially manifest themselves in different kinematic regions of beauty production phase 
space. If the highest virtuality part of the process occurs in the leading order BGF-like subprocess (left), the 
two b quarks will be almost back-to-back in the detector transverse plane. The two extra gluons can either 
be reabsorbed into the proton structure, recovering the original BGF graph, or manifest themselves as visible 
"parton shower" activity in the direction of the proton. Alternatively, if the dominant leading order subprocess 
is gluon exchange with one of the b quarks (right, "flavour excitation in the photon"), this b quark will recoil 
against a gluon jet. At sufficiantly large momentum transfer (rare at HERA), the second b quark can be treated as 
a "spectator", and will approximately follow the initial photon direction. At next-to-leading order, contributions 
to both processes are described by the same Feynman graph, but the two extreme kinematic cases (and all 
variants in between) are still included. If both heavy quarks are tagged, these different kinematic regions can 
be distinguished by measuring the momentum and angular correlations between the two quarks. 

Fig. 16 [25] shows the angular correlations between the two muons originating from different b quarks of 
a bb pair. Reasonable agreement is observed with QCD predictions. The predominantly back-to-back topology 
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Fig. 15: Example for higher order Feynman graph for beauty production. Different interpretations of the graph in terms of NLO or LO 
matrix elements plus parton showers are highlighted. Depending on the kinematics and the scheme chosen, part of the gluons could also 
be reabsorbed into the proton structure function definition, and/or the ybb vertex could be interpreted as part of the photon structure. 



comfirms the dominance of the BGF-Uke contribution. 
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Fig. 16: Differential cross section da / dA(j>'^'^ for dimuon events from bb decays in which each muon originates from a different 6(6) 
quark. The data (solid dots) are compared to the leading order + parton shower generators PYTHIA and RAPGAP (histogram) and to 
massive NLO QCD predictions (shaded band). 
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8 HERA II prospects 

Both the HERA collider and its detectors have been upgraded in 2001/2 to provide more luminosity with po- 
larized electron beams, and improve heavy flavour detection. This program is called HERA II. The luminosity 
accumulated so far already exceeds the HERA I luminosity. An integrated luminosity up to 700 pb^^ is ex- 
pected at the end of the HERA program in 2007. This enhances the statistics for many studies by almost an 
order of magnitude with respect to HERA I. The improved detectors offer further handles for improved heavy 
flavour measurements. HI has improved the forward coverage of its Micro- Vertex-Detector [68], and added 
a Fast Track Trigger [67]. ZEUS has implemented a Micro- Vertex-Detector (MVD) [65] for the first time for 
HERA II, and has added an upgraded forward tracking detector [66]. These improvements allow the application 
of measurement techniques which could not be used at HERA I, and can be used to improve the data quality, 
add additional statistics, and/or cover new phase space regions. 

New detectors require time to fully understand their systematics, but first preliminary results have already 
been obtained. Fig. 17 [64] shows the cross section for beauty production obtained using the new ZEUS MVD 
with the first 33 pb~^ of HERA II data, compared to the HERA I result. Good agreement is observed. 

The measurements which will profit most from the improved HERA II data sets include double dif- 
ferential measurements such as the beauty and charm contributions to the proton structure function F2, and 
multi-tag measurements to explicitly study quark-quark correlations. Statistical improvements of at least one 
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Fig. 17: Differential cross section as function of muon pr for dimuon -I- jet events in photoproduction. Preliminary results from the first 
33 pb^^ of HERA II data are compared to HERA I results and QCD predictions. 



order of magnitude can be expected when the increased luminosity and improved measurement techniques are 
combined. 

9 Conclusions 

Heavy flavour production at HERA is a very active field of research yielding multiple insights into the applica- 
bility of perturbative QCD. The problem of multiple scales complicates the perturbative expansions and limits 
the achievable theoretical precision. In general, QCD predictions agree well with the data, although indications 
for deviations persist in specific regions of phase space. Some of these might be attributable to missing NNLO 
or even higher order contributions. 

The overall reasonable agreement, as well as the self-consistency of the structure functions tested by or 
derived from heavy flavour production at HERA, enhances confidence in corresponding cross-section predic- 
tions at LHC, within their respective theoretical uncertainties. 
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Abstract 

We review selected aspects of the experimental techniques being prepared to study 
heavy flavour production in the four LHC experiments (ALICE, ATLAS, CMS and 
LHCb) and we present the expected performance for some of the most significative 
measurements. 

coordinators: A. Dainese, M. Smizanska, and C. Weiser 

1 Introduction 

author: M. Smizanska 

Unprecedently large cross sections are expected for heavy-flavour production in proton-proton collisions 
at the LHC energy of ^/s = 14 TeV. Next-to-leading order perturbative QCD calculations predict values of 
about 10 mb for charm and 0.5 mb for beauty, with a theoretical uncertainty of a factor 2-3. Despite these 
large cross sections, the LHC experiments, ALICE [1, 2], ATLAS [3], CMS [4], and LHCb [5], will have to 
deal with rejection of background coming from non-heavy flavour inelastic interactions for which the predicted 
cross is about 70 mb. The four experiments will work at different luminosity conditions. ATLAS and CMS 
are designed to work in a wide range of luminosities up to nominal 10^^ cm~^s~^, while the LHCb optimal 
luminosity will vary in the range (2-5) x 10^^ cm~^s~^ and ALICE is designed to work at 3 x 10'^'' cm~^s~^ 
in proton-proton collisions. Luminosity conditions in ATLAS, CMS and LHCb allow multiple interactions 
per bunch crossing, thus leading to requirements of even stronger identification and selection of heavy-flavour 
events already at trigger level. The first task will be the measurement of integrated and differential charm 
and beauty production cross sections in the new energy domain covered at the LHC. ALICE will play an 
important role, having acceptance down to very low transverse momentum, as we discuss in Section 4. These 
measurements can be performed within a relatively short period of running. Afterwards, the heavy-flavour 
studies will focus on less inclusive measurements addressing specific production mechanisms that allow to test 
higher order perturbative QCD predictions, as well as on rare decays of heavy-flavour hadrons, that may carry 
information on New Physics beyond the Standard Model. In order to meet these requirements dedicated and 
sophisticated trigger strategies have been prepared by the LHC experiments. 

2 Heavy flavour detection in the LHC experiments 

author: A. Dainese 

The four detectors that will take data at the LHC have different features and design requirements, but 
all of them are expected to have excellent capabilities for heavy-flavour measurements. Their complementarity 
will provide a very broad coverage in terms of phase-space, decay channels and observables. 

Experimentally, the key elements for a rich heavy-flavour program are track and vertex reconstruction 
and particle identification (PID). 

Open charm and beauty mesons have typical life-times in the order of a ps (ct values are about 125- 
300 /xm for D mesons and 500 /xm for B mesons) and the most direct detection strategy is the identification 
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Fig. 1 : Track impact parameter resolutions for thie four LHC experiments. Note that for ALICE, ATLAS and CMS the impact parameter 
is defined in the r<j> plane, while for LHCb it is defined in the rz plane. 



of single tracks or vertices that are displaced from the interaction vertex. The detector capability to perform 
this task can be characterized by the transverse impact parameter' (do) resolution. All experiments will be 
equipped with high position-resolution silicon-detector layers, including pixel detector for the innermost layers, 
for precise tracking and impact parameter measurement. Tracking is done in the central (pseudo)rapidity region 
for ALICE < 0.9), ATLAS and CMS (|r/| < 2.5), and in the forward region for LHCb (2 < ?? < 5). In Fig. 1 
we show the dp resolution, which is similar for the different experiments, and better than 50 /im for ~ 1-5- 
3 GeV. The inner detector systems of ATLAS, CMS and ALICE will operate in different magnetic fields: The 
ALICE magnetic field will vary within low values (0.2-0.5 T) leading to a very low px cutoff of 0. 1-0.2 GeV, 
while ATLAS (2 T) and CMS (4 T) have higher cutoffs of 0.5 and 1 GeV, respectively, but better px resolution 
at high pT (e.g., at pr = 100 GeV, bpr/pT ~ 1-2% for ATLAS/CMS at central rapidity and « 9% for 
ALICE). 

Both lepton and hadron identification are important for heavy-flavour detection. D and B mesons have 
relatively large branching ratios (BR) in the semi-leptonic channels, ~ 10% to electrons and ~ 10% to muons, 
and inclusive cross-section measurements can be performed via single leptons or di-leptons. Alternatively, 
high-pr leptons can be used as trigger-level tags to select B ^ J/^ + X candidate events, that provide more 
accurate cross section measurements. ALICE can identify electrons with > 1 GeV and \ri\ < 0.9, via 
transition radiation and dE/dx measurements, and muons in the forward region, 2.5 < < 4, which allows a 
very low pT cutoff of 1 GeV. CMS and ATLAS have a broad pseudorapidity coverage for muons, \r]\ < 2.4 and 
1 7? I < 2.7, respectively, but they have a higher pT cutoff varying between 4 and 6 GeV, depending on r]. Both 
CMS and ATLAS have high-resolution electro-magnetic calorimeters that will be used to identify electrons. 
Semi-leptonic inclusive measurements do not provide direct information on the D(B)-meson px distribution, 
especially at low px, because of the weak correlation between the lepton and the meson momenta. Therefore, 
for charm in particular, the reconstruction of exclusive (hadronic) decays is preferable. In this case, hadron 
identification allows a more effective rejection of the combinatorial background in the \ow-px region. ALICE 
disposes of vr/K/p separation via dE/dx and time-of-flight measurements for p < 3-4 GeV and \r]\ < 0.9. 

Figure 2 shows schematically the px vs. r] acceptances for charm (c) and beauty (b) hadrons in the four 
experiments, as expected for one year of running at nominal luminosity (note that the value of the luminosity is 

'We define as impact parameter the distance of closest approach to the interaction vertex of the track projection in the plane 
transverse to the beam direction. 
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1 year, pp, Vs = 14 TeV @ nominal lumin. 
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Fig. 2: Schematic acceptances in transverse momentum and pseudorapidity for open heavy flavour hadrons (indicated as 'Q-hadrons') 
in the four LHC experiments. The high-pT coverages correspond to one year (i.e. 7 months) of running at nominal luminosity (see 
beginning of this section). 



different for each experiment, as previously discussed). ATLAS and CMS have similar acceptances for beauty 
measurements; the minimum accessible px is relatively large because of the strong magnetic fields, which in 
turn, together with the high luminosity, allow to cover transverse momenta up to 200-300 GeV. The acceptance 
of LHCb, although centred at forward rapidity, has a significant overlap, with those of ATLAS and CMS. The 
acceptance of ALICE for beauty overlaps with ATLAS and CMS at central rapidity and with LHCb at forward 
rapidity. The moderate magnetic field allows measurements down to transverse momenta of about 2 GeV for 
B mesons in the forward muon arm and in the barrel, and down to about 1 GeV for D mesons in the barrel. 

3 Beauty triggers at the LHC 
3.1 ATLAS beauty trigger 

author: J. Baines 

The ATLAS trigger consists of three levels [6]. Level- 1 is implemented in hardware, whilst the higher 
level triggers (level-2 and the Event Filter, EE) are based on general-purpose processors. The level- 1 triggers 
are based on information from the calorimeter and muon trigger chambers. At higher trigger levels, information 
from the Inner Detector (ID) and precision muon detector is included. The size of the level-2 and EE processor 
farms is limited, which in turn limits the amount of data processing that can be performed in the trigger. The 
B -trigger must, therefore, have the flexibility to adapt selections both as the luminosity falls during a beam- 
coast and, over a longer time-scale, as the peak luminosity of the LHC increases. This is achieved by using a 
di-muon trigger at the start of higher luminosity fills and introducing additional triggers for lower luminosity 
fills or as the luminosity falls during a beam coast [7]. 

A di-muon trigger provides a very effective selection for a range of important channels, e.g. B[| ^ J /ip{^'^fi 
B K^*ij,fi and B — > p^fifi. The Level- 1 muon trigger is efficient down to a pT of about 5 GeV in the barrel 
region and about 3 GeV in the end-caps. However the actual thresholds used for the di-muon trigger will be 
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determined by rate limitations. For example, a pt threshold of 6 GeV would give a di-muon trigger rate of 
about 600 Hz after level-1 at a luminosity of 2 x 10^'^ cm^s~^. These triggers are mostly due to muons from 
heavy flavour decays plus some single muons which are doubly counted due to overlaps in the end-cap trigger 
chambers. The later are removed when the muons are confirmed at level-2 using muon precision chambers 
and ID information from inside the level-1 Region of Interest (Rol). At the EF tracks are refit, inside regions 
identified by level-2, and specific selections made on the basis of mass and decay length cuts. These consist of 
semi-inclusive selections, for example to select J/V'(/U+/i^) decays with a displaced vertex, and in some cases 
exclusive selections such as for B ^ The final trigger rate, after the EF, is about 20 Hz at a luminosity 

of 2 X 10^3 cm-^s-^ 

At lower luminosities, additional triggers are introduced based on a single muon trigger ( pT ~ 8 GeV) 
together with a calorimeter trigger. The calorimeter trigger identifies clusters of energy deposition in the elec- 
tromagnetic and hadronic calorimeter consistent with an electron or photon (EM Rol) or a jet (Jet Rol). For 
hadronic final states, such as Bg Dg^Tr"'" and Bg D^a^ track are reconstructed in the Inner Detector in 
Rol of about Ary x A(/) = 1.0 x 1.5. By limiting track reconstruction to the part of the ID lying within the 
Rol, about 10% on average, there is potential for up to a factor of ten saving in execution time compared to 
reconstruction in the full Inner Detector. Preliminary studies of efficiency and jet-cluster multiplicity have been 
made using a fast simulation which includes a detailed parameterization of the calorimeter. These studies in- 
dicate that a threshold on the jet cluster energy of Et > 5 GeV gives a reasonable multiplicity, i.e. a mean 
of about two Rol per event for events containing a muon trigger. This threshold would give a trigger that is 
efficient for Bg D~7r"'" events with a S-hadron px above about 15 GeV. 

Track reconstruction inside e/gamma Rol can be used to select channels such as B^J J /■ijj{e'^e^)K^, 
Bd K'^*7and Bg (pj. Preliminary studies show that a reasonable compromise between Rol multiplicity 
and electron efficiency might be obtained with a cluster energy threshold of Et > 2 GeV. This gives a mean 
Rol multiplicity of about one for events containing a muon trigger and is efficient for channels containing 
an electron with > 5 GeV. Following the ID track reconstruction further selections are made for specific 
channels of interest. These are kept as inclusive as possible at level-2 with some more exclusive selections at 
tiieEF. 

In LHC running, there will be competing demands for resources in the level-2 and EF trigger farms and 
for trigger band-width. By adopting a flexible strategy and making the maximum use of Rol information to 
guide reconstruction at level-2 and the EF, the B -physics coverage of ATLAS can be maximized. 

3.2 CMS beauty trigger 
author: R. Ranieri 

The Large Hadron Collider (LHC) will provide 40 MHz proton-proton colhsions at the centre of mass 
energy of 14 TeV. At the beginning a luminosity of 2x 10"^^ cm~^s~^ is expected, corresponding to 20 fb~^ col- 
lected per year. Assuming the bb production cross section to be 0.5 mb, 10^'^ b-physics events per year are 
foreseen: all kind of b-particles will be produced and studies will be performed not only in Bj^, but also in Bg 
meson system. A wide b-physics programme, including CP violation, Bg — Bg mixing and rare decays can 
therefore be covered by the CMS experiment. The apparatus will be equipped with a very precise tracking 
system made with silicon microstrip and pixel detectors [8,9]. 

The rate at which events can be archived for offline analyses is 100 Hz [10, 11]. The trigger thresholds are 
optimized for a wide physics discovery program with selection of high transverse momentum (px) processes. 
Low-pT events, as required for b-physics , are selected mainly by the first level muon trigger, then an exclusive 
reconstruction of few relevant benchmark channels can separate interesting events from the background. The 
b-physics programme could evolve with time following both the theoretical developments and the results which 
will be obtained in the next years by b-factories and Tevatron experiments. 

The lowest trigger level (Level-1) is based on the fast response of calorimeters and muon stations with 
coarse granularity. No information on secondary vertices is available, hence the Level-1 selection of b-physics 
events exploits the leptonic signatures from beauty hadron decays, therefore a single muon or a di-muon pair 
is required. The Level-1 output at start-up will be 50 kHz. Several studies have been done to optimize the 
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trigger thresholds in order to have the possibiUty of selecting most of the interesting physics signatures. A 
total of 3.6 kHz rate is dedicated to the Level-1 muon selection. It is obtained by requiring a single muon with 
Pt > 14 GeV or at least two muons with > 3 GeV. 

A further selection is made during the High-Level trigger (HLT) by using also the information from the 
tracking system. The CMS High-Level trigger is entirely based on a CPU farm with some thousand CPUs. 
Each processor analyses a single event; in principle offline event reconstruction can be performed, but in order 
to reduce the processing time fast track reconstruction has to be done. Some algorithms will be dedicated to 
the fast reconstruction and identification of physics processes, thus allowing to start the offline analysis directly 
from the online selection. They have to fulfill the HLT time constraint, hence they have to be able to analyze 
and accept (or reject) data within the time limits imposed by the HLT latency. To lower the execution time, 
which is due mainly to the processing of tracking system signals, track reconstruction is preferably performed 
only in limited regions of the space (regional track reconstruction) and stopped when a certain precision is 
reached in the measurement of some track parameters, such as transverse momentum and impact parameter 
(conditional track finding). Invariant mass of b-hadrons can thus be measured online with good resolution, 
allowing to select the searched event topologies. 

An additional trigger strategy, which relies on the possibility of lowering the trigger thresholds during 
the LHC beam coast or lower luminosity fills to collect more b-physics events is under study. 

The rare decay ^ ^ A'^y"" is triggered at Level-1 with 15.2% efficiency. At HLT, the two muons are 
required to be opposite charged and isolated, to come from a displaced common vertex and have an invariant 
mass within 150 MeV from the Bg mass. The estimated background rate is below 2 Hz and nearly 50 signal 
events are expected with 10 fb~^. 

The determination of Am^ and AF^, will be a valuable input for flavour dynamics in the Standard Model 
and its possible extensions. The measurement of Arug is allowed by the Bg^Ds~7r"'" decay followed by 
Ds~ — ^^7r~ and ^— >^K+K~. The Bg CP state at decay time is tagged by the charge of the pion associated 
to the Dg (in this case the 7r+). The only way to trigger on these hadronic events is to search for the muon 
coming from the decay of the other b quark in the event. In addition to the single muon Level-1 trigger, it was 
studied the possibility of a combined trigger with a low-pT muon and a soft jet. The CMS High-Level trigger 
algorithm reconstructs the charged particle tracks with only three points by using the precise pixel detector. 
Topological and kinematical cuts are applied to reconstruct the three resonances (p, Dg and B°. A 20 Hz output 
rate is achieved with about 1000 signal events in 20 fb^^. Since the overall possible rate on tape is 100 Hz, the 
bandwidth allocated to this channel probably could not exceed 5 Hz. If the fraction of events written to tape 
is scaled accordingly, more than 300 signal events are expected for 20 fb~^. In order to fully cover the range 
allowed by the Standard Model, about 1000 events are needed. 

The decay channel Bg ^ J /iptp is very important because it can not be studied with large accuracy before 
LHC and can reveal hints for physics beyond the Standard Model. Events with a couple of muons are passed 
to the HLT. The inclusive selection of J/V'^A*^M~ decays, obtained with mass requirements on the di-muon 
system, leads to a total of 15 Hz rate, 90% of which is made of J/^ from b quarks. With an additional amount of 
CPU time, perhaps sustainable by the HLT computing power, about 170 000 events are expected in 20 fb~^ with 
less than 2 Hz rate. 

3.3 LHCb beauty trigger 

author: H. Ruiz 

The LHCb detector [5] is optimized for exploiting the B-physics potential of LHC. Together with excel- 
lent vertexing and particle identification, an efficient trigger on a wide variety of B decays is one of the main 
design requirements of the experiment. 

The LHCb trigger system [12] is organized in three levels. The first one (LO) runs on custom electronics 
and operates synchronously at 40 MHz, with a 4 //s latency. The remaining two trigger levels (LI and HLT) run 
on a shared farm of 1400 commercial CPUs. A brief description of the three trigger levels and their performance 
follows. 

LO exploits the relatively high px of B decay products. High px candidates are identified both in the 
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calorimeter and in the muon system, with thresholds of about 3 and 1 GeV respectively. Complicated 
events that would consume unreasonable time at higher levels are promptly vetoed in two different ways. First, 
multiple primary vertex topologies are rejected by using two dedicated silicon layers of the vertex detector. 
Secondly, events with large multiplicity, measured at a scintillating pad layer, are vetoed. The input rate of 
events visible in the detector is about 10 MHz, with a bb content of 1.5%. LO reduces this rate by a factor of 10 
while increasing the bb content to 3%. The typical efficiency of LO is 90% for channels with dimuons, 70% for 
radiative decays and 50% for purely hadronic decays. 

At the 1 MHz input rate of LI it becomes feasible to use tracking information, allowing for the search 
for B vertex displacement signatures. Tracks are first searched at the vertex detector, and then confirmed in 
two dedicated tracking layers (trigger tracker or TT) which provide a rough estimation of the momentum of 
the tracks (Spt/pt ~ 25%). The generic LI decision is based on the presence of two tracks with an impact 
parameter higher than 0.15 mm with a sufficiently high value of log{pTi + Pt2)- Alternative selection criteria 
are applied, based on the presence of tracks matched to LO neutral calorimeter objects and muon candidates. 
The output rate of LI is 40 kHz with a bb content of 15%. The efficiencies are at the level of 90, 80 and 70% 
for channels with di-muons, only hadrons and radiative decays respectively. 

The HLT [13] consists of two sequential layers. The first one refines the LI decision with the all the 
tracking information from the detector, improving the pT measurement to the level of 5pt/pt ^ 1%- The 
rate is reduced to 13 kHz and the bb content is enriched to 30%. The second layer consists on a series of 
alternative selections. A first group of them aims for maximal efficiency on the base-line physics channels and 
the corresponding control samples, by making use of the complete reconstruction of the decay vertex and its 
kinematical properties. These selections fill 200 Hz of bandwidth, while providing efficiencies typically higher 
than 90%. The rest of selections aim for more generic signatures that will provide robustness and flexibility to 
the trigger system. In addition, the samples selected will be useful for calibration and systematic studies. The 
selections aim for generic J/ip and D* (600 Hz and 300 Hz respectively), and generic B decays (900 Hz). The 
latter is based on the detection of single muons with high pT and impact parameter. 

In total, 2 kHz of events will be written on tape, with an expected overall efficiency ranging between 
75% for channels with di-muons to 35% for purely hadronic final states. 

4 Measurements in preparation at the LHC and expected performance 

In the following we present, as examples, the expected performance for the detection of D and B mesons in 
ALICE^, and for the study of bb azimuthal correlation in ATLAS. We also include a summary of the capability 
of ALICE of quarkonia measurements {^p family and T family). 

4.1 Charm reconstruction in ALICE 

author: A. Dainese 

One of the most promising channels for open charm detection is the D° K^tt"*" decay (and charge 
conjugate) that has a BR of 3.8%. The expected yields (BR x dN/ dy at y = 0), in pp collisions at y/s = 14 TeV 
and in central Pb-Pb (0-5% cr*°*) at ^fs^ = 5.5 TcV are 7.5 x 10"'* and 5.3 x 10"^ per event, respectively [14]. 

The main feature of this decay topology is the presence of two tracks with impact parameters do ~ 
100 /im. The detection strategy [15] to cope with the large combinatorial background from the underlying 
event is based on the selection of displaced-vertex topologies, i.e. two tracks with large impact parameters and 
good alignment between the D° momentum and flight-line, and on invariant-mass analysis to extract the signal 
yield. This strategy was optimized separately for pp and Pb-Pb collisions, as a function of the transverse 
momentum, and statistical and systematic errors were estimated [16, 17]. The results, in terms of pt coverage 
and statistical precision, are found to be similar for the two colliding systems [16, 17]. 

Figure 3 (left) shows the expected sensitivity of ALICE for the measurement of the D° p^-differential 

^Given that ALICE is dedicated to the study of nucleus-nucleus collisions at the LHC, some of the presented results are relative to 
Pb-Pb collisions at i/s = 5.5 TeV per nucleon-nucleon coUisions. These results can be taken as lower limits for the performance in 
pp collisions, where the background contributions are much lower. 
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Fig. 3: ALICE sensitivity on d^a^^ /dprdy, in pp at 14 TeV, compared to tlie pQCD predictions obtained witli different sets of tfie 
input parameters rric [GeV], jip/l-io, I-Ir/i-Io and PDF set (fio is defined in the text). The inner bars represent the statistical error, the 
outer bars the quadratic sum of statistical and pT-dependent systematic errors. A normalization error of 5% is not shown. The panel on 
the right shows the corresponding 'data/theory' plot. 



cross section in pp collisions, along with NLO pQCD [18] calculation results corresponding to different choices 
of the charm quark mass and of renormalization and factorization scales. In the right-hand panel of the figure 
we present the ratio 'data/theory' ('default parameters/theory parameters') which better allows to compare the 
different pj-shapes obtained by changing the input 'theory parameters' and to illustrate the expected sensitivity 
of the ALICE measurement. The estimated experimental errors are much smaller than the theoretical uncer- 
tainty band. We note that the data cover the region at low transverse momentum where the accuracy of the 
pQCD calculation becomes poorer and where novel effects, determined by the high partonic density of the ini- 
tial state at LHC energies, may play an important role (see "Small-x effects in heavy quark production" section 
of this report). 

4.2 Beauty production measurements in ALICE 

author: A. Dainese 

The expected yields (BR x dN/dy at y = 0) for B ^ + X plus B ^ D e=^ + X) + X' in 
pp colhsions at ^ = 14 TeV and in central Pb-Pb (0-5% cr*°*) at ^/s^ = 5.5 TeV are 2.8 x 10"^ and 
1.8 X 10^^ per event, respectively [14]. 

The main sources of background electrons are: (a) decays of D mesons; (b) neutral pion Dalitz decays 
vr'^ — > 7e"'"e~ and decays of light mesons (e.g. p and uj); (c) conversions of photons in the beam pipe or in the 
inner detector layers and (d) pions misidentified as electrons. Given that electrons from beauty have average 
impact parameter c/q — 500 /im and a hard momentum spectrum, it is possible to obtain a high-purity sample 
with a strategy that relies on: electron identification with a combined dE/dx and transition radiation selection, 
which allows to reduce the pion contamination by a factor 10*; impact parameter cut to reject misidentified 
pions and electrons from sources (b) and (c); transverse momentum cut to reject electrons from charm decays. 
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Fig. 4: B production cross section vs. p™'" reconstructed by ALICE in 10^ central Pb-Pb events. Left: da^j^/dy at j/ = (normalized 
to one nucleon-nucleon collision) from single electrons in j?;! < 0.9; statistical (inner bars) and quadratic sum of statistical and pr- 
dependent systematic errors (outer bars) are shown; a 9% normalization error is not shown. Right: cr^ integrated in 2.5 < < 4 (not 
normalized to one nucleon-nucleon collision) from single muons and di-muons in 2.5 < rj < 4; only (very small) statistical errors 
shown. 



As an example, with do > 200 ^um and pT > 2 GeV, the expected statistics of electrons from B decays is 
8 X 10*^ for 10^ central Pb-Pb events, allowing the measurement of electron-level -differential cross section 
in the range 2 < pT < 18 GeV. The residual contamination of about 10%, acculated in the low-pT region, 
of electrons from prompt charm decays and from misidentified charged pions can be evaluated and subtracted 
using a Monte Carlo simulation tuned to reproduce the measured cross sections for pions and D° mesons. A 
Monte-Carlo-based procedure can then be used to compute, from the electron-level cross section, the B-level 
cross section da^{pT > p™™)/dy [17]. In the left-hand panel of Fig. 4 we show this cross section for central 
Pb-Pb collisions with the estimated uncertainties. The covered range is 2 < p™™ < 30 GeV. 

B production can be measured also in the ALICE forward muon spectrometer, 2.5 < 77 < 4, analyzing 
the single-muon pT distribution and the opposite-sign di-muons invariant mass distribution [17]. 

The main backgrounds to the 'beauty muon' signal are vr^, and charm decays. The cut pT > 
1.5 GeV is apphed to all reconstructed muons in order to increase the signal-to-background ratio. For the 
opposite-sign di-muons, the residual combinatorial background is subtracted using the technique of event- 
mixing and the resulting distribution is subdivided into two samples: the low-mass region, M^+^- < 5 GeV, 
dominated by muons originating from a single b quark decay through b c(-^ /^^)/^~ (BDsamc)» and the 
high-mass region, 5 < M^+^- < 20 GeV, dominated by bb — > fi'^lJ.~, with each muon coming from a 
different quark in the pair (BBdig). Both samples have a background from cc ijl^ ijl^ and a fit is done 
to extract the charm- and beauty-component yields. The single-muon pT distribution has three components 
with different slopes: K and vr, charm, and beauty decays. Also in this case a fit technique allows to extract 
a pt distribution of muons from B decays. A Monte Carlo procedure, similar to that used for semi-electronic 
decays, allows to extract B-Ievel cross sections for the data sets (low-mass //^//^, high-mass //^//^, and pT- 
biimed single-muon distribution), each set covering a specific B-meson pT > p^™ region, as preliminarly 
shown in Fig. 4 (right). Since only minimal cuts are applied, the reported statistical errors are very small and 
high-pT reach is excellent. Systematic errors are currently under study. 

4.3 Study of bb correlations in ATLAS 

author: Th. Lagouri 
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Fig. 5: Distribution of the opening angle Ac^(J/t/)-/^) between the }/ip from the decay Bd J/ipK" and the muon coming from the 
associated B hadron decay to muon, both direct b — » /tt and indirect b — » c — > /i. 



The ATLAS detector [19] is well engineered for studies of b-production, and together with the huge 
rate of b-quark production that wiU be seen at LHC, offers great potential for the making of novel precise b 
production measurements. Correlations between b and b quarks and events with more than one heavy-quark 
pair, bbbb, bbcc, bbss, that were difficult to access in previous experiments due to limited statistics, will be 
investigated in detail. A new technique has been developed in ATLAS for measuring correlations, and this will 
yield results that will shed new hght on our understanding of the QCD cross-section for bb-production. 

A detailed study investigated a possibility of bb correlations measurement using the A0(J/i/'-/x) distri- 
bution, the azimuthal separation of a J/ip and a muon [20, 21, 22]. This technique is expected to be superior to 
earlier methods used at the Tevatron Run-1 based on muon-muon or muon-b-jet correlations. The new method 
does not require separation cuts between the two objects. Such cuts were necessary to control the background, 
but they required a model-dependent extrapolation of the results to full azimuthal space [23]. Using a full sim- 
ulation of the Inner Detector and the Muon Spectrometer of the ATLAS detector [19] it is shown that such a 
distribution can be extracted from heavy flavour events at LHC. 

ATLAS studies were done for two channels selected to measure the azimuthal angle difference A0(bb) 
between b and b quarks: 

b^ Bd ^ J/V'H ^//)K°, h^n + X 

h^Bs^ 3 h^/i + X 

The numbers of events expected for 30 fb~^ as might be achieved after 3 years of running at a luminosity 
of 10^^ cm^s~^ are 4.8x10^ and 3.2x10^ respectively for these channels. No isolation cuts are needed to 
separate exclusively reconstructed B -decays from the muon produced in the semi-leptonic decay of the other 
B -particle in the event. The reconstruction efficiency remains high in topologies where the azimuthal angle 
difference A<p(J/ip-^) between and the muon is small. 

Special attention was devoted to background events in which the muon is produced from the decays K='=, 
tt"^ + X instead of b^ /i + X. The study showed that this background is not problematic in B^ decays, 

however it is important in the case of Bg meson. 

In summary, the results of the analysis suggest that backgrounds from K/tt decays are small, and that 
backgrounds from events containing 4 b quarks are relatively flat in A^(J/^/;-/li). The efficiency of the recon- 
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struction of muons with this technique is also relatively flat in A^(J/^-/u) and so we conclude that corrections 
to the measured A(j){J/ip-n) distribution are likely to be small. 



4.4 Quarkonia measurements in ALICE 

author: A. Morsch 

Heavy quarkonia states are hard, penetrating probes which provide an essential tool to study the earliest 
and hottest stages of heavy-ion collisions [24]. They can probe the strongly interacting matter created in these 




Fig. 6: Acceptance for J/?/) and T as a function of y and px for measurements in the di-muon channel and di-electron channels. To give 
an idea of the effect of the trigger, the acceptances are shown without (solid) and with (dashed) a sharp cut on the transverse momentum 
of single muons of 1 GeV/c (2 GeV/c) for J/?/) (T ) and for single electrons of 3 GeV. 



49 



reactions on short distance scales and are expected to be sensitive to the nature of the medium, i.e. confined or 
de-confined [25, 26]. The suppression (dissociation) of the heavy-quark resonances is considered as one of the 
most important observables for the study of the QGP at the LHC (see Ref. [2] for a recent review). 

In ALICE, quarkonia will be measured in the di-electron channel using a barrel (|ry| < 0.9) Transition 
Radiation Detector (TRD) [2] and in the di-muon channel using a forward Muon Spectrometer (2.5 < r/ < 
4) [2]. The complete spectrum of heavy-quark vector mesons (J/^, ijj', T, T', T") can be measured down to 
zero pt- In particular the good mass resolution allows to resolve the Upsilon family. 

The Muon Spectrometer uses a low-pT trigger threshold, pT > I GcV, on single muons for charmonia 
and a high-pr trigger, px > 2 GeV, for bottonia detection. The TRD can trigger on single electrons with 
Pt > ^ GeV, which results in a minimum transverse momentum of triggered charmonia of 5.2 GeV. Electron 
identification combined with the excellent vertexing capabilities of the inner tracking system allows ALICE to 
distinguish direct charmonium production from secondary charmonium production through B decays. 

The energy density dependence will be studied by varying the impact parameters and by studying in 
addition to the heaviest collision system (Pb-Pb) also intermediate mass and low mass A-A systems. To 
determine the primary production cross-section of the resonances and the amount of pre-resonance absorption, 
corresponding measurement have to be performed for pA and pp collisions. 

Table 1: Main characteristics of quarkonia detection witii ttie TRD and tlie Muon Spectrometer in ALICE. 





Muon Spectrometer 


TRD 


Acceptance 

Mass Resolution J/^ 

Mass Resolution T 


2.5 < r/ < 4 
72 MeV 
99 MeV 


|?7| < 0.9 
34 MeV 
93 MeV 



Table I shows the main quarkonia detection characteristics of the TRD and the Muon Spectrometer and 

the acceptances for J/^|) and T as a function of y and px are shown in Fig. 4.4. In one year of pp running at 
{L) = 3 X 10^'^ cm~^s~^ ALICE will detect several 10® J/tps and several 10^ Ts in the di-muon channel. A 
T statistics of 10^-10^ can be obtained in the di-electron channel. For the J/tp we expect « 10^ untriggered 
low-pT and 10^ high-p^ triggered events. 
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Abstract 

We study small-x effects on heavy flavor production at the LHC in two approaches in- 
cluding nonlinear, saturation-motivated, terms in the parton evolution. One approach 
is based on coUinear factorization, the other on kr factorization. The prospects for 
direct experimental study of these effects in pp collisions at the LHC are discussed. 

coordinators: A. Dainese, H. Jung, and R. Vogt 

1 Introduction 

HERA data are used to constrain the small x, moderate parton densities in two approaches. In the first, 
HERA F2 data are refit using DGLAP evolution with the first nonlinear recombination terms. Recombination 
slows the evolution so that, after refitting the data, the gluon distribution is enhanced relative to that obtained by 
DGLAP alone. The resulting set of parton densities produces charm enhancement in pp collisions at the LHC. 
On the other hand, assuming kx factorization, the unintegrated gluon distribution is determined from the HERA 
F2 data, the only inclusive HERA measurements which directly accesses the gluon density. Saturation effects 
are then included, reducing the small x gluon densities with little distinguishable effect on F2. This approach 
leads instead to heavy flavor suppression at the LHC. After a short general introduction, both approaches and 
their predicted effects on heavy quark production are discussed in detail. Direct experimental study of these 
effects in pp collisions at the LHC may be able to differentiate between the two approaches. 

2 Small-x partons, absorption and the LHC 

authors: A.D. Martin and M.G. Ryskin 

2.1 Partons densities at low xl 

Almost nothing is known about the behaviour of partons at low x. There are essentially no data available for 
X < 10~^ with in the perturbative region and there is no reliable theory to extrapolate down in x. 

In the Dokshitzer-Gribov-Lipatov-AItarelU-Parisi (DGLAP)-based [1, 2, 3, 4] global analyses, the small- 
X behaviour is driven by input distributions at a starting scale Q = Qo- Usually these 'input' distributions are 
written in the form x~^(l — x)^ where A and i] are free parameters fit to the data. So one can say nothing 
without data in the x region of interest. Moreover, there may be large low-x contributions to the gluon of 
the form (a^ ln(l/x))" - the so-called Balitsky-Fadin-Kuraev-Lipatov (BFKL) effects [5, 6, 7, 8], beyond the 
DGLAP approximation. 
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Fig. 1: The gluon-gluon splitting function, Pgg, with fixed and running coupling in the LL and resummed NLL BFKL approximations, 
compared with the LO DGLAP behaviour. The figure is taken from Ciafaloni et al. [9, 10, 11, 12, 13, 14]. The subscript B refers to 
scheme B which ensures energy-momentum conservation in the splitting. 

Thus it seems better to discuss low-x behaviour in terms of BFKL-evolution. However there are also 
problems here. The next-to-leading logarithm (NLL) corrections to the leading order (LO) BFKL (CCFM) 
amplitude are known to be very large and one needs to resum such corrections to obtain a relatively stable result. 
We cannot justify the perturbative QCD approach at low so that the solution of the BFKL equation need 
to be matched to some non-perturbative amplitude at Q = Qo- This non-perturbative distribution (analogous 
to the 'input' in the DGLAP case) is not known theoretically. Either it has to be fit to low x data or some 
phenomenological model such as a Regge-based parametrization has to be used. 

After a reasonable resummation of the NLL corrections in the region where the starting virtuality Qq is 
not close to the final value of Q,Q > Qq, the resummed BFKL amplitude turns out to be similar to that resulting 
from DGLAP evolution [9, 10, 11, 12, 13, 14]. For example, the preasymptotic effects on the resummation of 
the gluon-gluon splitting function are so large that the NLL BFKL power growth only sets in for z < 10^^, as 
can be seen from Fig. 1. Moreover, the behaviour of the convolution Pgg ® g/g, normalized to g, in the NNLO 
DGLAP and NLL approximations is exactly the same down to z ~ 10~^ [15]. 

Thus, in practice, the DGLAP and BFKL based approaches are rather close to each other in the HERA 
kinematic regime. In both cases, the main problem is the low-x behaviour of the amplitude at Q = Qo where 
we need to phenomenologically determine possible non-perturbative contributions, power corrections and so 
on. 

2.2 The puzzle of the x"^ behaviour 

Since the BFKL amplitude grows as a power oi x, A oi x"^, it will violate unitarity as x — > 0. Indeed, 
even after the NLL resummation, the expected power, A ~ 0.3, is rather large. Thus, we first discuss absorption 
effects which tame the violation of unitarity. The upper limit of the small x behaviour of the parton distributions 
a = (7, g is given by the extrapolation 



below xq = 0.001. The distributions are reliably determined from global parton analyses at x > xq. 

On the other hand, it is reasonable to expect that at Q ^ Qo ~ 1 GeV the behaviour will reflect that 
of hadron-hadron interactions: A = 0.08 [16]. Most likely the lower value of A is explained by absorp- 
tive/screening effects in hadron-hadron collisions. So, for extrapolation down to x ~ 10^"^ — 10^^ we may 




(1) 
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Fig. 2: The behaviour of the gluon and sea quark distributions at = 2 GeV^ found in the CTEQ6.1M [17] and MRST2004 NLO [18] 
global analyses. The valence-like behaviour of the gluon is evident. 



regard A = 0.08 as a lower limit since, in DIS, we expect smaller absorptive effects than those in hadron-hadron 
interactions. 

However, present global analyses, which do not allow for absorption effects, reveal that at Q ~ 1 — 1.5 
GeV and low x, the sea quarks have a Pomeron-like behaviour, xq ~ x~^''^, whereas the gluon distribution has 
a valence-like behaviour, xg ~ ^/x. This different behaviour is evident from Fig. 2, which shows the behaviour 
of the gluon and sea quark distributions, xS = 2x{u + J + s) for = 2 GeV^. Such a result looks strange 
from the Regge viewpoint where the same vacuum singularity (Pomeron) should drive both the sea quarks and 
the gluons since the same power is expected for sea quarks and gluons, \g = Xg. 

This difference demonstrates that the actual situation is even more complicated. It is worth noting that a 
simultaneous analysis of inclusive and diffractive DIS data indicates that, after accounting for screening effects 
and allowing for some power corrections, it may be possible to describe the HERA data with Xg = Xg = [19]. 
The absorptive effects, estimated from fitting the diffractive DIS data, enlarge the input gluon distribution at 
low x. 

It may initially seem strange that accounting for absorptive effects gives a larger gluon density^ at low 
X and Q^. The point is that the only way to describe the data, which are sensitive to absorptive effects, within 
the framework of DGLAP evolution without absorption, is to choose a very low 'input' gluon density in order 
to mimic the screening corrections 'hidden' in the data. Indeed, there is a tendency for the gluon distribution to 
even become negative at low x and Q^. On the other hand, allowing for absorption during DGLAP evolution 
(with the help of the Gribov-Levin-Ryskin (GLR) and Mueller-Qiu (MQ), GLRIVIQ, equations [22, 23]) the 
same data may be described with a larger and definitely positive input gluon density at Q = <3o- 



2.3 Estimates of absorptive effects: GLRMQ to BK 

The saturation of parton densities (A = 0) may be obtained using the Balitski-Kovchegov (BK) [24, 25] equa- 
tion, based on the BFKL equation, as well as the aforementioned GLRMQ equations. The latter equation is 
based on DGLAP evolution. These equations sum the set of fan diagrams which describe the rescattering of 
intermediate partons on the target nucleon. The screening caused by these rescatterings prohibits the power 
growth of the parton densities. 

The GLR equation for the gluon may be written symbolically as 

- Pgg <S) g + Pgq q - / — [y^i^jQ )] ■ (2) 



dlnQ^ ^ " ' '"^^ 16i?2Q2 J y 



'The same result was obtained in Ref. [20, 21] - note, however, it was based on LO evolution and the large NLO correction to Pqg 



will change the g, g relationship. 
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The nonlinear shadowing term, oc — [g]^, arises from perturbative QCD diagrams which couple four gluons to 
two gluons so that two gluon ladders recombine into a single gluon ladder. The minus sign occurs because 
the scattering amplitude corresponding to a gluon ladder is predominantly imaginary. The parameter ii is a 
measure of the transverse area ttR^ where the gluon density is sufficiently dense for recombination to occur. 

The BK equation is an improved version of the GLR equation. It accounts for the more precise triple- 
pomeron vertex (first calculated in Ref. [26, 27, 28]) and can be used for the non-forward ampUtude. The GLR 
equation was in momentum space, whereas the BK equation is written in coordinate space in terms of the dipole 
scattering amplitude A^(x, y, Y) = A'xy(l"). Here x and y are the transverse coordinates of the two t-channel 
gluons which form the colour-singlet dipole and Y = ln(l/x) is the rapidity. The BK equation reads 



For small dipole densities, A^, the quadratic term in the brackets may be neglected and Eq. (3) reproduces the 
conventional BFKL equation. However for large N, that is AT —s- 1, the right-hand side of Eq. (3) vanishes and 
we reach saturation when N = 1. The equation sums up the set of fan diagrams where at small Y the target 
emits any number of pomerons (i.e. Unear BFKL ampUtudes) while at large Y we have only one BFKL dipole. 

Starting from the same initial conditions, the solution of the BK equation gives /ewer small-x partons 
than that predicted by its parent linear BFKL/CCFM equation^. 

In principle, it would appear more appropriate to use the BFKL-based BK equation to describe the parton 
densities at low x. Unfortunately, however, the BK equation is only a model. It cannot be used for numerical 
predictions. We discuss the reasons below. 

2.4 Status of the BK equation 

The Balitski-Kovchegov (BK) equation [24, 25] is an attempt to describe saturation phenomena. However it is 
just a 'toy model' and cannot, at present, be used to reliably estimate absorptive effects at small x. The reasons 
are as follows: 

• The BK equation is based on the LO BFKL/CCFM equation. We know that the NLL corrections are 
large. We need to know the NLL corrections not only for the linear part of the evolution, but also for the 
nonlinear term. 

• Even neglecting the NLL corrections, we need to match the solution to a boundary condition at rather 
low Q^. This boundary condition is not theoretically known. 

• It sums a limited set of diagrams. The selection of diagrams (the fan graphs) was justified in the region 
where absorptive effects are relatively small. When these corrections become important, as in the satu- 
ration region, one has to allow for many other graphs whose contributions become comparable to the fan 
diagram contributions^ . 

• To solve the BK equation we need an initial condition at fixed x and all Q^. These conditions are not well 
enough known. In particular, the maximum (saturation) value of the gluon density depends on the radius: 

xg{x, q^) oc R^q^. At the moment, the radius i? is a free parameter. It may be small — the so-called 'hot 
spot' scenario. Moreover, we should account for the possibility of dissociation of the target proton^. The 
contribution coming from the dissociation is expected to have a smaller R. 

2.5 Relevance to, and of, the LHC 

How do the uncertainties at low x affect the predictions for the LHC? Fortunately for inclusive production 

of possible massive new particles with M > 100 GeV, the partons are sampled at x values and scales M 
reliably determined from NLO and NNLO global analyses. For illustration, we discuss W production which 

^Analogously, starting from the same input (and not fitting the input to the data) the GLR equation gives fewer small-a; partons than 
that predicted by the parent Hnear DGLAP equation. 

'Unfortunately the problem of summing all relevant diagrams has not been solved, even in the simpler case of Reggeon field theory. 
We know that these channels provide more than 30 - 40% of measured at HERA. 



dY 




(3) 
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has been studied in detail [29, 30, 31]. Central W production (yw = 0) at the LHC samples partons at 
X = Mw/\/s = 0.006. However to predict the total cross section, aw, we need to integrate over rapidity, 
important for \yw\ ^ 4 so that aw has some sensitivity to partons as low as x ~ 10~^. The total uncertainty 
on the NNLO prediction of aw has been estimated to be ±4% [29]. Therefore W production at the LHC can 
serve as a good luminosity monitor. To reduce the uncertainty in the prediction of aw will require a better 
theoretical understanding of low x partons. 

Of course, if the new particles are sufficiently massive, M > 1 TeV, and produced by gluon fusion, then 
the uncertainties due to the PDFs will be larger. However, there are situations where the scale is considerably 
lower such as exclusive double-diffractive Higgs production which depends on the unintegrated gluon at ~ 5 
GeV^ with x ~ Mh/\/s ~ 0.01. The absorptive effects are also expected to be small here. 

Turning the discussion around, is it possible for the LHC experiments to determine the behaviour of 
partons in the x region below 10~^ at low scales? One possibility is fj,'^fi~ Drell-Yan production in which 
events are observed with the //"'"/Lt" invariant mass as low as possible and the rapidity as large as possible. For 
example, for M^^ = 4 GeV and = 3, we sample quarks at a; = 1.4 x 10~^. This process predominantly 
samples the sea quark distributions. To study the small x behaviour of the gluon at low scales we may consider 
Xc production, or prompt photon production driven by the subprocess gq ^yq. 

In practice, rather than Xc, it may be better to study pp J ftp X as a function of yj/^. This process is 
also sensitive to the gluon distribution through the subprocesses gg ^ J/ip g, gg —>■ X ^ ^/'^ 1- There are 
also contributions from the subprocesses gg bb with b ^ J/ijj, and qq J/ip. The analysis of such data 
will be considerably helped by the detailed observations of prompt J/ijj and J/tp from b in central production 
at the Tevatron [32]. In fact, the first ever NLO global parton analysis [33] used J/ijj data as a function of 
rapidity to constrain the gluon distribution. 

The LHCb detector covers the rapidity region of 2 < r/ < 5 [34], and may be able to perform some 
of the above measurements. There is another possibility. Since LHCb will operate at a luminosity of 2 x 
10^^ cm~2s^^, there will usually be a single collision per bunch crossing and thus practically no 'pile-up' 
problems. Installing a forward detector at LHCb would offer the possibility of observing asymmetric events 
with one very large rapidity gap to probe the region of xjp < 10~^. 



3 Including nonlinear terms in gluon evolution: the GLRMQ and BK approaches 
3.1 GLRMQ approach 

authors: K.J. Eskola and VJ. Kolhinen 

The DGLAP [1, 2, 3, 4] evolution equations describe the scale evolution of the parton distribution func- 
tions (PDFs) well in the region of large interaction scale, ^ 4 GeV^ [35, 17, 36]. However, toward small 
values of x and Q"^, the gluon recombination terms start to play an increasingly important role. The inclusion 
of correction terms which arise from fusion of two gluon ladders leads to nonlinear power corrections to the 
DGLAP evolution equations. The first of these nonlinear corrections are the GLRMQ terms. 

With the GLRMQ corrections, the gluon evolution equation becomes 



dxg{x,Q^) dxg{x,Q'^ 



DGLAP 2 



dlnQ^ 51nQ2 
We model the two-gluon density in the latter term on the right-hand side as 



r %'G('Hy,Q'). (4) 

Jx y 



x^G^^\x,Q^) = -i^[x5(x,Q2)]2^ (5) 

where i? = 1 fm is the radius of the proton (we comment further on this later). The corrections to the sea quark 
distributions are 



dxq{x,Q'^) dxqix.Q"^^ 



2 



DGLAP 20 



ainQ2 d\^Q 
We have assumed that the higher-twist gluon term, Ght [23], is negUgible 
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Fig. 3: Calculated F2 (x, Q^) values compared with the HI data. 



Since these coirection terms are negative, they slow down the evolution of the PDFs. Due to the 1 /Q^ 
dependence, they also die out in the evolution so that at large scales Eqs. (4) and (6) relax into the linear DGLAP 
equations. 

In order to study the interplay between the nonlinear corrections and the initial PDFs and observe 
the nonlinear effects in fits to the DIS data, in Ref. [37] we compared the structure function F2{x,Q'^) = 
J2q '=^q[^Q{^j Q'^) + xq{x, Q'^)], calculated with the nonhnearly-evolved PDFs, to the HERA HI data [38]. As 
reference distributions we used the CTEQ5L and CTEQ6L PDF sets at large scales. We chose these sets be- 
cause the CTEQ collaboration uses only the large scale, > 4 GeV^, data in their fits, thus avoiding some of 
the possible nonlinear effects appearing in the small x,Q^<4 GeV^ region in their initial distributions. 

At small X, sea quarks dominate F2 and the gluon distribution dictates its scale evolution. At leading 
order (LO), the DGLAP contribution can be approximated as dF2{x, Q'^)/dlnQ'^ « {10as/27TT)xg{2x, Q"^) 
[39]. Larger xg{x, Q^) values correspond to faster F2{x, Q"^) evolution. The scale evolution of ^2(0;, Q"^) at 
small X computed with CTEQ5L is too fast due to a large small-x small-Q^ gluon distribution. The newer 
CTEQ6L set has much smaller gluon distribution in this region (see Fig. 4 (left)), giving a slower evolution and 
hence a good fit to the HI data. 

Our goal in Ref. [37] was to determine whether the good fit to the data could be maintained using 
the GLRMQ-corrected DGLAP scale evolution together with initial scale PDFs differing from CTEQ6L. We 
constructed a new set of PDFs using the CTEQ5L and CTEQ6L distributions piecewise as baselines at scales 
~ 3 — 10 GeV^ where the linear terms dominate the evolution and evolved them nonlinearly to lower Q^. 
We then interpolated between the sets in x and assumed a power-like dependence at small-x for gluons and 
sea quarks. These initial PDF candidates were then evolved to higher scales and compared to the data. This 
iterative procedure was repeated until a sufficiently good fit to the data was found. 
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Fig. 4: Left: initial gluon distributions at Qq = 1.4 GeV^. Right: evolution of gluon distributions for several fixed values of x shows 
that the effect of the nonlinear terms vanishes as increases. 



As a result, we obtained a new set of initial PDFs^, called EHKQS, which, when evolved using the 
nonlinear DGLAP+GLRMQ evolution equations, produced an equally good or even better fit to the HI data 
relative to CTEQ6L, shown in Fig. 3. At ~ 1.4 GeV^ and x ~ 10~^, a good fit to the HERA data requires 
the nonlinear evolution to be compensated by a larger gluon distribution than obtained with DGLAP alone. The 
enhancement is a factor of ~ 6 relative to CTEQ6L, as shown in Fig. 4 (left). The dependence of EHKQS 
is compared to CTEQ6L and CTEQ5L in Fig. 4 (right) for several different values of x. 

We used R = 1 fm as the free proton radius in the two-gluon density term. We did not repeat the 
calculations with different R but, depending on the transverse matter density of the free proton, some ~ 20% 
uncertainty in R can be expected. Since the nonlinear contributions decrease as R increases, a larger R would 
lead to reduced enhancement of the small x gluons below ~ 10 GeV^. Thus, minimizing the of the fit 
with respect to is a future task. 

3.2 BK approach 

authors: H. Jung, K. Kutak, K. Peters, L. Motyka 

A theoretical framework capable of describing the QCD evolution of parton densities taking gluon rescat- 
tering (corresponding to nonlinear effects) into account is the Balitsky-Kovchegov (BK) equation [24, 40, 41, 
42, 25]. The equation, based on the BFKL approach [6, 43, 7], may be used to determine the unintegrated gluon 
density. The BK equation resums the BFKL pomeron fan diagrams with the triple pomeron vertex derived in 
the high energy limit of QCD. In the doubly logarithmic limit, the BK equation reduces [25] to the colUnear 

'in fact, we produced three new sets of initial distributions, differing by the charm quark mass and parton distribution at the initial 
scale, see Ref. [37] for more details. All sets produced equally good fits to the HERA data. 
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Gribov-Levin-Ryskin (GLR) equation [22]. It is the non-collinear limit, however, which gives the dominant 
contribution to the triple pomeron vertex [44, 45]. We conclude that GLR approach misses an essential part of 
the nonlinear gluon dynamics. 

The solution to the BK equation, constrained by the low-x HERA data will be used to extrapolate the 
parton densities to the LHC kinematical domain. Extrapolation is necessary as the LHC may probe very low 
values of x, down to 10~^ for M = 10 GeV and 77 ~ 9, where unitarity corrections may be important even at 
relatively large scales of a few GeV^. Last but not least, unitarity corrections may also break kr factorization. 
We will also discuss which processes may be affected. 

This section is organized as follows. First we give a brief description of the formalism used to determine 
the gluon evolution. Within this formalism, we fit the HERA charm structure function, , data, the most 
relevant inclusive measurement directly sensitive to the gluon density. Using further assumptions about the 
sea quarks, F2 can also be described well. The resulting gluon density is then used to compute heavy quark 
production and to investigate the nonlinear effects. First we estimate bb production at CDF and DO. Then, 
cross sections for heavy quark production at various LHC experiments are estimated, tracing the impact of the 
unitarity corrections. Finally, conclusions are given. 

The standard framework to determine parton evolution is the colhnear DGLAP formahsm. It works 

rather well for inclusive quantities but, for more exclusive processes, the /cy -factorization scheme is more ap- 
propriate because both the longitudinal and transverse components of the gluon momenta are considered. In this 
framework, the process-independent quantity is the unintegrated gluon distribution, connected to the process- 
dependent hard matrix element via the /cy-factorization theorem. Linear evolution of the unintegrated gluon 
distribution may be described by one of the small x evolution equations using the fcy -factorization scheme, the 
BFKL and CCFM [46, 47, 48, 49] equations. These equations are based on resummation of large logarithmic 
pQCD corrections, a" ln™(l/a;), and are equivalent at the leading logarithmic level. 

The very small x kinematic region is also the regime where the growth of the gluon density must be 
tamed in order to preserve unitarity. Recently, a successful description of unitarity corrections to DIS was 
derived within the color dipole formulation of QCD. This is the Balitsky-Kovchegov (BK) equation which 
describes the BFKL evolution of the gluon in a large target, including a nonlinear term corresponding to gluon 
recombination at high density. 

In our analysis, we determine the unintegrated gluon distribution from the BK equation unified with the 
DGLAP equation following KMS (Kwieciriski, Martin and Stasto) [50, 51, 52, 53]. We use the abbreviation 
KKMS (Kutak, Kwieciriski, Martin and Stasto) [52, 53] for the unified nonlinear equation. The linear part 
of this equation is given by the BFKL kernel with subleading ln(l/a;) corrections, supplemented by the non- 
singular parts of the DGLAP splitting functions. Thus resummation of both the leading \nQ^ and ln(l/a;) 
terms are achieved. The subleading terms in ln(l/,7;) are approximated by the so-called consistency constraint 
and the running coupling constant. The nonlinear part is taken directly from the BK equation, ensuring that the 
unitarity constraints are preserved. One expects that this framework provides a more reliable description of the 
gluon evolution at extremely small x, where ln(l/a;) S> 1 and the unitarity corrections are important, than does 



We give a short review of the KKMS equation, starting from the impact parameter dependent BK equa- 
tion. The equation for the unintegrated gluon density, h{x,k'^,b), at impact parameter b from the center of the 
target, becomes 



the BFKL equation at LLx accuracy, extended by the negative recombination term. The (dimensionless) unin- 
tegrated gluon distribution is obtained from h{x,k'^,b) by integration over b. 



DGLAP. 



dh{x, k'^,1 
dlnl/x 



b) _ asN, 




(6) 




(7) 
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A comment about the impact parameter treatment is in order. In Eq. (7), we assume that the evolution is local in 
b. However, the complete BK equation results in some diffusion in the impact parameter plane. This diffusion 
effect may be neglected if the target is much larger than the inverse of the saturation scale. In this scheme, 
the impact parameter dependence enters through the initial condition at large xq, h{xo, k'^,b) = /{xq, k'^)S{b) 
where /(xq, fc^) is the unintegrated gluon distribution. Note that, due to nonlinearities, the b dependence of 
h{x,k'^,b) does not factorize from x and k at low x. 

The input profile function is assumed to be Gaussian, S{b) = exp(—6^/i?^)/7ri?^, where -R^ corresponds 
to the square of the average nucleon radius. Since the size of the target, R, sets the magnitude of the initial 
parton density in the impact parameter plane, h{xo, k'^,b), the unitarity corrections depend on R. At smaller R, 
gluons are more densely packed in the target and the nonlinear effects are stronger. 

References [52, 53] proposed to combine Eq. (6) with the unified BFKL-DGLAP framework developed 
in Ref. [50]. In this scheme, the (linear) BFKL part is modified by the consistency constraint [54, 55], result- 
ing in the resummation of most of the subleading corrections in ln(l/a:) which arise from imposing energy- 
momentum conservation on the BFKL kernel [56, 57, 58, 59]. In addition, we assume that the strong coupling 
constant runs with scale k'^, another source of important NLLx corrections. Finally, the non-singular part of the 
leading order DGLAP splitting function and quark singlet distribution were included in the x evolution. The 
final unproved nonUnear equation for the unintegrated gluon density is 



h{x,k'^,b) = U°\x,k'^,b)+ 



I a.(fc^)iVc ,.2 ri dz r [ Mf,fc'^&)e(^-fc'^)-Mf,fc^b) , fe(f.fc^fe) 1 



(P^A^) - ^) 4 h{l,k'^ b) + P,,(.)S (f , k'^b) 



+ 
+ 



(1 - k'^f £ f [l^ f^Mk'') In (^) h{z, k'\b) 



(8) 



The second Une of the equation corresponds to the BFKL evolution. The theta function, — fe'^), reflects 
the consistency constraint that generates the dominant part of the subleading BFICL corrections. The third line 
corresponds to the DGLAP effects generated by the part of the splitting function, Pgg{z), that is not singular 
in the limit z ^ and also by the quarks where S(a;, A;^, 6^) corresponds to the impact-parameter dependent 
singlet quark distribution. The nonlinear screening contribution following from the BK equation is given in 
the last term. The inhomogeneous contribution, defined in terms of the integrated gluon distribution, carries 
information about the transverse profile of the proton, 

m{x,k\b) = ^^S{b) j\zPgg{z)-^g , (9) 

at Aig = 1 GeV^. The initial integrated density at k^ is parameterized as 

xg{x,kl) = N{l-xY (10) 

where p = 2.5. The size of the dense gluon system inside the proton is assumed to be i? = 2.8 GeV~^, in 
accord with the diffractive slope, ~ 4 GeV~^, of the elastic J/ip photoproduction cross section at HERA. 
In this process, the impact parameter profile of the proton defines the t dependence of the elastic cross section, 
Bd ^ R^ /% by Fourier transform. In the 'hot-spot' scenario, the radius can be smaller, R=\.h GeV'^ We 
also use the hot spot value to compare with measurements and make predictions for the LHC. 

Equation (8) was solved numerically both in the linear approximation and in full. The method for solving 
Eq. (8) was developed in Refs. [50, 52]. In Fig. 5, the effects of linear and nonlinear evolution on the uninte- 
grated gluon distribution are given as a function of x for /c^ = 5 and 50 GeV^. Nonlinear evolution leads to 
sizeable suppression at the smallest x values. While the nonlinear effects are small in the HERA x range, they 
may be important at the LHC. In the following sections, we address the importance of these nonlinear effects. 

The initial distribution in Eq. (10) was obtained by fitting the HERA measurements [60, 61] using 
the Monte Carlo Cascade [62, 63] for evolution and convolution with the off-shell matrix elements. We 
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Fig. 5: The unintegrated gluon distribution obtained from Eq. (8) as a function of x for different values of k^. The solid lines correspond 
to the solution of the nonlinear equation with R = 2.8 GeV~^ while the dashed Unes correspond to the linear part. 

find per degree of freedom of 0.46 (1.17) for HI (ZEUS). The fits were repeated both with the standard 
KMS evolution without the nonlinear contribution and with extended KMS evolution including the nonlinear 
part. The predicted F2 is equivalent for both linear and nonlinear evolution, independent of R. Thus nonlinear 
evolution is only a small effect at HERA, even in the hot-spot scenario with R = 1.5 GeV~^. 




Fig. 6: The charm structure function, F^, [61] at = 4 GeV^ (a) and F2 [64] at = 6 GeV^ (b) obtained for KKMS evolution 
with different values of R. 

In Fig. 6(a) we compare the measured [61] to our prediction at = 4 GeV^. We have determined 
our initial distribution from F2 since it is the only inclusive measurement at HERA directly sensitive to the 
gluon distribution. However, we can also describe F2 [64] by making further assumptions about the sea quark 
distribution, following the KMS approach. The agreement with the data, shown in Fig. 6(b), is also quite good. 
There is only a small effect for > 5 GeV^, even in the hot-spot scenario with i? = 1.5 GeV~^. 

Next, this constrained gluon density was used to calculate gg bb production at the Tevatron as a 
cross check of the fit and the evolution formalism. We use nih = 4.75 GeV and a renormalization scale in 
as of = Ami + Pt- predicted cross section was then compared to both CDF [65, 66] and DO [67] 
measurements. The predictions agree well with the data. 

In Fig. 7(a) the cross section for B decays to J/ijj is shown as a function the J/^ pr [65, 66]. The 
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Fig. 7: Bottom production, measured by CDF, is compared to predictions using CASCADE with linear and nonlinear KKMS evolution, 
including two values of R. (a) The pr distribution of B meson decays to J/tp- (b) The azimuthal angle, A(j), distribution of bb pair 
production smeared by the experimental resolution. 



KKMS gluon density fits tlie data well in all three scenarios with deviations only appearing for px > 12 GeV. 
It is interesting to note that the approach described here gives even better agreement than the NLO coUinear 
approach [68]. 

In Fig. 7(b), the azimuthal angle distribution between the b and b quarks, Acf), is given. The A0 and 
bb pt distributions are correlated since A0 < 180° corresponds to higher pair px- Since the /cr-factorization 
formula allows the incoming gluons to have sizable transverse momenta, the calculated Acf) distribution agrees 
very well with the data for A(/) > 60° with only smearing due to the experimental resolution. The enhancement 
of the data relative to the calculations at low A(p requires further study. 

Both plots compare linear (solid histograms) and nonlinear KKMS evolution (dotted and dashed his- 
tograms) for R = 1.5 GeV~^ and 2.8 GeV^^ respectively. The nonlinear part of the evolution also has no 
impact in this kinematic region. 

4 Phenomenological applications: lieavy quark production at the LHC 
4.1 GLRMQ approach 

author: R. Vogt 

Since the HERA F2 data can be described by both hnear DGLAP and nonlinear DGLAP+GLRMQ 
evolution, as shown in Fig. 3 of Section 3.1, additional independent probes are needed. Here, we discuss 
how charm quark production in pp collisions could probe the gluon enhancement predicted in Section 3. 1 and 
described in detail in Ref. [20, 21]. Charm production is an ideal choice since the charm mass is low and its 
production is dominated by gluons. Assuming factorization, the inclusive differential charm cross section is 

dapp^ccxiQ^, Vs) =^ fi{xi,Q^) ® fj{x2,Q^) daij_^^^k}{Q'^, xi,X2) (11) 

where o"jj^cc{fc}(Q^) ^i) ^2) are the perturbatively calculable partonic cross sections for charm production at 
scales ~ Aq^^j-,, xi and X2 are the parton momentum fractions and fi{x, Q"^) are the proton parton 

densities. We assume that the renormalization and factorization scales are equal. Only the leading order gg and 
qq channels are considered here. 
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The values of the charm quark mass and scale used in the calculations are chosen to give good agreement 
with the total cross section data at NLO: m = 1.2 GeV and Q"^ = 4m^ for standard DGLAP-evolved NLO 
PDFs such as CTEQ6M [69] and MRST [70]. Nearly equivalent agreement may be obtained with m = 1.3 
GeV and = m? [71, 72]. Both choices assure that the PDFs are evaluated above the minimum scales. While 
scales proportional to m are used in the total cross section, inclusive calculations of distributions also depend 

on the transverse momentum scale, pt, so that mr = \^w?' +])\ is used instead [73]. 

To illustrate the effects of the nonlinear EHKQS distributions [37] of Section 3.1 on charm production at 
the LHC, we show 

^ da(EHKQs)/d, ^ da{^^ms)/dvT 

da {CTEQ6ih)/dy da{CTEQ6ih)/dpT 

in Fig. 8 where y is the charm quark rapidity. The results are calculated for the maximum LHC pp, pPh and 
Pb+Pb energies, = 14 (solid), 8.8 (dashed) and 5.5 (dot-dashed) TeV respectively. The results for m = 1.2 
GeV and = are on the left-hand side while those with m = 1.3 GeV and = are on the 
right-hand side. 

The change in the slope of R{y) occurs when one x drops below the minimum value of the EHKQS 
set where further nonUnearities become important, a;^?f^^^ = 10~^, and enters the unconstrained x region. 
The minimum x of CTEQ61L is lower, x^™^^^^ = 10"^ While the EHKQS gluon distribution is fixed at 
its minimum for x < xf^^'^^, the CTEQ61L distribution continues to change until a^^™^^^^- In inclusive 
kinematics with an identified charm quark and fixed xt = 2mr / VS, the unconstrained x-region contributes 
to charm production in the region 

yi = In (^1/xT - \J l/x\ - l/Xmin ^ < |y| < In [l/xT + \J 1/xf, - l/Xmin ^ • (13) 

The upper limit is close to the phase space boundary. Expanding the lower limit, yi , in powers of / Xmin ^ 1 , 
yi « ln[mT / (xrainVS)] > ln[m/(xniin\/S')]. If m = 1.2 GeV, the small x region contributes to charm 
production at |y| > yi = 2.2, 2.6 and 3.1 for s/S = 14, 8.8 and 5.5 TeV, respectively. The rather sharp 
turnover in R{y) indicates where the x < 10"^ region begins to contribute. For \y\ > yi and > A GeV^, 
as X decreases, the CTEQ61L gluon distribution increases considerably above that of the EHKQS distribution. 
Thus R{y) < 1 at large rapidities when = Am\. At midrapidity R{y) is insensitive to the EHKQS 
extrapolation region, x < x™^^^^. Since R{y) is integrated over px, it not only reflects the enhancement at 
rriT = m because oc m\ and the px distribution peaks around ~ 1 GeV. When = m\, the ratios are 
broad because the CTEQ61L gluon distribution is relatively flat as a function of x for ~ 2 — 3 GeV^. The 
enhancement decreases and broadens with decreasing energy. 

Since the rapidity distributions are rather flat, there are still important contributions to the px distributions 
from the extrapolation region, up to ~ 30% at \fS = 14 TeV for m = 1.2 GeV and = 4m^. Thus the 
sensitivity of R{pt) to the unconstrained region should be kept in mind. At the largest ^/S, the contribution 
from the x < 10~^ region is greatest and if > 4m^, xgi^'^^^^^^(x, Q^) > xg^^^'^^{x, Q^). Because the 
contribution from the region x < 10^^ decreases with ^/S, at low px R{pt) decreases with energy. In contrast, 
for = rn^, xg^^^^^{x, Q"^) > xg^'^^^'^^^{x, Q^) and the enhancement decreases with energy. 

Because the DGLAP gluon distributions are already well constrained by HERA data, they cannot absorb 
additional large effects. Therefore we conclude that, if a low-pr enhancement in the charm cross section 
relative to the DGLAP-based result is observed in future experiments, it is a signal of nonlinear effects on the 
PDF evolution. 



4.2 BK approach 

authors: H. Jung, K. Kutak, K. Peters, L. Motyka 

Since the Tevatron measurements are well described using the unintegrated parton densities constrained 
by HERA and convoluted with the off-shell matrix elements, the same approach may be used for heavy quark 
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Fig. 8: We present R{y), (a) and (c), and R{pt), (b) and (d), in pp collisions at ^/S ~ 14 (solid), 8.8 (dashed) and 
5.5 (dot-dashed) TeV. The left-hand side shows m = 1.2 GeV and = 4m^, the right-hand side m = 1.3 GeV and 




PT,bb(GeV) PT,bb(GeV) Pt.d (GeV) 



Fig. 9: (a) and (b) show hb production as a function of pair pr without cuts in 3 < |?7| < 5 (a) and in the ATLAS/CMS acceptance (b). 
The -D° meson pT distribution in the ALICE acceptance is shown in (c). 



64 



production at the LHC at e.g. y/s = 14 TeV. As discussed previously, see also Fig. 5, heavy quark production 
at this energy is already in the region where saturation effects may be relevant. In the kinematic regions, such 
as at the LHC, where nonlinear evolution may become important, the cross section will be suppressed due to 
the negative sign of the nonlinear term in Eq. (8). 

First, we compute the bb production cross section at 14 TeV without any experimental cuts. In Fig. 9(a) 
the bb differential cross section is shown as a function of pair pr in the forward region, 3 < |r7| < 5. We 
compare linear evolution (solid histogram), nonlinear evolution with R = 1.5 GeV~^ (dashed histogram) and 
R = 2.8 GeV~^ (dotted histogram). The grey band shows the uncertainty in the linear result due to the b quark 
mass. We take a central value of 4.75 GeV (the solid histogram) and vary mi, from 4.5 to 5 GeV. The bb pair 
results are shown since the pair pT is most sensitive to the gluon kx and thus to the saturation effects. In the 
hot-spot scenario, saturation effects are visible for pT^^ < 5 GeV. These saturation effects grow with rapidity, 
increasing the suppression to a factor of 3 — 4 in the higher rapidity regions. For larger R, the saturation effects 
are not very significant. 

In Fig. 9(b), the bb production cross section is computed within the ATLAS and CMS acceptance (pr > 
10 GeV and \r]\ < 2.5 for both the b and b quarks, see Ref. [74]). With these cuts, the observed suppression 
due to nonlinear effects nearly vanishes. This result suggests that kx factorization can safely be applied in the 
central r] region. Thus saturation effects due to nonlinear gluon evolution are seen only for px^- < 10 GeV and 
at high T]. This regime is accessible with upgraded ATLAS/CMS detectors or in LHCb where the b quark px 
can be measured to 2 GeV for 1.9 < < 4.9. In this kinematic regime, the hot-spot scenario predicts a factor 
of two suppression of the bb cross section. 

Similarly, we investigated cc production at ALICE. In ALICE, it will be possible to measure the 
down to pt ^ 0.5 GeV in |r/| < 0.9. The result is shown in Fig. 9(c) with rUc = 1.5 GeV. In the hot-spot 
scenario (dashed curve), a factor of two suppression occurs at ~ 1 GeV. 

5 Perspectives for experimental observation at LHC 

authors: A. Dainese and R. Vogt 

5.1 Introduction 

In Section 4.1, charm production in pp collisions at the LHC was suggested as a promising way to study 
the effects of nonlinear evolution on the parton densities. Due to gluon dominance of charm production and 
the small values of x and probed, x Ri 2 x 10~^ and 1.69 — 6 GeV^ at midrapidity and transverse 

momentum'' pT ~ 0, charm production at the LHC is sensitive to the gluon enhancement arising from nonlinear 
evolution. The resulting charm enhancement was quantified in Ref. [20, 21] by the LO ratios of the differential 
cross sections computed with the nonlinearly-evolved EHKQS PDFs [37], obtained from DGLAP+GLRMQ 
evolution, relative to tiie DGLAP-evolved CTEQ6IL PDFs. 

The enhancement of the nonlinearly-evolved gluons increases as x and decrease. Consequently, 

the charm enhancement increases with center of mass energy, ^/S. Thus the maximum enhancement at the 
LHC will be at \/S = 14 TeV and small charm quark transverse momentum. The sensitivity of the charm 
enhancement to the value of the charm quark mass, m, as well as to the choice of the factorization, 
and renormalization, Q\, scales was studied in Ref. [20, 21] assuming Q"^ = Q\ = Q\ oc m\ where 
m|, = Px + The most significant charm enhancement occurs when m and Q'^ jm^ are both small. A 
comparison of the NLO total cross sections with low energy data shows that the data prefer such small m and 
combinations [71, 72]. The largest enhancement is obtained with m = 1.3 GeV and = roj., see Fig. 8 
in Section 4.1. 

In Section 4.1, only charm enhancement was described. Neither its subsequent hadronization to D 

mesons nor its decay and detection were considered. In this section, we address these issues to determine 
whether the charm enhancement survives hadronization and decay to be measured in the ALICE detector [75]. 

*Here we use pt for the transverse momentum of the charm quark and for the transverse momentum of the D meson. 
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The calculation described in Section 4. 1 was to leading order since the EHKQS sets are evolved accord- 
ing to the LO DGLAP+GLRMQ equations using a one-loop evaluation of ag. Thus these LO distributions 
should generally not be mixed with NLO matrix elements and the two-loop a^. However, the charm quark total 
cross section is increased and the pT distribution is broadened at NLO relative to LO [76]. Thus, to determine 
whether or not the enhancement is experimentally measurable, we assume that the enhancement is the same at 
NLO and LO and employ a NLO cross section closest to the calculation of the enhancement in Ref. [20, 21]. 

As described in Ref. [76], the theoretical K factor may be defined in more than one way, depending on 
how the LO contribution to the cross section is calculated. In all cases, the ©(af ) contribution to cross section 
is calculated using NLO PDFs and the two-loop evaluation of Ug. If the LO contribution is also calculated 
using NLO PDFs and a two-loop ag, this is the "standard NLO" cross section. It is used in most NLO codes, 
both in the global analyses of the NLO PDFs and in evaluations of cross sections and rates [76]. The K 
factor formed when taking the ratio of the "standard NLO" cross section to the LO cross section with the NLO 
PDFs [76], k'^\ indicates the convergence of terms in a fixed-order calculation [77]. On the other hand, if the 
LO contribution to the total NLO cross section employs LO PDFs and the one-loop Ug, we have a cross section 
which we refer to here as "alternative NLO". The K factor calculated taking the ratio of the "alternative NLO" 
cross section to the LO cross section with LO PDFs [76], indicates the convergence of the hadronic cross 
section toward a result. If -K'q^^ > i^^2^\ convergence of the hadronic cross section is more likely [77]. This is 
indeed the case for charm production [76]. We also note that kI^^ is a much weaker function of energy than 
k'^K Since, in the absence of nonlinear NLO PDFs, the "alternative NLO" cross section is more consistent 
with the calculated enhancement, we use this cross section to calculate the NLO D meson rates and px spectra. 
We note also that, in both cases, the pt distributions have the same slope even though K2'\ for the alternative 
NLO cross section, is somewhat smaller Thus, using anon-standard NLO calculation will not change the slope 
of the Pt distributions, distorting the result. 

The LO and NLO calculations used to obtain the full NLO result in both cases can be defined by modi- 
fication of Eq. (11) in Section 4.1. We define the full LO charm production cross section as 

d<yto= E /?°(^i,QF)®/i^°(^2,Ql.)®da^°,,K^^g|),xi,X2) (14) 

where the superscript "LO" on ddij^cc indicates the use of the LO matrix elements while the superscript "IL" 
indicates that the one-loop expression of is used. The LO cross section typically used in NLO codes employs 
the NLO PDFs and the two-loop (2L) Ug so that 

dc^to= E ff''''{xi.Ql)<^fr''{x2Ml)®da\^_^^{af{Ql),xi,X2). (15) 

i,j=<l,l,9 

In either case, the NLO contribution, ©(ctf ) for heavy quark production, is 

rf^O(a3)= E /f'^°(^i>QF)®/r°(^2,Q|)C3 5^ daf^Sc-^{af{Ql),Ql,x,,X2) (16) 

i,j=(l,g,g k=0,q,q,g 

where the superscript "NLO" on daij^cck indicates the use of the NLO matrix elements. The additional sum 
over k in Eq. (16) includes the virtual (k = 0) and real (k = q,q or g depending on i and j) NLO corrections. 
In the calculations of da^ and dcr^^^s), we use the value of Aq^j-, given for the NLO PDFs and work in the 
MS scheme. The standard NLO cross section is then 

while our "alternative NLO" cross section is defined as 

dc^^io = dcrjx) + dao^a^) ■ (18) 

Since the enhancement in Ref. [20, 21] was defined using da^ only, the best we can do is to use the alternative 
NLO cross section in our analysis, as described below. 
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We now discuss how the enhancement is taken into account in the context of the NLO computation. 
We calculate the LO inclusive charm px distribution, cPa/dprdy, with the detected charm (anticharm) quark 
in the rapidity interval Ay with \y\ < 1, motivated by the pseudorapidity acceptance of the ALICE tracking 
barrel, {rjl < 0.9. The rapidity, y2, of the undetected anticharm (charm) quark is integrated over. The charm 
enhancement factor R{pt, Ay) is then 

, f , d^crfEHKQS) 
dy / dy2- 



m A ^ ./At/ J dpTdydy2 

JAy J dpTdydy2 

Next, we assume that the enhancement calculated at LO is the same when calculated at NLO. This is the 
only reasonable assumption we can make to test whether the enhancement can be detected with ALICE which 
will measure the physical p^ distribution. The alternative NLO cross section is therefore the closest in spirit to 
the LO computation in Ref. [20, 21]. Thus, the enhanced NLO charm pT distribution is 

R{PT, Ay) da^io {Ay)/dpT . (20) 

Our results are obtained with the same parameters used in Section 4. 1, m = 1.2 GeV and = Arrv^ as well as 
m = 1.3 GeV and = m^. These two choices are the baseline results against which other parameter choices 
will be compared to see if the enhancement is detectable. 



5.2 From charm to D enhancement 

To make a more realistic D meson distribution, we have modified the charm pT distribution by the heavy quark 
string fragmentation in PYTHIA [78]. Charm events in pp collisions at = 14 TeV are generated using 
PYTHIA (default settings) with the requirement that one of the quarks is in the interval \y\ < 1. The charm 
quarks are hadronized using the default string model. Since c and c quarks fragment to D and D mesons^, 
respectively, in each event related (c, D) and (c, D) pairs can easily be identified^. These pairs are reweighted 
to match an arbitrary NLO charm quark pT distribution, dN^-^Q/dpr- If d^PYTKiA/ ^Pt is the charm pT 
distribution given by PYTHIA, each (c, D) pair is assigned the weight 

WW = /J!^'-°'t (21) 

where pr is the transverse momentum of the charm quark of the pair. Therefore, the reweighted final-state 
D distribution corresponds to the one that would be obtained by applying string fragmentation to the NLO 
c-quark distribution. The resulting D distribution is significantly harder than that obtained using the Peterson 
fragmentation function [79]. The enhancement survives after fragmentation although the D enhancement is 
somewhat lower than that of the charm because for a given p^, the D spectrum receives contributions from 
charm quarks with pT ~ , where the charm enhancement is smaller. 



5.3 Sensitivity to the enhancement 

Figure 10 shows the double-differential cross section, d^ao/dp^^dy, in \y\ < 1 as a function of the trans- 
verse momentum. The points represent the expected "data" measured by ALICE, obtained from the alternative 
NLO cross section scaled by the enhancement factor R{pt, Ay) defined in Eq. (19), and modified by string 
fragmentation. The solid and dashed curves are obtained by applying string fragmentation to the alternative 
NLO and standard NLO cc cross sections, respectively. Thus, the "data" points include the enhancement while 
the curves do not. The horizontal error bars indicate the bin width, the vertical error bars represent the statistical 
error and the shaded band gives the pr-dependent systematic error. The 5% pT-independent systematic error on 
the normalization is not shown. (See Ref. [80] for a discussion of the error analysis. The standard NLO cross 

^Herei) s D+,D°. 

^Events containing charm baryons were rejected. 
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Fig. 10: Comparison of the simulated ALICE data generated from R{pT , Ay)da^l^Q with the ahernative (soHd) and 
standard (dashed) NLO calculations. The effect of string fragmentation is included in the "data" points as well as in the 
curves. The left-hand side shows the result for m = 1.2 GeV and = 4m|, while the right-hand side is the result for 
m = 1.3 GeV and = m^. The error bars on the data represent the statistical error and the shaded band represents the 
pr-dependent systematic error. The 5% normaUzation error is not shown. 

section, Eq. (17), and the C(af) contribution to the alternative NLO cross section, Eq. (16), were calculated 
using the HVQMNR code [81] with CTEQ6M and A^^^ = 0.326 GeV. The LO contribution to the alternative 
NLO cross section, Eq. (14), was calculated using the CTEQ61L PDFs. Fragmentation was included as de- 
scribed in Section 5.2. The enhancement, the difference between the data and the solid curves for pl^ ^ 3 GeV, 
is more pronounced for the larger mass and lower scale, on the right-hand side of Fig. 10. 

There is a significant difference between the alternative and standard NLO distributions. Part of the 
difference is due to the one- and two-loop evaluations of since af^ < a]^. However, the most important 
contribution is the large differences between the LO and NLO gluon distributions, especially at low scales [80]. 

In order to address the question of the experimental sensitivity to the effect of nonlinear gluon evolution 
on low-pT charm production, we consider, as a function of p^, the ratio of the simulated data, including the 
enhancement, to alternative NLO calculations using a range of m and with PYTHIA string fragmentation. 
We denote this ratio as "Data/Theory" and try to reproduce this ratio with NLO calculations employing recent 
linearly-evolved PDFs and tuning m and Q^. 

Since the enhancement has disappeared for p!^ 5 GeV, we refer to this unenhanced region as high 
p^. The pf^ region below 5 GeV, where the enhancement is important, is referred to as low p^. If no set of 
parameters can describe both the high- and low-p^ components of the distribution equally well, and, if the set 
that best reproduces the high-py part underestimates the low-p^ part, this would be a strong indication of the 
presence of nonlinear effects. 

The Data/Theory plots are shown in Fig. 11. The points with the statistical (vertical bars) and pT- 
dependent systematic (shaded region) error correspond to the data of Fig. 10, including the enhancement, 
divided by themselves, depicting the sensitivity to the theory calculations. The black squares on the right-hand 
sides of the lines Data/Theory = 1 represent the 5% -independent error on the ratio coming from the cross 
section normalization. This error is negligible relative to present estimates of other systematic uncertainties 
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Fig. 11: Ratio of the generated ALICE data relative to calculations of the alternative NLO cross sections with several sets 
of parameters and PYTHIA string fragmentation. The left-hand side shows the result for m = 1.2 GeV and = Arn^ 
while the right-hand side is the result for m — 1.3 GeV and = m^. 



(~ 13%). 

On the left-hand side, the thick solid curve with m = 1.2 GeV and = 4m|. best agrees with the 
high-p^ ratio by construction since i? « 1 at large p^. It also shows the effect of the enhancement well beyond 
the error band for < 2 GeV. Better agreement with the data over the entire range can be achieved only by 
choosing a charm quark mass lower than 1.2 GeV, below the nominal range of charm masses, as illustrated by 
the dashed curve for m = 1.1 GeV. Higher masses with = produce much larger Data/Theory ratios 
than the input distribution. The ratio with m = 1.3 GeV and = m\ (dot-dot-dashed curve) gives a much 
larger ratio at low and drops below « 1 for > 8 GeV. 

We also present the ratio using the MRST parton densities (MRST2001 LO [36] in Eq. (14) and MRST2002 
NLO [82] in Eq. (16)) with m = 1.2 GeV and = Am^. We find that this result also agrees reasonably well 
with the CTEQ6 results for the same m and Q^. Thus, the enhancement seems to be rather independent of the 
PDF. The CTEQ61L and the MRST2001 LO distributions are similar at low x, suggesting that non-linearly 
evolved PDFs based on MRST2001 LO would produce an enhancement like that of Ref. [20, 21]. 

On the right-hand side of Fig. 1 1 the thick solid curve, employing the same parameters as the data, gives 
the best agreement at high p^. We note that even though the results with = 4m|. and m < 1.3 GeV lie 
closer to the data at low and within the combined statistical and systematic error at higher pl^, the curves 
with these parameters have the wrong slopes for ^ 8 GeV. The statistical sensitivity is expected to be good 
enough to distinguish the difference in curvature. The results with the MRST PDFs do not alter the conclusions. 



5.4 Conclusions 

We have studied whether the EHKQS gluon distributions [37] could generate an observable D meson enhance- 
ment in pp collisions at the LHC. Using the EHKQS LO PDFs and LO matrix elements for charm quark 
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production and PYTHIA string fragmentation for D meson hadronization, the enhancement indeed survives to 
the D mesons. 

The D meson enhancement, however, drops rapidly with transverse momentum. Therefore, D measure- 
ment capability at small is necessary to verify the effect experimentally. The ALICE detector can directly 
reconstruct K~tt^. We have demonstrated that, in the most optimistic case, the enhancement can be 

detected above the experimental statistical and systematic errors. When the charm mass is somewhat smaller, 
m = 1.2 GeV, but the scale is larger, = Arn^, it is more difficult to detect the enhancement over the 
experimental uncertainties. 
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Abstract 

The fragmentation of heavy quarks into hadrons is a key non-perturbative ingredient 
for the heavy quark production calculations. The formalism is reviewed, and the 
extraction of non-perturbative parameters from e~^e~ and from ep data is discussed. 

coordinator: M. Corradi 



1 Introduction 

author: M. Cacciari 

When we try to describe in QCD the production of a hadron we are always faced with the necessity 
to take into account the non-perturbative hadronization phase, i.e. the processes which transform perturbative 
objects (quarks and gluons) into real particles. In the case of light hadrons the QCD factorization theorem [1,2, 
3, 4, 5, 6] allows to factorize these non-perturbative effects into universal (but factorization-scheme dependent) 
fragmentation functions (FF): 

In this equation, valid up to higher twist corrections of order A/pr (A being a hadronic scale of the order 
of a few hundred MeV and pt for instance a transverse momentum), the partonic cross sections dai/dpx for 
production of the parton i are calculated in perturbative QCD, while the fragmentation functions Di^h{x; fi) 
are usually extracted from fits to experimental data. Thanks to their universality they can be used for predictions 
in different processes. The artificial factorization scale is a reminder of the non-physical character of both 
the partonic cross sections and the fragmentation functions: it is usually taken of the order of the hard scale pT 
of the process, and the fragmentation functions are evolved from a low scale up to /x by means of the DGLAP 
evolution equations. 

This general picture becomes somewhat different when we want to calculate the production of heavy- 
flavoured mesons. In fact, thanks to the large mass of the charm and the bottom quark, acting as a cutoff 
for the collinear singularities which appear in higher orders in perturbative calculations, one can calculate the 
perturbative prediction for heavy quark production. Still, of course, the quark hadron transition must be 
described. Mimicking the factorization theorem given above, it has become customary to complement the 
perturbative calculation for heavy quark production with a non-perturbative fragmentation function accounting 
for its hadronization into a meson: 

dan , X (' dx da^''* ^ ^ ^^^^ ^ ^ 



dp' 



It is worth noting that at this stage this formula is not given by a rigorous theorem, but rather by some sensible 
assumptions. Moreover, it will in general fail (or at least be subject to large uncertainties) in the region where 
the mass m of the heavy quark is not much larger than its transverse momentum px, since the choice of the 
scahng variable is not unique any more, and 0{m/pT) corrections cannot be neglected. 
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Basic arguments in QCD allow to identify the main characteristics of the non-perturbative fragmentation 
function Dq'_^^{x). In 1977 J.D. Bjorken [7] and M. Suzuki [8] independently argued that the average fraction 
of momentum lost by the heavy quark when hadronizing into a heavy-flavoured hadron is given by 

~ 1 - A . (3) 

m 

Since (by definition) the mass of a heavy quark is much larger than a hadronic scale A, this amounts to saying 
that the non-perturbative FF for a heavy quark is very hard, i.e. the quark loses very little momentum when 
hadronizing. This can also be seen with a very simplistic argument: a fast- traveling massive quark will lose 
very little speed (and hence momentum) when picking up from the vacuum a light quark of mass A to form a 
heavy meson^ . 

This basic behaviour is to be found as a common trait in all the non-perturbative heavy quark FFs derived 
from various phenomenological models. Among the most commonly used ones we can cite the Kartvelishvih- 
Likhoded-Petrov [12], Bowler [13], Peterson-Schlatter-Schmitt-Zerwas [14] and ColUns-Spiller [15] fragmen- 
tation functions. These models all provide some functional form for the D^^jj{x) function, and one or more 
free parameters which control its hardness. Such parameters are usually not predicted by the models (except 
perhaps on an order-of-magnitude basis), and must be fitted to the experimental data. 

During the '80s many such fits were performed, and these and similar functions were also included in 
many Monte Carlo event generators. Eventually, some 'best' set of parameter values (for instance for the PSSZ 
form) was determined [16] and subsequently widely used. 

These first applications, given the limited accuracy of the available data, tended to overlook two aspects 
which have become more important in recent years, when the accuracy of the data has vastly improved: 

• A non-perturbative FF is designed to describe the heavy quark — > hadron transition, dealing with events 
mainly populated by soft gluons of energies of a few hundred MeV. However, if a heavy quark is produced 
in a high energy event it will initially be far off shell: perturbative hard gluons will be emitted to bring it 
on-shell, reducing the heavy quark momentum and yielding in the process large collinear logarithms (for 
instance of the form a'g \og^\pT/m) in a transverse momentum differential cross section). Of course, the 
amount of gluon radiation is related to the distance between the heavy quark mass scale and the hard scale 
of the interaction, and is therefore process-dependent. One can (and it was indeed done) either fit different 
free parameters at different centre-of-mass energies (or transverse momenta), or try to evolve directly the 
non-perturbative FF by means of the DGLAP equations, hence including into it the perturbative collinear 
logarithms. However, this is not what non-perturbative fragmentation functions are meant for, and doing 
so spoils the validity of the relation in eq. (3). 

• Since only the final heavy hadron is observed, both the non-perturbative FF and the perturbative cross 
section for producing the heavy quark must be regarded as non-physical objects. The details of the fitted 
non-perturbative FF (e.g. the precise value(s) of its free parameter(s)) will depend on those of the per- 
turbative cross sections: different perturbative calculations (leading order, next-to-leading order, Monte 
Carlo, ...) and different perturbative parameters (heavy quark masses, strong coupling, ...) will lead to 
different non-perturbative FFs. These in turn will have to be used only with a perturbative description 
similar to the one they have been determined with. Hence the hmited accuracy (and hence usefulness) of 
a 'standard' determination of the parameters [17]. 

The first point was addressed by Mele and Nason in a paper [18] which deeply changed the field of 
heavy quark fragmentation, and essentially propelled it into the modem era. Mele and Nason observed that, 
in the limit where one neglects heavy quark mass terms suppressed by a large energy scale, a heavy quark 
cross section can be factored into a massless, MS-subtracted cross section for producing a light parton, and a 
process-independent-^, perturbative heavy quark fragmentation function describing the transition of the massless 

'More modem and more rigorous derivations of this result can be found in [9, 10, 1 1]. 

"Mele and Nason extracted this function from the e^e^ cross section, convincingly conjecturing its process independence, which 
was successively established on more general grounds in [19] 
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Fig. 1: Power-law fits to the heavy quark pr distributions at LHC (left) and HERA (right) obtained with the NLO programs MNR and 
FMNR. The resulting exponents are iV = 4.5/3.8 for charm^eauty at LHC and iV = 5.5/5.0 for c/b at HERA. 



parton into the heavy quark: 

7 pGVt TG.S 

^^(fr,m) = E/-^F;'') A,0(x;,.m) + 0(^) . (4) 
dpT J X dpT V X / \pt J 

The key feature of this equation is that it is entirely perturbative: every term can be calculated in perturbative 
QCD. The perturbative fragmentation functions Di^Q{x; fi, m) (not to be confused with the non-perturbative 
one D^^^{x)) can be evolved via DGLAP equations from an initial scale of the order of the heavy quark mass 
up to the large scale of the order of pT- This resums to all orders in the strong coupling the collinear logarithms 
generated by the gluon emissions which bring the heavy quark on its mass shell, leading to a more accurate 
theoretical prediction for dag/dpT- 

Once a reliable perturbative cross section for the production of a heavy quark is established, one is simply 
left with the need to account for its hadronization. For this purpose one of the functional forms listed above 
can be used for the non-perturbative FF, and implemented as in eq. (2), but using the improved, resummed 
cross section given by eq. (4). Since most of the the scaling-violation logarithms are accounted for by the 
evolution of the perturbative FF, the non-perturbative one can now be scale-independent and only contain the 
physics related to the hadronization of the heavy quark. It will always, however, depend on the details of the 
perturbative picture used. 

2 Extraction of heavy quark fragmentation parameters from e+e^ and their impact on HERA and 
LHC 

author: M. Corradi 



2.1 Importance of (x)'^p 

According to the factorization of the fragmentation functions (FF), the differential cross section da/dpT for the 
production of a heavy hadron H can be written as the convolution of the perturbative heavy quark differential 
cross section da^^^^ /dpx and the non-perturbative fragmentation function D"p(x): 

da , ^ f dx , dcrP®''* 



f'j!iD"i./^ (?!.). (5) 

J X dpT \ X / 



dpi 

This convolution neglects mass terms 0{mQ/pT) and non-perturbative terms (!?(Aqcd /?tiq). 

The heavy quark px distribution behaves at large px like a power law da^^^^/dpT = Cp^^ . Figure 1 
shows power-law fits to the p'r distributions of heavy quarks at LHC and in photoproduction at HERA as 
obtained from the NLO programs MNR [20] and FMNR [21]. For pr > 10 GeV N was found to range from 
3.8 (b at LHC) to 5.5 (c at HERA). Combining this power-law behavior with eq. (5), the hadron p^ distribution 
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Function 


D{x) 




p 3.rflrn. c tcrs 


X range 


Kartvelishvili 


(1 - x)x" 




a = 2/5 — 3 


[0,1] 


Peterson 


1 /, 1 

1 (l _ i _ 

X \ X 




) 


[0, 1] 






Gauss 


exp(-(2^ 


f) 


a — 6 /2 


[— oo, oo] 


Box 


const. 






[1 - 25, 1] 


Triangular: 


X — XQ 




XO = 1 — 35 


[1 - 35, 1] 



Table 1: Test functions used in Fig. 2. The functions assume a value different from zero in the range given by the third column. 





<x>=0.666 




Fig. 2: Effect of the convolution of the heavy quark transverse momentum distribution with different test functions for different values 
of the non-perturbative FF (x) = 0.9 (left), 0.8 (center), 0.666 (right). For each (a;), the upper plot shows the test functions, the middle 
plot shows the perturbative pt distribution obtained with the MNR program for beauty at LHC and the hadron pr distributions after 
the convolution with the test functions. The lower plot shows the ratio of the different hadron pr distributions to the result obtained 
with the Peterson one. 



is given by 



where D'^ 



da 

dpT 



{P7 



dxx^-^ D''^{x)Cp[ 



-N 



dpT 



(6) 



Jdx 



X 



Z)"P(x) is the N Mellin moment of the non-perturbative FF. 



The hadron distribution is therefore governed by the 4^^ - 5^^ Mellin moments of Z)'^p(x). It is interesting 
to translate the Mellin moments into more intuitive central moments 



^J'n 



(x))"D'"P(x) for n>2 



(7) 



where {x) = f dx xD^p(x) is the mean value. The first Mellin moments, written in terms of (x) and are: 

£>! = 1 , 1)2 = (x) , 1)3 = {x)"^ + ^2 , £'4 = {xf + 3fl2{x) + /U3. 

In heavy quark fragmentation, the mean value of D'^p(x) can be written as (x) = 1 — 5 where 5 = 
C(^qcd/"iQ) is small [11]. For any positive function with (x) = 1 — S, defined in the interval [0,1], the 
central moments are limited by 5, \fin\ < 6. In practice, reasonable heavy quark fragmentation functions are 
concentrated in a small region around 1 — 6 and therefore the higher central moments are small. To be specific, 
if the function is different from zero in a region of size ±KS (with K = 0{l)) around 1 — 6 then |/U„| < {K6)"'. 
This means that the Mellin moments of reasonable FFs are given, to a good approximation, by the mean value 
to the — 1 power: 

(8) 

{x)''-' + §^f^2{xr-' + 0{6^).For 
a reasonable FF and a perturbative distribution falling with the power —N, Eq. 6 and 8 give 



DN = {xf-' + 0{6^) 
The expansion to 6"^ involves the second central moment fi2- Dj^ = 



da 
dpT 



(pt) 



dp^ 



(9) 
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Therefore the effect of the non-perturbative FF is to introduce a shift in the normalisation that depends on 
the average x, while the details of the shape of D{x) have negligible effect. To check that this reasoning 
works with realistic fragmentation functions and realistic perturbative pT distributions, various functions with 
the same (x) but different shapes have been tested in convolution with the perturbative px spectrum for b 
production at LHC obtained with the NLO program MNR. The test functions considered are the Peterson [14] 
and Kartvelishvili [12] fragmentation functions, a Gaussian distribution with a = 1 — fi, a. flat and a triangular 
distribution. Table 1 gives more detail about these functions. Three average values were chosen: (x) = 
0.9, 0.8, 0.666. Figure 2 shows the result of this test. For each average value, the convolutions are very similar, 
even if the test functions are very different. For (x) = 0.9, 0.8, which are typical beauty or charm values, the 
hadron spectra agree within few %. For the extreme value of (x) = 0.666, the results for the Peterson and 
Kartvelishvili functions give very similar hadron spectra while the less realistic Gaussian and Box shapes differ 
at most 10% from Peterson at large px and the extreme Triangular function shows deviations up to ~ 20%. 

In conclusion the relevant fragmentation parameter for the inclusive hadron spectra at pp and ep colliders 
is the mean value of the non-perturbative FF. The next part will discuss, on the basis of e+e^ data, what 
values of {x)^^ are relevant for different calculations. 



2.2 Extraction of (x)°p from e+e data 

In e+e^ interactions it is convenient to express the factorization ansatz, given for the heavy-hadron pT in eq. (5), 
in terms of the heavy-hadron momentum normalized to the maximum available momentum: Xp = p^/Pmax' 

where p^^^ = ^ (i^cms)^ - rnj^- 

which corresponds to the following relation for the mean values: (xp) = {x)'^'^ (x)^^^^ where {xp) is the mean 
hadron Xp, (x)"p is the mean value of the non-perturbative FF and (x)^^"^^ = f dx x ^^j— is the mean 
value of the perturbative distribution. Then, taking (xp) from experimental data and (x)p'^''* from a particular 
perturbative calculation, it is possible to extract the value of (x)'^p valid for that calculation as 

(x)°P = (xp)/(x)P-*. (10) 



Two perturbative calculations will be considered to extract (x)p°'^*: a fixed-order (FO) next-to-leading 
order (i.e. 0{as)) calculation and a calculation that includes also the resummation of next-to-leading loga- 
rithms (NLL) and Sudakov resummation, both obtained with the HVQF program [19]. From the point of view 
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Charm (D*) 




Measured 


Value 


ISR corr. 


X ^ ^ Xp 


{Xp)""" 


measurement 


(GeV) 


variable 




(%) 


(%) 




OPAL [22] 


92 


X'^E/ 


516+1!1^1^2 ± 010 


+0.4 


—0.4 


0.516 ± 0.012 


ALEPH [23] 


92 


V'EI 


n 4878 + 0046 + OOfil 


4-0 4 


—0.4 


488 + 008 


DELPHI [24] 


92 


{xe) 


0.487 ±0.015 ±0.005 


+0.4 


-0.4 


0.487 ±0.016 


TASSO [25] 


36.2 


(xp.) 


0.58 ± 0.02 


+6.7 


-1.8 


0.61 ±0.02 


DELCO [26] 


29 




31+S i2 


+6.3 




0.55 ± 0.03 


ARGUS [27] 


10.5 


v^V 1 


0.64 ± 0.03 


+4.2 




0.67 ± 0.03 


CLEO [28] 


10.5 


{Xp) 


0.611 ±0.007 ±0.004 


+4.2 




0.637 ± 0.008 


BELLE [29] 


10.58 


(Xp) 


0.61217 ± 0.00036 ± 0.00143 


+4.2 




0.6379 ± 0.0016 




Beauty (S""*) 


-^cms 


Measured 


Value 


ISR corr. 


Xe —* Xp 


(Xp)'=°" 


measurement 


(GeV) 


variable 




(%) 


(%) 




OPAL [30] 


92 


(xe) 


0.7193 ± 0.0016 ± 0.0038 


+0.3 


-0.9 


0.715 ±0.004 


SLD [31] 


92 


{xe) 


0.709 ± 0.003 ± 0.005 


+0.3 


-0.9 


0.705 ± 0.006 


ALEPH [32] 


92 


(xe) 


0.716 ± 0.006 ± 0.006 


+0.3 


-0.9 


0.712 ± 0.008 


DELPHI [33] 


92 


(xe) 


0.7153 ±0.0007 ±0.0050 


+0.3 


-0.9 


0.711 ±0.005 


JADE [34] 


36.2 


(xe) 


0.76 ± 0.03 ± 0.04 


+5.4 


-3.5 


0.77 ± 0.06 


DELCO [35] 


29 


(xe) 


0.72 ± 0.05 


+4.8 


-4.7 


0.72 ± 0.05 


PEP4-TPC [36] 


29 


(xe) 


0.77 ± 0.04 ± 0.03 


+4.8 


-4.7 


0.77 ± 0.07 



Table 2: Experimental results on the average fragmentation function in e^e~ collisions for D* mesons and weakly decaying beauty 
hadrons. The table reports, for each experiment, the published variable and the corrections applied to obtain {xp)'^°^^ . All the mea- 
surement have been corrected for initial state radiation (ISR). Measurements reported in terms of (xe) have been corrected to (xp). 
In the case of ARGUS the average has been calculated from the full distribution. In the case of TASSO the error on the average was 
re-evaluated using the full distribution since the published error seems incompatible with the data. DELCO reports a fit with a Peterson 
distribution that has been translated into {xp)'^"" . Systematical and statistical uncertainties, where reported separately, have been added 
in quadrature. The ALEPH beauty measurement refers to and B'^ mesons only (i.e. excluding Bg and A;,), a MC study shows that 
this correspond to underestimating {xp)^°''^ by ~ 0.1% only, which is negligible. 

of fragmentation, the FO calculation only considers the emission of a gluon from one of the two heavy quarks 
generated in the e+e" collision while the NLL calculation includes the evolution of the FF from the hard in- 
teraction scale down to the scale given by the heavy quark mass. The parameters used for the FO and NLL 
models are rUc = 1.5 GeV, = 4.75 GeV, Aqcd = 0.226 GeV and the renormalisation and factorization 
scales IJ.R = i^F = Ecms- The starting scale for FF evolution in the NLL model was chosen to be mq. The 
theoretical uncertainty was obtained by varying independently the normalisation and factorization scales by a 
factor 2 and 1/2 and taking the largest positive and negative variations as the uncertainty. 

The experimental data are also compared to the PYTHlA 6.2 Monte Carlo program [37] which con- 
tains an effective resummation of leading-logarithms based on a parton-shower algorithm and which is inter- 
faced to the Lund fragmentation model. In this case the MC model gives directly {xp), while (x)"^^^^ has 
been obtained taking the heavy quark at the end of the parton shower phase. The quark masses have been 
set to rric = 1.5 GeV and nih = 4.75 GeV, and all the parameters were set to the default values except 
for specific fragmentation parameters explained below. Three sets of fragmentation parameters were cho- 
sen for charm: the default fragmentation (Lund-Bowler), a longitudinal string fragmentation of the Peterson 
form with e = 0.06 (mstj(ii)=3, parj(54)=-o.06) and the Lund-Bowler fragmentation with parameters re- 
tuned by the CLEO collaboration [28] (PARJ(41)=0.178, PARJ (42) =0 . 393, PARJ (13) =0 . 627). The tWO SCtS cho- 

sen for beauty are the default Lund-Bowler fragmentation and the Peterson fragmentation with e = 0.002 
(mstj(ii)=3, parj(55)=-o.oo2). Figurc 3 shows (x)^^^^ from the perturbative calculations as a function of the 
centre of mass energy -Ecms for charm and bottom. 

2.3 Charm 

Charm fragmentation data are available from various e+e^ experiments. The most precise are those at the 
pole at LEP (ALEPH [23], OPAL [22], DELPHI [24]) and near the T(4s) (ARGUS [27], CLEO [28], 
BELLE [29]). Less precise data are available in the intermediate continuum region from DELCO [26] at 
PEP and TASSO [25] at PETRA. Measurements in which the beauty component was not subtracted have 
been discarded [38, 39, 40]. The experimental data are reported in Table 2. Only measurements relative to 
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Fig. 4: Average fragmentation function as a function of the center of mass energy for cliarm (left) and beauty (right). The plots show 
the experimental results and the curves from NLL and FO theory with a non-perturbative fragmentation obtained using the data at the 
energy. The curves from Pythia 6.2 with different fragmentation choices are also shown. The experimental points at the T(4s) 
and resonances are shown slightly displaced in the horizontal axis for better legibility. 



the 15*^(2010) meson are considered, to avoid the comphcations due to cascade decays that are present for 
ground state mesons. Charm quarks originating from gluon sphtting rather than from the virtual boson from 
e+e~ annihilation may be relevant at LEP energies. This contribution is anyway already subtracted in the 
published data considered here, and it is consistently not considered in the perturbative calculations. Most of 
the experiments published the mean value of the x distribution. The only exception is ARGUS, for which the 
mean value was computed from the published distribution. Some of the experiments give the results directly 
in terms of Xp, others in terms of the energy fraction xe = '^Eh / Ecms- The latter has been corrected to Xp 
using the Pythia MC. The difference between {xp) and (xe) can be as large as 12% at ii^cms = 10.5 GeV 
and reduces to less than 1% at ii^cms = 92 GeV. Since the low-mass measurements are already given in terms 
of Xp, the applied corrections from xe to Xp was always small. QED corrections are also needed to compare 
the experimental data to the QCD predictions. The initial state radiation (ISR) from the electrons has the effect 
of reducing the energy available for the e+e^ annihilation and therefore to reduce the observed value of (xp). 
A correction, obtained by comparing the Pythia MC with and without ISR, was applied to the data to obtain 
{xpy°". The correction is ~ 4% at E'cms = 10.5 GeV, is largest in the intermediate region and is negligible at 
^cms = 92 GeV 

Only LEP data at i?cms = 92 GeV were used to extract (x)'^p since the factorization of the non- 
perturbative FF could be spoiled by large 0{mQ/Ecms) terms at lower energies. Table 3 reports the LEP 
average {xpY°^^ , the perturbative results at 92 GeV and the resulting (x)'^p for NLL and FO calculations as 
well as (x) and (x)p^'^* from Pythia. Figure 4 (left) shows {xp) obtained by multiplying the perturbative 
calculations with the corresponding (x)'^p, compared to the experimental data and to the Pythia MC with 
different fragmentation parameters. 

With the non-perturbative (x)'^p = 0.849 ± 0.018 obtained at LEP energies, the NLL calculation can 
reproduce all the data within a quite small theoretical uncertainty. The FO calculation is instead too flat to 
reproduce the data even considering its large theoretical uncertainty band. The non-perturbative fragmentation 
(x)°P obtained at LEP energy for the FO calculation is quite small (0.65 ib 0.04) since it compensates the effect 
of the FF evolution that is missing in the perturbative part. Therefore FO calculations with (x) extracted at 
LEP energy undershoot drastically the data at the T(4s). 

The Pythia MC with the Lund-Bowler fragmentation reproduces the data reasonably well. The result 
with default parameters is slightly above the data while the result with the parameters tuned by the CLEO col- 
laboration is slightly below. Both are compatible within the experimental uncertainty with all the experimental 
values with the exception of the very precise measurement from Belle from which they differ anyway by less 
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92 GeV 








Data 










NLL 






W.OOiJ _1_ \J .\J\J t 




FO 






n 7f| -L n (iQ 




Pythia Def . 




0.500 


0.640 




Pytiiia CLEO 




0.181 


0.610 




Py iiua Pet. ( : 


= O.Oli 


0. 190 


0.(> 10 






Beauty (S"-") i 


g> 92 GeV 








Data 




0.7114 ± 0.0026 






NLL 






0.768 ±0.010 


0.927 ±0.013 


FO 






0.83 ±0.02 


0.85 ± 0.02 


Pythia Def. 




0.686 


0.773 




Pythia Pet. e : 


= 0.002 


0.710 


0.773 





Table 3: Average fragmentation functions at the Z resonance for charm (top) and beauty. The table shows the average of the ex- 
perimental data, the results from the NLL and FO calculations and from the Pythia MC with different fragmentation parameters. 
For the NLL and FO calculations (xp)"^ is obtained by dividing the average from the experimental data by the perturbative result 

{Xp) ^ = {Xp) / {Xp)^ 

than 2%. PYTHIA with the Peterson fragmentation with e = 0.06 reproduces well the LEP data but is too low 
at lower energies. 

2.4 Beauty 

In the case of beauty we consider fragmentation measurements for the mix of weakly decaying hadrons B"^^. 
Precise measurements are available only at the peak (SLD [31], ALEPH [32], OPAL [30], DELPHI [33]). 

Lower energy measurements from PEP (PEP4-TPC [36], DELCO [35]) and PETRA (JADE [34]) have larger 
uncertainties. As for charm, corrections have been applied for ISR and to convert {xe) to {xp). The data 
are shown in Table 2 and the results in Table 3 and Figure 4 (right). Since precise data are available only at 
a single energy, it is impossible to test the energy beaviour of the theoretical predictions. As in the charm 
case, the energy dependence of Pythia and NLL theory are similar, while the FO prediction is much more 
flat, suggesting that also for beauty the non-perturbative fragmentation obtained for FO at the could not be 
applied at lower energy. Pythia with Peterson fragmentation with e = 0.002 reproduces the data, while the 
default Lund-Bowler fragmentation is too soft. 

2.5 Effect on predictions for heavy quark production at HERA and LHC 

Going back to the heavy-hadron production in ep and pp collisions, eq. 9 shows that the uncertainty on the 
differential heavy-hadron cross section da/dpT is related to the uncertainty on the average non-perturbative 
fragmentation by 

A{da/dpT) = NA{{x)''P), 

where — iV is the exponent of the differential cross section. 

The state of the art calculations for photo- and hadro-production (FONLL [41, 42]) include NLO matrix 

elements and the resummations of next-to-leading logarithms. The appropriate non-perturbative fragmentation 
for FONLL is therefore obtained with the NLL theory which has the same kind of perturbative accuracy [43]. 
Since the NLL calculation gives a good description of e'^e~ data, it seems appropriate to use the value and the 
uncertainty of (x)'^p as obtained from e+e" data at the Z^ peak. The relative error for the D* fragmentation 
is A(a:)"P/(a:)"P = 2% which translates into an uncertainty of 9% on charm production at large px at LHC 
(A'" = 4.5) of 9%. For beauty, the relative uncertainty A(x)°p/(x)°p = 1.4% translates into an uncertainty 
on large-pT .B-hadron production at LHC (N = 3.8) of 5.3%. These uncertainty are smaller or of the order 
of the perturbative uncertainties of the calculation. Nevertheless, it should be noted that this approach is only 
valid for large transverse momenta. At small transverse momenta the factorization ansatz breaks down and 
large corrections of order thq/pt may appear. Therefore, for the low -pT region, the uncertainty on the px 
distribution is large and difficult to evaluate. 
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x*^ / 


e(min : max) 


Q;(min : 


nicLx) 


\**^ // / 


Act/ct 
(LHC) 


Act/ct 
(HERA) 


FONLL D* 


0.849 ± 0.018 


0.0040(0.0027 : 


0.0057) 


10(9 : 




2.1% 


9% 


12% 


FONLL S""* 


0.927 ±0.013 


0.00045(0.00026 


; 0.00072) 


24(20 


: 30) 


1.4% 


5% 


7% 


FO-NLO D* 


0.75 ±0.10 


0.02(0.004 : 


0.08) 


5(3 : 


10) 


13% 


60% 


70% 


FO-NLO B""^ 


0.89 ± 0.04 


0.0015(0.0004 


: 0.004) 


15(10 


: 25) 


4.5% 


20% 


22% 



Table 4: Proposed value and uncertainty on {x)"p to be used with FO-NLO and FONLL programs for photo- and hadro-production 
of D* mesons and weakly decaying B hadrons. The corresponding value and range for the Peterson e and for the KartveUshviU a 
parameters are also reported. The last colimms show the corresponding relative uncertainty on da/dpr at LHC (assuming a negative 
power N = 4.5/3.8 for charm/beauty) and HERA {N = 5.5/5.0 for c/b). 

For processes such as DIS and for particular observables FONLL calculations are not available. The 
best theory available in this case is the fixed order NLO theory. In this case the situation is complex since the 
equivalent FO calculation for e+e~ does not reproduce the experimental data. The proposed solution is to vary 
(x)^^ from the same value obtained in the NLL case (that would be correct at low pT, where the FF evolution is 
irrelevant) to the value obtained at the Z° energy (that would be valid at px ~ 100 GeV). Therefore we consider 
for charm (a;)'^P = 0.075 ± 0.010 and for beauty (a;)''P = 0.089 ± 0.004. When these values are transported 
to heavy-hadron production at LHC, the corresponding uncertainties on da/dpx at large px are 60% for charm 
and 20% for beauty. Therefore the NLO fixed order calculations cannot be used for precise predictions of the 
charm (and to a lesser extent beauty) production at pp and ep colliders. 

In the FO-NLO and FONLL programs the hadron distributions are obtained by reducing the quark mo- 
menta according to a given fragmentation functions. Typical fragmentation functions used in these programs 
are the Peterson and Kartvelishvili forms. Table 4 summarises the proposed values and uncertainties for 
to be used with FO-NLO and FONLL calculations and reports the corresponding values and ranges for the 
Peterson and Kartvelishvili parameters. Similar ranges are used in the calculations presented in the section on 
"Benchmark cross sections" in these proceedings. 
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3 Measurements of the charm quark fragmentation function at HERA 

authors: J. Bracimk and G. Grindhammer 

The differential cross section for the inclusive production of a heavy hadron H from a heavy quark h 
can be computed in perturbative QCD (pQCD) as a convolution of a short-distance cross section a{ph) with a 
fragmentation function D^{z): 



The quantity z is the fractional momentum of the heavy quark h which is transferred to the heavy hadron 
H, and the normalized fragmentation function D^{z) gives the probability to observe the hadron H with a 
momentum fraction z. 

The precise definition of {z) is in some sense arbitrary. Due to the short and long-distance processes 
involved, the fragmentation function contains a perturbative and a non-perturbative component. Since the 
former can be calculated only up to some order in the strong coupling, the non-perturbative component in 
practice will have to absorb some of the missing higher order corrections. The calculable perturbative part can 
be absorbed into the definition of a{ph). Since for heavy quarks perturbative gluon emission do not lead to 
coUinear divergencies, the perturbative evolution is well defined, and it is possible to absorb them into a{ph) 
and to perform perturbative evolution down to a scale of the heavy quark mass rrih- In this case the non- 
perturbative fragmentation function D^{z) accounts for the transition of an almost on-shell quark h into a 
heavy hadron H. 

According to the QCD factorization theorem, the non-perturbative fragmentation functions (FF) depend 
neither on the type of the hard process nor on the scale at which the heavy quark h is originally produced. This 
implies universality of FF and allows - if valid - to extract fragmentation functions from data for one particular 
reaction (usually e+e^ annihilation) and to use them to predict cross sections for other reactions (e.g. in pp 
and ep-coUisions). In order to be able to check the reliability of pQCD predictions, it is necessary to check the 
universality of FF. 

In practice, different theoretically motivated functional forms for D^{z) are used, depending on one 
more free parameters which are fitted to data. Among frequently used expressions are those by Peterson et 
al. [14] and by Kartvehshvih et al. [12]. 

From Equation 1 1 it is clear that D^{z) cannot be measured directly, since all observables are convoluted 
with the perturbative cross section. In case of ep and pp scattering there are additional convolutions with the 
parton density functions of one or two interacting hadrons. However, there are some observables which are 
more sensitive to D^{z) then others. 

In e^e~ , a convenient way to study fragmentation is to study the differential cross section of a heavy 
meson as a function of a scaled momentum or energy z. A customary experimental definition^ of z is z = 

^Sometime there are slightly different definitions of z [28] in case of heavy meson production close to threshold. 




(11) 
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Fig. 6: Normalized differential cross section as a function of zjct as measured by ZEUS in photoproduction for jets with an associated 
D*-meson with \rjjct \ < 2.4 and -Etjci > 9 GeV. 

Eh / E\ycam^ where i^bcam is the energy of the beams in the center-of-mass system. In leading order, i.e. 
without gluon emissions, it is also the energy of the charm and anti-charm quark and is equal to D^{z). In 
contrast to e+e~ annihilation the choice of a fragmentation observable in ep collisions is more difficult. Two 
different observables have been used so far, both of them having the feature that in leading order QCD, the 
z-distributions are equal to D^{z). 

In the case of what is called here the jet method, the energy of the charm quark is approximated by the 
energy of the charm-jet, tagged by a L'*-meson, which is considered to be part of the jet. The scaling variable 
is then defined as zjet = {E + pi,)d* /{E + p)]ct- 

The idea of the so called hemisphere method (see Figure 5) is to exploit the special kinematics of charm 
events in ep collisions. The dominant charm production process has been shown to be boson-gluon fusion. If 
such an event is viewed in the photon-proton center-of-mass frame, the photon puts its full energy into the hard 
subprocess, while the proton interacts via a gluon, which typically carries only a small fraction of the proton 
momentum. As the result, both quarks produced, c and c, move in the direction of the photon. Assuming no 
initial gluon and no gluon radiation, their transverse momenta are balanced (see Fig. 5, left). 

This can be seen best by projecting the quark momenta onto the plane perpendicular to the 7-p axis. In 
this plane it is possible to distinguish rather efficiently between the products of the fragmentation of the charm 
quark and its antiquark. The momenta of all particles are projected onto the plane and the thrust axis in this 
plane is found (see Fig. 5, right). The plane is then divided into two hemispheres by the line perpendicular to the 
thrust axis. All particles, lying in the hemisphere containing the L)*-meson are marked and their three-momenta 
and energy are summed-up to give the hemisphere's momentum and energy, which is used to approximate the 
momentum and energy of the respective charm/anti-charm quark. The scaling variable Zhem is then defined as 

Zhem = (£^+PL)D'/Ehcm(^+P)- 

The ZEUS collaboration has provided preliminary results [44] on a measurement of normalized differ- 
ential cross sections of L'*-mesons as a function of zjot- The measurement was done in photoproduction, in 
the kinematic range 

Q2 < 1 GeV^ and 130 < < 280 GeV. The D* -mesons were reconstructed using the 
'golden channel' D* —>■ D^tTs —>■ KtttTs and were required to be in the central rapidity region \r]\ < 1.5 and 
to have > 2 GeV. Jets were reconstructed using the inclusive k± algorithm. They fulfill the conditions 
|77jct| < 2.4 and -ETjet > 9 GeV. The jets were reconstructed as massless jets. The beauty contribution to the 
L>*-meson cross section, which amounts to about 9%, was subtracted using the prediction of PYTHIA. The 
scaling variable was calculated as zjet = {E + pl)d* /(2-Ejet)- The cross section as a function of zjet is shown 
in Fig. 6. The uncertainties due to choice of the model used to correct for detector effects, and the subtraction 
of the beauty component were the largest contributions to the total uncertainty. 

The HI collaboration has recently presented preliminary results [45] on the normalized differential cross 
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Fig. 7: Normalized differential cross section of _D*-meson as a function of zjot and Zhcm in DIS as measured by HI. 



section also of D*-mesons as a function of both Zhom and zjet- Their measurement was performed in the 
kinematic range 2 < < 100 GeV^ and 0.05 < y < 0.7. The D* -mesons were reconstructed using the 
'golden channel' with \ri\ < 1.5 in the central rapidity region and px > 1-5 GeV. The jets were reconstructed 
using the inclusive k± algorithm in the photon-proton center of mass frame, using the massive recombination 
scheme. The jets were required to have i^'xjet > 3 GeV. The scaling variables were calculated as zjet = {E + 
Pl)d*/{E + p)jet and Zhem = {E + pl)d*/ Z]hem(-^ + P) ^"^^ shown in Fig. 7. The resolved contribution 
was varied between 10 and 50% and the beauty contribution as predicted by the model was varied by a factor of 
two. The resulting uncertainties are part of the systematic error of the data points. For these distributions, the 
contribution of D*-mesons coming from the fragmentation of beauty, as predicted by RAPGAP, was subtracted. 
It amounts to about 1.3% for the hemisphere method and 1.8% for the jet method. The dominant systematic 
errors are due to the model uncertainty and the signal extraction procedure. 



Parametrization 




Hemisphere 
Method 


Jet 
Method 


Suggested 
range 


Peterson 


e 


0.0i»_o.004 


U.UOU_Q QQ5 


0.014 < e < 0.036 


Kartvelishvili 


a 


-3-^-0.6 


^•^-0.5 


4 < a < 6.8 



Table 5: Extracted fragmentation parameters for zjct and Zhcm from HI. 

Both collaborations used the normalized z-distributions to extract the best fragmentation parameters for 
a given QCD model. 

In case of ZEUS, PYTHIA was used together with the Peterson fragmentation function. The MC was 
fit to the data using a -minimization procedure to determine the best value of e. The result of the fit is 
e = 0.064 ±0.006l[!:j]^^. 

The HI collaboration used RAPGAP 3.1 interfaced with PYTHIA 6.2. The contribution due to D*- 
mesons produced in resolved photon processes (in DIS), which amounts to 33% as predicted by the model, 
has been included in addition to the dominant direct photon contribution. The Peterson and Kartvelishvili 
parametrizations were both fitted to the data. The results are shown in Table 5. 

The parameter of the Peterson fragmentation function as measured by ZEUS and H 1 do not agree with 
each other. This may be due to the different phase-space regions covered by the two measurements (photo- 
production versus DIS, i^Tjet > 9 GeV versus -^Tjet > 3 GeV ) and most importantly, the parameters were 
extracted for two different models'^. More detailed investigations are needed to resolve this question. 

''while ZEUS has used the default parameters for PYTHIA, HI has taken the tuned parameter values of the ALEPH coUabora- 
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The fragmentation function parameters extracted by HI with the hemisphere and the jet method differ 
by less than 3 a. At the present level of statistical and systematic errors it is not possible to exclude a statistical 
fluctuation. On the other hand, the potential discrepancy may be a sign of deficiencies in the modelling of the 
hadronic final state in RAPGAP. 




Fig. 8: Comparison of the z-distributions from CLEO, OPAL and ALEPH (left) and Belle and ZEUS (right) with the one from the 
hemisphere method from HI. All distributions are normalized to unit area from z = 0.4toz = l. 

The measured Zhem distribution of HI is compared to data from the ALEPH [23], OPAL [22] and 
CLEO [28] collaborations in Fig. 8 (left) and to ZEUS [44] and Belle [29] in Figure 8 (right)^ The results 
of HI are in rough agreement with recent data from CLEO and Belle, taken at at 10.5 and 10.6 GeV, cor- 
responding roughly to the average energy of the system at HI. Differences beyond the measurement errors 
can be observed. However, this may be due to the somewhat different definitions used for the fragmentation 
observable z, different kinematics, different processes, or it may be a sign of the violation of universality. 



PARAMETRIZATION 


ALEPH 


OPAL 


BELLE 


ZEUS 


HI: Zhem 


Hi: Zjot 


Peterson e 


0.034 ± 0.0037 


0.034 ± 0.009 


0.054 


064+" "^^ 


U.0i8_o 


Q 030+0. UUti 


Kartvelishvili a 




4.2 ±0.6 


5.6 






^■3-0. 5 



Table 6: Extracted fragmentation parameters from e+e annihilation data by ALEPH [23], OPAL [22] and BELLE [29] and from ep 
data by ZEUS [44] and HI [45]. 

While the z distributions don't need to agree, the fragmentation parameters, which are extracted from 
them, should agree. This can be expected only, if a model with consistent parameter settings is used which 
provides an equally good description of the different processes at their respective scales. The values of the 
Peterson fragmentation parameter, as extracted by different experiments within the PYTHIA/JETSET models, 
are summarized in Table 6. Contrary to expectations, discrepancies between various experiments can be seen. 
A consistent phenomenological analysis of these data is therefore needed in order to resolve the reasons for the 
discrepancies. 

The measurement of the charm fragmentation function at HERA provides an important test of our un- 
derstanding of heavy quark production. We may hope that HERA II data and a phenomenological analysis of 
existing data will bring new insights in this area. 

tion [46] 

^Data points were taken from the figure in [29] and [44]. 
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Abstract 

Reference heavy-flavour cross sections at HERA and LHC have been computed fol- 
lowing different theoretical approaches and the results have been compared. 

coordinators: M. Corradi, A. Dainese 

1 Introduction 

This section presents a comparison of cross sections for HERA and LHC calculated according to different 
theoretical approaches. Different programs were used to calculate the same reference cross sections, using, 
as far as possible, the same input parameters and a consistent method to evaluate uncertainties. In this way it 
is possible to identify processes and kinematical regions in which different approaches give the same answer 
and regions where they differ. Unified criteria to evaluate the theoretical uncertainty should also allow to 
understand what approach is expected to be more precise. Moreover these calculations, which incorporate 
up-to-date parameters and PDF parametrisations, can be used as a reference for experiments and for further 
theoretical predictions. The cross sections presented here, are available in computer-readable format from the 
web page http : //www-zeus . desy . de/~ corradi /benchmarks, where figures in color can also be 
found. 

2 Programs 

A fist of the programs used for the cross section calculations is given below. For further details see the references 
and the theoretical review on heavy quark production in these proceedings. 

• MNR [1] is a fixed-order (FO) NLO program for heavy-flavour hadro-production, it was used for LHC 

cross sections; 

• FMNR [2, 3] is an extension of the previous program to photoproduction, it was used for photoproduction 
at HERA; 

• HVQDIS [4, 5] is a FO-NLO program for heavy-flavour production in deep-inelastic scattering (DIS), it 
has been used for DIS at HERA; 

• FONLL [6, 7] provides matched massive-massless calculations with NLO accuracy and resummation of 
large pt logarithms. It is available for hadro- and photo-production and was used for HERA photopro- 
duction and LHC cross sections; 

• GM-VFNS [8, 9, 10, 11] is a calculation in the generalised massive variable flavour number scheme. It 
has been used for charmed hadron px spectra at LHC and in photoproduction at HERA; 

• Cascade 1.2009 [12] is a full Monte Carlo program based on unintegrated parton densities and Kt 
factorisation. It has been used to calculate cross sections for Photoproduction and DIS at HERA and for 
LHC; 
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Parameter 


program 


central value 


lower/upper 


— T~[r 

■'^^QCD 


all 


0.226 GeV 


fix 


rUc 


all 


1.5 GeV 


1.3/1.7 GeV 


rub 


all 


4.75 GeV 


4.5/5.0 GeV 


p-PDF 


all-CASCADE 


CTEQ6.1 [15] 


MRST2002 [16]/Alekhin [17] 




Cascade 


CCFM AO 


- 


7-PDF 


FMNR, FONLL 


AGF [18] 


GRV [19] 




all 


0.235 


fix 


c fragmentation: 


(F)MNR,HVQDIS 


Pet. [20] Cc = 0.021 


0.002/0.11 




FONLL 


BCFYr = 0.1 


0.06/0.135 




GM-VFNS 


[9] 


- 




Cascade, Rapgap 


Pet. €c = 0.075 


Def./ec = 0.05 


b fragmentation: 


(F)MNR,HVQDIS 


Pet. Cb = 0.001 


0.0002/0.004 




FONLL 


Kart. a = 29.1 


25.6/34.0 




Cascade, Rapgap 


Pet. Cb = 0.002 


Def./cb = 0.005 



Table 1: The table shows input parameter to the different programs with the corresponding lower and upper values used for the 
uncertainty: Aqcd, the quark masses, the proton and photon parton densities, the fraction of c quarks decaying into a D* meson, and 
the parameters used for fragmentation. The fragmentation form are abbreviated to Pet. for Peterson, Kart. for Kartvelishvili, Def for 
the default Pythia fragmentation 

• Rapgap 3 [13] is a multi-puipose MC program for ep collisions, it implements heavy-flavour production 
through the boson-gluon-fusion process 7*5 QQ at leading order. It has been used for DIS at HERA. 
Both Cascade and Rapgap use Pythia [14] routines for fragmentation. 

3 Parameters and uncertainties 

The different calculations were compared using the same input parameters and, where possible, with total 
uncertainty bands computed in a consistent way. The total uncertainty band includes the effect of the uncertainty 
on the input parameters and on the missing higher orders in the perturbative expansion. 

3.1 Perturbative uncertainty 

The perturbative uncertainty was obtained by varying the renormalisation and factorisation scales independently 
in the range 0.5/xo < A^F, /^i? < 2;Uo, while keeping 1/2 < hr/ ijlf < 2, were /io is the nominal value, typically 
set to the transverse mass p\ + niQ or to 4m^ + in the DIS case. The largest positive and negative variations 
were taken as the perturbative uncertainty band. 

3.2 Input parameters 

The uncertainty from the input parameters was obtained by varying each parameter around the central value. 
An effort was made within the working group to find the best central value and uncertainty for the input 
parameters. The values used for the perturbative parameters Aq^j-,, rric, mj, as well as the parton distribution 
functions (PDF) for the proton and for the photon are reported in Table 1 . 

For practical reasons, rather than using the full treatment of the PDF uncertainty, few different parametri- 
sations were tried and it was checked that the choice of the PDF set always gives a small contribution to the 
total uncertainty band. In the case of Cascade, the CCFM AO parametrisation was used as the central value 
while the PDF parametrisations A0+ and AO-, obtained from fits to DIS data with different renormaUsation 
scales, were used in conjunction with the variation of the renormaUsation scale. 

Since the different programs have different perturbative contents, different parameters for the non- 
perturbative fragmentation function were used. The values were chosen in order to correspond to the same 
average fragmentation in e+e" colhsions as explained in the section on heavy quark fragmentation in these 
proceedings. Table 1 reports the fragmentation form and the corresponding parameters used in the different 
programs. 
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Fig. 1: Cross sections for charm photoproduction at HERA {Q^ < 1 GeV'^, 0.2 < y < 0.8). The differential cross sections as a 
function of the pr of the c quark for rapidity |y | < 2 and as a function of the rapidity of the c quark for pr > 2.5 GeV are shown 
in (a), (b). Plots (c) and (d) show similar cross sections for the production of a D* meson. The cross sections are shown for FMNR 
(shaded band), FONLL (empty band with dashed lines), GM-VFNS (empty band with dotted lines) and CASCADE (empty band with 
full lines). 



In the FONLL calculation for charm, the BCFY [21] fragmentation parameter r was varied in conjunction 
with the variation of the charm mass since different values of r are obtained from e+e~ data for different 
rric [22]. Similarly for beauty, the Kartvelishvili [23] parameter a was varied in conjunction with the variation 
of the b mass [23]. For GM-VFNS, the fragmentation functions and fractions were taken from [9]. 

The total uncertainty band was obtained from the sum of the uncertainties added in quadrature coming 
from the parameter variations and the perturbative uncertainty. 

4 Results 

4.1 HERA Photoproduction 

The results for HERA Photoproduction are given as ep cross-sections for 0.2 < y < 0.8 (y is the Bjorken 
variable while Y is the rapidity in the laboratory frame) and < I GeV^. The beam energies have been set 
to Eg = 27.52 GeV, Ep = 920 GeV with the proton beam going in the positive rapidity direction. 

Figure 1 shows the differential cross sections as a function of the charm quark transverse momentum 
(a) and pseudorapidity (b). In (c) and (d) the same cross sections are given for the charmed D* meson. A 
meaningful comparison can be performed only for the hadron variables, which are the real physical observables, 
since the quark level may be defined differently in different approaches. The FO calculation (FMNR) shows a 
large uncertainty (~ 60%) at the hadron level due to the related uncertainty on the fragmentation parameters. 
The resummed programs FONLL and GM-VFNS have much smaller uncertainty and are within the FMNR 
uncertainty band. The central values from FMNR and FONLL coincide at low transverse momenta. GM- 
VFNS, instead, tends to grow unphysically at low pt{D*). As can be seen in (c), the quark-level disagreement 
between FO (FMNR) and FONLL calculations is consistently removed at the hadron-level. The unintegrated- 
PDF Monte Carlo CASCADE tends to be above the other calculations, in particular at large px- In the case 
of beauty (Fig. 2) the uncertainty bands are smaller (~ 20% for FMNR), CASCADE and FMNR are in good 
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Fig. 2: Cross sections for beauty photoproduction at HERA (Q^ < 1 GeV^, 0.2 < y < 0.8). The differential cross sections in pr and 
rapidity of tiie b quark are shown in (a), (b). Plots (c) and (d) show the cross sections for the production of a weakly-decaying B hadron 
as a function of pt{B) and Y{B). The cross sections are shown for FMNR (shaded band) and CASCADE (empty band with full lines). 
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Fig. 3: Breakdown of the different components of the FMNR uncertainty for da/dpr for charmed (a) and beauty (b) hadrons in photo- 
production at HERA. The plots show the ratio of the upper/lower side of each uncertainty to the nominal value. The following sources 
of uncertainty are shown: quark mass (mg), parton density parametrisation (PDF), fragmentation parameter and the perturbatlve 
uncertainty from scale variations. 
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Fig. 4: Charm cross sections in DIS at HERA. Each plot shows the distribution of logjQ(a;) in a different range for the inclusive 
cross-section, the cross-section for a D* meson in the "visible" range pt{D*) > 1.5 GeV, |F(i?*)| < 1.5 and the cross-section for a 
muon from charm decay in the range prifJ^) > 3 GeV, \Y(iJ.)\ < 2. The thick curves show the central value from HVQDIS, the thin 
curves represent the uncertainty band from CASCADE and the shaded area shows the uncertainty band from Rapgap. The fluctuations 
in the muon cross sections are due to the limited statistics. 

agreement. Due to the large b mass, the resummed calculation FONLL (not shown) is expected to be similar 
to the fixed-order one (FMNR). For both beauty and charm, FMNR and FONLL show a shoulder at positive 
rapidities (b, d) due to the "hadron-like" component of the photon that is not present in CASCADE. 

Figure 3 shows the different components of the FMNR uncertainty band for charm and beauty. The un- 
certainties for quark production are typically dominated by the perturbative scale uncertainty with the exception 
of the low transverse momentum region (jjt rnq) where the uncertainty from the quark-mass can dominate. 
For hadron production, the fragmentation dominates the FMNR uncertainty at large px- The PDF uncertainty 
was found to be small. Resummed calculation have smaller uncertainty bands due to the smaller perturbative 
and fragmentation contributions at large px- 

4.2 HERA DIS 

Heavy quark production in DIS is not available in the matched massive-massless approach (except for total cross 
sections). Therefore the DIS comparison was limited to the FO-NLO program HVQDIS, the unintegrated-PDF 
MC Cascade and the Rapgap Monte Carlo. The DIS cross sections at HERA are reported as da / dlogiQ{x) 
for different bins of and are intended at the Born level, without electroweak corrections. Figure 4 shows, 
for each bin, the inclusive charm cross-section, the cross section for observing a D* meson in the "visible" 
range pt{D*) > 1.5 GeV, \Y{D*)\ < 1.5 and for observing a muon in the range prifJ-) > 3 GeV, |F(/u)| < 2. 
The three calculations are compatible at intermediate values of x (^ 10^^). At large x and low Q^, CASCADE 
and Rapgap drop to zero much faster than HVQDIS. At low x Rapgap is significantly larger than HVQDIS 
while both are within the uncertainty band given by CASCADE. A similar behavior is seen for beauty (Fig. 5). 
The uncertainty on HVQDIS, not given here, is expected to be small (~ 10 — 20% for beauty [24]). The 
high-x discrepancy between HVQDIS and the other two calculations seems therfore to be beyond the program 
uncertainties and deserves further investigations. 
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Fig. 5: Beauty cross sections in DIS at HERA. Each plot shows the distribution of logj^Q(a;) in a different range for the inclusive 
cross-section, the cross section for a hadron containing a b quark in the "visible" range Pt{B) > 3 GeV, jy(i3)j < 2 and the 
cross section for a muon from beauty decay in the range prifJ.) > 3 GeV, < 2. The thick curves show the central value 

from HVQDIS, the thin curves represent the uncertainty band from CASCADE and the shaded area shows the uncertainty band from 
Rapgap. The muon distributions are not given for HVQDIS. 
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Fig. 6: Cross sections for cliarm production in pp collisions at the LHC with ^/s — 14 TeV. The differential cross sections in pr for c 
quark in the two rapidity ranges \Y\ < 2.5 and 2.5 < |y | < 4 are shown in the upper panels. The lower panels show the cross sections 
for the production of a D* meson as a function of pt{D*) in the same rapidity ranges. 

4.3 LHC 

For LHC, we computed the cross sections in pp collisions at -y/i = 14 TeV. 

Figures 6 and 7 show the single inclusive cross sections as a function of pT, at quark (upper panels) and 
hadron (lower panels) level, for charm and beauty, respectively. Two rapidity intervals are considered: \Y\ < 
2.5, approximately covering the acceptance of the barrel detectors of ATLAS {\r]\ < 2.5), CMS (|?7| < 2.5), 
and ALICE (|?7| < 0.9); 2.5 < |y| < 4, approximately covering the acceptance of LHCb (2 < < 5) and of 
the ALICE muon spectrometer (2.5 < rj < 4). 

For charm, we compare the fixed-order NLO results from MNR to the results from the CASCADE event 
generator, from the GM-VFNS calculation and from the FONLL calculation. The agreement is in general 
good, in particular in the \ow-pT region; at high-p-r CASCADE predicts a larger cross section than the other 
calculations, especially at forward rapidities. The FONLL central prediction is in agreement with that of the 
FO NLO calculation at low px, while deviating from it at high px, where it gives a smaller cross section. 

For beauty, we compare FO NLO (MNR), FONLL and CASCADE. Again, there is agreement at low 
Pt, where, as expected, the FONLL result coincides with the MNR result. At high px, both CASCADE and 
FONLL predict a larger cross section than the MNR central values, but all models remain compatible within the 
theoretical uncertainties. At forward rapidities, for beauty as for charm, CASCADE gives a significantly larger 
cross section than MNR. 

Figure 8 shows the breakdown of the uncertainties for hadron production as obtained with MNR. The 
perturbative component dominates at LHC. Only the fragmentation component for charm hadron production 
becomes comparable in size to the perturbative one at large px- 

4.4 Q-Q correlations 

The azimuthal separation between the two heavy quarks A(j){QQ) and the transverse momentum of the quark- 
antiquark system pt{QQ) are particularly sensitive to higher-order effects since at leading order their distri- 
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Fig. 7: Cross sections for beauty production in pp collisions at the LHC with ^/s = 14 TeV. The differential cross sections in pr for b 
quark in the two rapidity ranges |y | < 2.5 and 2.5 < |y | < 4 are shown in the upper panels. The lower panels show the cross sections 
for the production of a beauty hadron as a function of pr in the same rapidity ranges. 




Fig. 8: Breakdown of the different components of the uncertainty on da/dpr for charmed (a) and beauty (b) hadrons at LHC as 
obtained from MNR. The plots show the ratio of the upper/lower side of each uncertainty to the nominal value. The following sources 
of uncertainty are shown: quark mass (mg), parton density parametrisation (PDF), fragmentation parameter and the perturbative 
uncertainty from scale variations. 
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Fig. 9: Q-Q correlations for charm at LHC: pr of the cc pair (left) and azimuthal angle between the c and the c (right). For both 
cross sections, the following kinematic cuts are applied: \Y'^\ < 2.5, |y^| < 2.5, > 3 GeV, > 6 GeV. 




bB [GeV] bb [rad] 

Fig. 10: Q-Q correlations for charm at LHC: pr of the bb pair (left) and azimuthal angle A(j) between the 6 and the 6 (right). For both 
cross sections, the following kinematic cuts are applied: \Y''\ < 2.5, \Y''\ < 2.5, > 3 GeV, Pt > & GeV. 

butions are delta functions peaked at /S.(f){QQ) = it and Pt{QQ) = 0. The distribution of these variables is 
therefore a direct probe of QCD radiation and is well suited for comparing different calculations. 

Figures 9 and 10 show the heavy-quark pair px distribution and the quark-antiquark relative azimuthal 
angle distribution for charm and beauty at LHC, respectively. For both distributions, the two quarks of the 
pair are required to have |y| < 2.5; also minimum px selections are applied to mimic the effect of realistic 

experimental cuts (pj, > 3 GeV and > 6 GeV). In the region near A(f){QQ) = vr and pt{QQ) = 0, where 
the cancellation of soft and coUinear divergencies occur, the fixed-order NLO calculation gives an unphysical 
negative cross section with next to a large positive peak. A larger binning would be needed to average this 
behavior and produce a more physical results. The CASCADE MC, has a more realistic behavior. Both calcu- 
lations have a non-zero value at A(I}{QQ) = related to "gluon-splitting" events. A similar result was found 
for HERA as shown in Figure 1 1. This kind of distribution is expected to be well described by programs that 
merge NLO matrix elements to the parton-shower MC approach such as MC@NLO [25]. 

5 Conclusions 

Heavy-flavour cross sections for HERA and LHC, obtained with fixed-order NLO programs, with matched 
massive/massless calculations and within the i^i^-factorisation approach have been compared. Similar results 
are found for photoproduction at HERA and for the LHC. As expected the resummed calculations were found 
to be compatible with the fixed-order results but have smaller uncertainties at large pT- Resummed calculations 
for charm in two different schemes (GM-VFNS and FONLL) are anyway somewhat incompatible both at 
HERA and LHC, suggesting that their uncertainty may be underestimated. The i^Tr-factorisation program 
Cascade predicts larger cross sections than the other approaches at large px at LHC and for charm at HERA. 
The comparison for DIS was limited to FO-NLO and a MC program with leading order matrix elements. Large 
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Fig. 11: Azimuthal Q-Q correlations in photoproduction at HERA for charm (a) and beauty (b). One of the two quarks was required to 
be in the "visible" region pt{Q) > 2.5 GeV and \Y{Q)\ < 2. 

discrepancies, which deserve further investigations, were found in this case. A comparison with experimental 
data would be needed for further understanding of the quality of the available calculations. 
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